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PART  A  -  DIRECTOR’S  OVERVIEW 


Beginning  May  1,  1987,  a  new  JSEP  three-year  program  was  initiated  at  Berkeley.  Two  major  themes 
were  identified:  High  Speed  Wide  Band  Elements  for  High  Frequency  Electronics  and  New  Architec¬ 
tures  for  Parallel  Computation.  The  first  contains  four  work  units  and  a  supplement  and  the  second, 
three  work  units. 

In  an  effort  to  coordinate  the  separate  parts  of  the  program  and  to  encourage  collaboration  between  the 
various  investigators,  a  couple  of  internal  reviews  have  been  held.  On  September  21,  1987,  there  was  a 
one-day  conference,  with  presentations  by  14  Research  Assistants  then  supported  by  the  program.  It 
was  decided  after  that  meeting  that  a  specialists’  meeting  within  the  two  team  areas  would  be  fruitful. 
On  February  5,  1988,  a  half-day  meeting  was  held  to  cover  all  the  activity  under  new  architectures  for 
parallel  computation.  A  similar  review  of  the  high-frequency  thrust  will  be  held  in  the  near  future.  We 
find  these  meetings  useful  and  essential  for  managing  the  nudges  in  the  direction  of  the  program. 

At  present  there  are  three  or  four  more  Research  Assistants  supported  by  this  program  than  can  be  con¬ 
tinued  with  the  combination  of  inflation  and  reduction  in  funding.  Thus,  it  is  critical  to  identify  the 
most  fruitful  areas  for  emphasis.  The  joint  service  electronics  program  continues  to  be  an  important 
factor  in  the  ability  of  this  laboratory  to  respond  to  important  new  research  opportunities.  The  partial 
funding  of  the  supplement  on  Electro  Optic  probing  and  the  GaAs  devices  and  materials  research  of  our 
new  faculty  member.  Professor  Steve  Smith,  are  positive  examples. 

In  the  following  we  list  a  number  of  significant  accomplishments  which,  because  of  the  recent  start  of 
the  new  three-year  program,  include  research  begun  under  the  three-year  proposal  ending  April  30, 
1987.  Finally,  in  Part  C,  we  describe  the  work  under  the  individual  work  units  and  cite  the  publications 
under  JSEP  sponsorship. 


PART  B  -  ACCOMPLISHMENTS 


High  Spe«d/Wide  Band  Silicon  and  GaAs  Devices 
Professor  Chenming  Hu  and  Professor  Ping  Ko 

As  part  of  a  project  to  construct,  model  and  compare  very  high  performance  Si  and  GaAs  devices, 
deep-submicron  Si  MOSFETs  have  been  investigated. 

Using  a  novel  optical  lithography  process,  these  extremely  high-gain  MOSFETs  with  varying 
oxide  thicknesses  and  channel  lengths  have  been  fabricated.  Figure  1  shows  a  Scanning  Electron 
Microscope  (SEM)  picture  of  a  transistor  cross-section  with  an  effective  channel  length  0.2  Uum,  an 
oxide  thickness  of  86  A,  and  junction  depth  of  0.2  um.  Devices  with  effective  channel  lengths  as  small 
as  0.15  um  and  oxide  thickness  as  thin  as  36  A  have  been  successfully  fabricated  and  operated. 


Figure  2  displays  a  plot  of  linear  transconductance  (a  measure  of  device  speed)  versus  the 
effective  channel  length.  MOSFETs  with  36  A  gate  oxides  and  0.15  um  effective  channel  lengths  exhi¬ 
bit  record  room-temperature  transconductances  of  460  mS/mm.  For  these  devices  intrinsic  cut-off  fre¬ 
quency  is  73  GHz.  NMOS  ring  oscillators  have  been  fabricated  with  delays  as  small  as  22  ps/stage. 
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Figure  2.  Plot  of  linear  transconductance  versus  the  effective  channel  lei 


Defect-Engineering: 


Proximity  Gettering  with  Mega-eiectron-volt  Carbon  Implantation 
Hing  Wong  and  Professor  Nathan  Cheung 

The  effectiveness  of  impunity  gettering  in  semiconductors  can  be  enhanced  by  placing  the  getter¬ 
ing  centers  close  to  the  active  device  region  while  maintaining  crystal  perfections  within  the  device 
region.  This  proximity  gettering  idea  is  illustrated  in  Figure  1.  With  Mega-electron-volt  (MeV)  ion 
implantation,  a  buried  gettering  layer  can  be  readily  formed  without  creating  significant  defects  in  the 
top  active  device  region.  Recently,  we  have  demonstrated  that  a  strong  proximity  gettering  layer  with 
excellent  thermal  stability  can  be  formed  with  MeV  carbon  implantation  in  silicon. 

MeV  carbon  ions  with  doses  ranging  from  5X1  (T14cm-2  to  5XKTt6cm-2  were  implanted  into  sili¬ 
con  wafers  contaminated  with  gold.  The  Au  was  intentionally  diffused  into  the  wafers  as  a  marker 
impurity.  The  wafers  were  then  annealed  in  an  N2  ambient  at  1000*0  for  one  hour.  Secondary  Ion 
Mass  Spectroscopy  (SIMS)  was  used  to  measure  the  depth  profiles  of  C  and  the  gettered  Au.  Shown  in 
Figure  2  (a)  is  the  SIMS  profile  of  2X10*16cnT2  2.4  MeV  carbon  implant  after  annealing.  A  large 
amount  of  Au  (5X10l3c/n-2)  was  trapped  in  the  C  implanted  region.  Complexing  of  the  C  with  the 
background  oxygen  is  suggested  by  the  O  peak  at  the  same  depth  as  the  carbon.  A  dose  dependence 
study  showed  that  the  gettering  efficiency  of  implanted  carbon  is  linear  with  dose.  The  thermal  stability 
of  the  carbon  gettering  effect  was  tested  for  annealing  up  to  12  hours  at  1000*  C.  The  amount  of  get¬ 
tered  Au  increased  with  annealing  time,  indicating  that  gettering  is  still  in  diffusion  limited  regime 
under  such  annealing  conditions. 

Cross-section  transmission  electron  microscopy  (XTEM)  was  used  to  study  the  structure  of  the 
carbon  implanted  gettering  layer.  Figure  2  (b)  is  an  XTEM  micrograph  obtained  on  the  same  wafer 
whose  SIMS  profiles  are  shown  in  Figure  2  (a).  Unlike  ions,  implanted  carbon  does  not  create  any  visi¬ 
ble  extended  defects,  such  as  dislocations  or  large  precipitates.  The  only  feature  in  the  carbon 
implanted  layer  is  a  vague  dark  band  with  a  thickness  of  about  0.5 um.  This  finding  suggests  that  the 
gettering  effect  of  carbon  is  through  point  defects  or  their  clusters,  not  through  dislocations  or  precipi¬ 
tate  interfaces. 

Four  MeV  carbon  ions  were  also  implanted  into  diode  structures  to  check  the  device  quality  of 
the  top  layer  after  MeV  implantation.  No  detrimental  degradation  in  diode  leakage  current  was  found. 
With  a  dose  of  1015cm"2  and  a  one-hour  annealing  at  950®  C,  the  leakage  increase  was  negligible  at  5 
volts  reverse  bias.  Even  at  a  dose  of  I0i6cm~2,  the  leakage  current  was  less  than  1  On  A/cm2  at  5  volts 
reverse  bias. 


Figure  Captions: 


figure  1.  The  effectiveness  of  gettering  can  be  enhanced  by 

putting  the  gettering  layer  close  to  the  active  device  region. 

figure  2.  (a)  SIMS  profiles  and  (b)  XTEM  micrograph  obtained  on 
a  sample  implanted  with  2X101<cm-2  2.4  MeV  carbon 
and  annealed  at  1000°  C  for  1  hour.  A  large  amount 
of  Au  (about  5X10l3cm-2)  was  trapped  in  the  carbon 
implanted  region.  No  extended  defects  were  observed. 

In  (b)  the  top  lum  of  silicon  was  sputtered  off  during 
TEM  sample  preparation. 
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GaAs  probing:  Surface  Properties  to  3-D  Field  Mapping 


•  V» 


kv*. 


v-'* 

rjv. 

<■ 

WO 

:S 


Yu-Hwa  Lo  and  Professor  S.  Wang 


Description 

Electro-optic  probing,  pioneered  and  demonstrated  first  in  LiNb03  by  Mouron  and  later  extended 
to  GaAs  by  Bloom,  has  the  potential  of  becoming  a  new  powerful  tool  for  materials  and  device 
research.  We  are  exploring  the  use  of  the  technique  to  study  materials  and  device  properties  which  are 
either  inaccessible  to,  or  difficult  for,  other  probing  techniques.  We  have  carried  on  experiments  in 
three  hitherto  unexplored  areas:  (1)  probing  surface  properties  of  a  semiconductor,  measuring  stand¬ 
ing  wave  ratio  along  a  coplanar  waveguide,  and  (3)  three  dimensional  mapping  of  internal  fields  in  dev¬ 
ices. 


1.  Probing  Surface  Properties:  t 

P 

The  surface  property  plays  an  important  role  in  device  performance.  One  problem  often  encoun¬ 
tered  in  processing  semi-insulating  (SI)  GaAs  substrates  is  thermal  conversion  of  the  conductivity  type. 
Backgating  in  GaAs  integrated  circuits  has  been  attributed  to  uie  change  in  EL2  defect  center  density. 
However,  it  is  almost  impost’ ble  to  detect  the  change  in  surface  conductivity.  Figure  1  shows  (a)  the 
configuration  of  a  semi-insulating  sample  in  relation  to  the  probing  beam,  (b)  the  measured  d^V.^/dx2 
signal  for  an  unannealed  sample,  and  (c)  the  measured  d2V,_^/dx2  signal  for  an  annealed  (at  700°  C) 
sample,  where  V,_*  is  the  electro-optic  signal.  From  the  direction  of  the  signal  movement,  the  surface 
type  can  be  ascertained.  The  surface  has  changed  from  semi-insulating  in  (b)  to  slightly  p  type  in  (c). 

2.  Standing-wave-ratio  Measurement: 

For  mm- save  integrated  circuits,  it  is  important  to  be  able  to  measure  the  standing  wave  ration 
along  waveguides  to  yield  information  about  circuit  elements  which  we  design.  Figure  2  shows  (a)  the 
experimental  set  up,  and  (b)  the  electro-optic  signal  observed  along  a  coplanar  waveguide.  The  E-O 
signal  shows  clearly  a  standing-wave  pattern.  The  experiment  was  first  carried  out  at  8.21  GHz  and  is 
now  extended  to  20  GHz.  The  electro-optic  measurement  of  SWR  should  complement  the  S-parameter 
measurement  in  characterizing  mm-save  integrated-circuit  elements. 
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3.  3-D  Field  Mapping: 


To  understand  the  physics  of  complicated  HI-V  compound  devices  and  to  set  up  physical  models 
for  prediction  of  device  performance,  it  is  essential  to  know  the  internal  potential,  field  and  charge  dis¬ 
tribution.  There  is  no  technique  at  present  available  to  provide  such  crucial  information  in  either  bulk 
semiconductors  or  devices.  Our  work  is  to  establish  the  feasibility  of  providing  the  information  by 
electro-optic  probing.  The  electro-optic  signal  senses  the  field  component  in  the  direction  of  light  pro¬ 
pagation.  By  changing  the  angle  of  incidence,  information  is  obtained  about  all  three  field  components. 
An  example  is  shown  in  Figure  3  for  a  simple  structure  consisting  of  two  electrodes  on  a  GaAs  surface. 
The  electro-optic  signal  expected  may  be  simulated  as  shown  in  Figure  4.  The  more  difficult  deconvo¬ 
lution  problem  (for  obtaining  the  orthogonal  fields  from  the  E-O  signal)  is  presently  under  investigation. 
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Miniaturized  of  Josephson  Logic  Circuits 
Herbert  Ko  and  Professor  T.  VanDuzer 

As  part  of  the  study  of  the  impact  of  high-Te ,  high-energy-gap  superconductors  on  digital  circuit 
performance,  the  effect  of  scaling  on  the  performance  of  Josephson  logic  circuits  to  smaller  sizes  has 
been  studied.  Theoretical  work  has  shown  the  importance  of  the  ratio  of  the  maximum  zero-voltage 
current  of  a  Josephson  junction  to  to  its  capacitance  as  the  determining  factor  in  setting  the  speed.  This 
ratio  is  set  in  tunnel  junctions  by  the  tunneling  current  density  (which  has  a  maximum  practical  value 
of  about  10  A/cm2  in  niobium  and  lead-alloy  unctions).  Detailed  simulations  of  two  logic  families  were 
carried  out.  The  minimum  cell  OR  gate  delay  is  about  6  ps  and  the  cell  delay  is  about  20  ps  in  both 
cases,  including  the  delay  of  propagation  to  the  next  gate. 

Calculations  also  were  made  on  a  model  of  a  logic  chip  in  which  the  chip  is  divided  into  cells, 
each  of  which  contains  two  OR  gates  driving  an  AND  gate.  Each  cell  also  has  space  for  ten  east-west 
transmission  lines  and  ten  running  north-south.  The  optimization  of  speed  and  circuit  density  can  be 
performed,  taking  into  account  practical  constraints  on  the  parameters  such  as  maximum  tunneling, 
current  density,  minimum  junction  current  (noise  limitation),  and  minimum  practical  dielectric  thickness 
and  permittivity.  The  optimum  junction  size  is  about  1  urn2.  Thus,  new  circuit  fabrication  processes 
were  developed  to  achieve  such  a  small  dimension.  A  technique  was  developed  that  provided  good 
control  of  the  junction  dimensions  and,  therefore,  tunneling  currents. 

Experiments  were  done  to  confirm  the  simulation  results  for  single  OR  and  AND  gates.  To  per¬ 
form  the  experiments,  an  on-chip  superconductor  sampler  circuit  was  fabricated  along  with  the  associ¬ 
ated  external  electronics  required  to  make  measurements  with  picosecond  resolution.  The  measure¬ 
ments  made  with  the  sampler  have  confirmed  the  scaling,  with  closest  agreement  on  the  faster  circuits. 
Previous  measurements  of  gate  speeds  have  been  done  on  chains  of  gates  so  the  results  were  averages 
and  the  bias  of  the  individual  gates  could  not  be  separately  controlled.  The  present  work  is  the  first  in 
which  the  delay  of  a  single  gate  was  studied. 
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Statistical  Design  of  Integrated  Circuits 


David  Riley  and  Professor  Alberto  Sangiovanni-Vincentelli 

Research  in  statistical  design  of  integrated  circuits  in  the  last  year  has  been  aimed  at  demonstrat¬ 
ing  the  feasibility  of  a  key  element  of  a  new  problem  formulation  for  statistical  design  that  we  have 
developed. 

Statistical  design  of  ICs  is  concerned  with  finding  optimal  circuit  parameter  values,  taking  into 
consideration  random  disturbances  inherent  in  wafer  fabrication  that  can  cause  the  performance  of  some 
of  the  circuits  produced  to  be  significantly  degraded. 

In  the  beginning  stage  of  the  project,  it  was  decided  that  several  basic  assumptions  made  in  con¬ 
ventional  formulations  for  statistical  IC  design  were  restrictive  when  applied  to  real  industrial  design, 
particularly  for  the  design  of  ICs  to  be  sold  on  the  open  market.  With  particular  attention  to  this  case, 
we  have  developed  a  new  problem  formulation  which  generalizes  conventional  formulations  in  a 
number  of  major  respects,  and  which  forms  the  basis  for  a  new  statistical  design  methodology.  Each 
item  of  the  following  list  compares  a  significant  restrictive  assumption  of  conventional  formulations 
with  the  related  relaxed  assumption  of  the  new  formulation. 

1.  The  number  of  electrical  performance  categories  of  the  circuits  produced  is  two 
(good/bad);  the  number  of  performance  categories  is  arbitrary. 

2.  The  design  parameters  optimally  determined  in  the  formulation  are  device  dimensions; 
they  include  all  parameters  having  a  first-order  effect  on  the  economic  value  of  the  circuits. 

3.  The  limits  defining  the  performance  categories  are  prescribed;  they  are  among  the  design- 
able  parameters  of  the  formulation. 

4.  In  the  most  common  problem  formulations,  parametric  yield  maximization,  the  criterion 
quantity  is,  in  equivalent  economic  terms,  the  ratio  of  revenue  to  cost;  it  is  a  form  of  profit  (revenue 
minus  cost)  appropriate  to  the  problem. 

5.  The  only  type  of  random  effect  degrading  the  economic  value  of  the  circuits  produced  is 
parametric  fluctuations  through  deviations  in  device  parameter  values;  defect  phenomena  can  also 
degrade  the  economic  value. 

Work  on  the  new  formulation  has  entailed  development  of  a  combined  technical  and  economic 
model  which  specified  the  following: 
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1.  The  criterion  quantity  (the  profit  to  be  expected  from  the  design),  the  required  constraint 
quantities,  the  design  parameters  and,  the  random  variables  required  in  the  statistical  modeling. 

2.  The  functional  dependence  of  the  first  two  in  Item  1  on  the  latter  two. 

Development  of  a  realistic  model  has  been  made  possible  through  extensive  collaboration  with 
Harris  Semiconductor.  The  model  has  been  documented  in  Ref  [1]. 

Note  that  the  design  parameters  include  fabrication  test  limits,  mask  dimensions,  back-end  fabri¬ 
cation  assignments  (e.g.  to  package  types  or  reliability  treatments),  final  product  performance 
specifications  and,  selling  prices. 

If  the  new  model  is  to  be  useful  in  industrial  practice,  it  must  be  feasible  to  obtain  all  the  data 
required  by  the  model.  Demonstration  of  this  feasibility  should  precede  development  of  the  software 
implementation  of  the  methodology. 

The  raw  cost  data  needed  for  the  cost  model  is  readily  available  to  the  producer.  However,  the 
revenue  model  requires  the  modeling  of  the  demand  function  of  the  IC  being  designed  -  -  the  expected 
quantity  demanded  as  a  function  of  the  products  specs  and  price.  The  demand  depends  on  the  behavior 
of  potential  users  in  choosing  among  commercially  available  realizations  of  the  type  of  the  IC  being 
designed.  Data  to  estimate  this  choice  behavior  is  not  readily  available  to  the  producer.  It  is  demon¬ 
strating  the  feasibility  of  estimating  the  demand  function  that  has  been  the  aim  of  the  work  in  the  last 
year. 

The  demand  model  is  based  on  the  postulation  of  an  appropriate  functional  form,  and  on  the  sta¬ 
tistical  estimation  of  parameters  of  the  function.  Data  for  this  estimation  is  gathered  in  a  survey  of 
engineers  throughout  American  industry  who  select  the  type  of  IC  of  interest  for  incorporation  into 
electronic  systems  they  design.  Additional  data  obtained  in  the  survey  allows  evaluation  of  alternative 
functional  forms.  The  model  has  been  developed  in  a  collaboration  with  the  Berkeley  Economics 
Department 

A  survey  aimed  at  estimating  the  demand  function  for  a  particular  type  of  IC  (linear  buffer 
amplifiers)  was  begun  in  September  of  last  year  to  demonstrate  the  feasibility  of  our  approach.  A 
mixed  phone  and  mail  survey  has  been  carried  out  with  the  assistance  of  part-time  student  employees. 
About  175  responses  have  been  obtained  to  date,  and  the  final  number  of  responses  is  expected  to 
exceed  200.  Statistical  estimation  of  parameters  for  the  demand  function  based  on  responses  obtained 
to  date  has  been  carried  out  successfully  and,  using  software  specially  developed  for  the  project,  prelim¬ 
inary  demand  estimation  computations  have  been  completed  and  appropriate  graphical  displays  pro¬ 
duced. 


Note  that  elements  of  our  demand  modeling  approach  have  application  to  two  very  significant, 
non-open- market  scenarios:  1)  the  procurement  of  an  1C  by  a  systems  design  division  of  a  firm  from  an 
IC  design  division  of  the  same  firm  and,  2)  the  procurement  of  an  IC  by  the  government  agencies 
through  competitive  bidding.  This  is  because  the  approach  includes  a  technique  for  modeling  the  per¬ 
formance  (and  price,  if  applicable),  and  tradeoff  judgments  of  individuals,  or  of  a  collection  of  individu¬ 
als,  in  an  electronic  system  design  group. 

Reference: 

(1]  D.C.  Riley  and  A.  Sangiovanni-Vincentelli,  "Models  for  a  New,  Profit-based  Methodology  for  Sta¬ 
tistical  Design  of  Integrated  Circuits",  IEEE  Transactions  on  Computer-Aided  Design  of  Integrated  Cir¬ 
cuits  and  Systems,  vol.  CAD-5,  No.  1.  pp.  131-169,  January  1986. 


Timing-Driven  Layout  for  Custom  Super  Chips 


Michael  A3.  Jackson,  Ren-Song  Tsay  and  Professor  E.S.  Kuh 

This  research  addresses  the  topic  of  timing-driven  layout  for  synchronous  digital  MOS  systems. 
Here  we  attempt  to  constructively  influence  the  timing  of  the  chip  by  manipulating  the  topology  of  the 
cell  placement  and  the  interconnect  routes.  Our  method  consists  of  two  main  states:  The  first  phase 
involves  a  pre-layout  timing  analysis  reminiscent  of  timing  verification,  a  "fast"  placement,  and  a  series 
of  optimization  problems  formulated  as  linear  programming  problems,  while  the  second  phase  requires 
layout  algorithms  that  use  the  results  of  the  optimization  step. 

The  pre-layout  timing  analysis  identifies  critical  paths  through  blocks  of  combinational  logic,  and 
the  slack  times  to  each  output  node.  In  parallel  with  this,  the  "fast"  replacement  yields  an  estimate  of 
interconnect  lengths  based  on  the  connectivity  of  the  design.  The  results  of  these  two  analyses  are  then 
formulated  as  linear  programming  problems.  Solutions  to  the  linear  programming  problems  yield  an 
interconnect  length  (spatial  constraint)  for  every  net  that  guarantees  the  desired  performance,  once  each 
spatial  constraint  is  satisfied  during  layout. 

The  second  phase  of  the  methods  consists  of  placement  and  routing  algorithms  that  accept  spatial 
constraints  and  generate  layouts  satisfying  them.  This  new  methodology  has  been  implemented  in  the 
routing  tools  of  BEAR,  a  new  building-block  layout  system  under  development  at  U.C.  Berkeley.  In 
the  event  that  all  constraints  are  not  satisfied,  sophisticated  interactive  techniques  arc  used  to  resolve 
unfulfilled  constraints.  Implementation  of  these  ideas  within  the  floor-planning/placement  algorithms  is 
currently  under  development  and  preliminary  experimental  results  with  the  routing  tools  are  very 
encouraging. 

In  a  separate  but  related  project,  we  have  developed  a  method  to  place  cells  for  superchips.  The 
method  is  based  on  the  solution  of  linear  sparse  equations  and  successive  partitioning.  It  is  extremely 
fast  and  the  quality  of  the  result  is  excellent.  A  triple-layer  100K  gate  array  with  26K  instances  was 
placed  in  50  minutes  using  VAX  S65Q.  Sufi  to  be  implemented  is  timing-driven  algorithm  described 
above. 
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PART  C  -  INDIVIDUAL  WORK  UNITS 

TH£ME  I  •  High-speed  Wide-band  Elements  for  High-frequency  Electronics 

HFD.l.  High  Speed/Frequencies  GaAs  and  Si  Device  Research 
C.  Hu 

The  current  resolution  limit  of  conventional  optical  lithography  has  required  that  deep-submicron 
MOS  devices  be  fabricated  exclusively  by  electron  beam  and  X-ray  lithography.  Both  techniques  are, 
however,  expensive  and  complicated.  We  have  developed  a  photoresist-thinning  process  that,  when 
used  in  conjunction  with  conventional  optical  lithography,  permits  the  controlled  definition  of  deep- 
submicron  features. 

N  channel  MOS  transistors  with  0.15-pm  effective  channel  length  and  showing  excellent  device 
characteristics  have  been  successfully  fabricated  using  conventional  g-line  optical  lithography  and  this 
thinning  technique. 

We  have  found  the  TV  characteristics,  gain,  and  short  channel  effects  in  these  0.15  to  0.5  pm 
MOSFET  are  well  described  by  the  models  we  have  developed  for  1  pm  devices.  The  basic  scaling 
guidelines  developed  for  the  1-pm  device  are  also  applicable  to  deep-submicron  devices.  However, 
lowering  the  power  supply  from  5V  to  3.3V  fails  to  solve  the  hot-electron  reliability  problems;  some 
form  of  hot-electron  resistant  structure  will  be  needed. 

We  have  discovered  that  when  a  thin-gate  oxide  MOSET  is  biased  in  the  off-state,  a  large  drain 
leakage  current  can  be  observed  between  drain  and  substrate.  Gate  current  is  also  observed  at  higher 
drain  voltages.  These  currents  are  important  to  short-channel  thin-oxide  devices  and  some  high-voltage 
devices,  such  as  EPROM  and  EEPROM  cells.  A  proposed  model  based  on  the  band-to-band  tunneling 
theory  can  explain  the  experimental  results  and  suggest  guidelines  for  the  design  of  thin-oxide  VLSI 
MOSFETs.  We  attribute  the  gate-induced  drain  leakage  current  to  the  band-to-band  tunneling  in  silicon 
in  the  gate/drain  overlap  region,  where  a  deep  depletion  region  is  formed.  Electrons  and  holes  are  gen¬ 
erated  by  the  tunneling  of  valence  band  electrons  into  the  conduction  band  and  collected  by  the  drain 
and  substrate  separately.  This  effect  sets  a  new  limit  on  submicron  device  scaling. 
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As  part  of  the  continuing  effort  to  formulate  GaAs  device  models  for  circuit  simulation,  we  are 
undertaking  improvements  to  our  models  and  examining  fresh  approaches.  The  principal  goals  are  to 
develop  a  robust  extraction  system  based  on  a  PC,  and  to  implement  the  model  in  SPICE3.  The  current 
model  displays  excellent  agreement  with  both  HEMT  and  MESFET  characteristics  and  provides  a  one- 
region  description  incorporating  threshold  shift,  subthreshold  slope  variation,  transconductance  degrada¬ 
tion,  and  breakdown.  It  also  includes  a  single-region  Schottky  diode  and  capacitance  descriptions. 

To  simulate  numerically  the  phenomenon  of  negative  conductivity  in  GaAs,  we  have  developed  a 
program  based  on  an  implicit  difference  scheme.  The  simulator  can  analyze  the  field  and  carrier  con¬ 
centration  behavior  in  a  device  for  a  variety  of  doping  profiles  and  boundary  conditions,  and  provides 
information  about  the  device  impedance,  microwave  power,  and  efficiency  through  a  Fourier  analysis  of 
resulting  current  and  voltage  waveforms.  The  program  has  been  tested  by  comparison  with  known 
results.  It  is  being  used  to  examine  the  feasibility  of  a  proposed  high-frequency  oscillator  based  on  the 
simultaneous  propagation  of  multiple  domains  through  GaAs. 

Publications 

T.  C.  Ong,  P.  K.  Ko,  and  C.  Hu,  "Modeling  of  substrate  current  in  P-MOSFETs,"  IEEE  Electron  Device 
Letters,  vol.  EDL-8,  no.  9,  September  1987,  pp.  413-416. 

J.  Chen,  T.  Y  Chan,  I.  C.  Chen,  P.  K.  Ko,  and  C.  Hu,  "Sub-breakdown  drain  leakage  current  in  MOS- 
FETs."  IEEE  Electron  Dev.  Lett.,  vol.  EDL-8,  no.  11,  November  1987. 

T.  Y.  Chan,  J.  Chen,  P.  K.  Ko,  and  C.  Hu,  "The  impact  of  gate-induced  drain  leakage  current  on  MOS- 
FET  scaling,"  Technical  Digest  of  International  Electron  Devices  Meeting  (IEDM),  Washington, 
EXT,  December  1987. 

S.  Holland,  I.  C.  Chen,  and  C.  Hu,  "Ultra-thin  silicon-dioxide  breakdown  characteristics  of  MOS  dev¬ 
ices  with  n+  and  p+  polysilicon  gates,"  IEEE  Electron  Device  Letters,  vol.  EDL-8,  no.  12, 
December  1987,  pp.  572-575. 

M.  C.  Jeng,  J.  Chang,  A.  T  .  Wu,  T.  Y.  Chan,  J.  Moon,  G.  May,  P.  K.  Ko,  and  C.  Hu,  "Performance 
and  hot-electron  reliability  of  deep-submicron  MOSFET’s,"  Technical  Digest  of  international 
Electron  Devices  Meeting  (IEDM),  Washington,  DC,  December  1987. 


HFD.2.  GaAs/Si  Devices  and  Materials  Studies 
S.  Wang 

Two  important  issueds  facing  GaAs/Si  heteroepitaxy  are  (1)  dislocation  density  and  (2)  internal 
stress.  We  have  designed  a  series  of  experiments,  including  selected-area  growth,  rapid-thermal  anneal¬ 
ing,  and  plasma  hydrogenation,  to  test  their  effects  on  dislocation  density  and  used  TEM  to  discern  the 
changes  in  threading  dislocations  and  stacking  faults.  We  have  also  performed  Raman  and  photo- 
luminescence  (PL)  experiments  to  determine  the  internal  stress  as  a  function  of  distance  from  the  inter¬ 
face  in  the  former  and  as  a  result  of  differential  thermal  expansion  in  the  latter. 

Since  the  results  on  Raman  and  PL  experiments  have  already  been  published  (see  publications  1- 
4,  they  are  discussed  first  A  bevel  was  etched  in  a  GaAs  epitaxial  film  grown  on  a  Si  substrate  so  that 
the  Raman  spectrum  of  the  GaAs  layers  can  be  measured  as  a  function  of  distance  from  the  GaAs/Si 
interface.  The  amount  of  strain  and  disorder  has  been  estimated  from  the  GaAs  longitudinal  optical  pho¬ 
non  line  shape  and  frequency  to  decrease  from  a  value  of  5.5xl(T3  at  a  distance  500  A  from  the  inter¬ 
face  to  a  value  around  lxlO-3  at  a  distance  2000  A  from  the  interface.  Therefore,  the  misfit  strain 
caused  by  lattice  mismatch  is  largely  relaxed  in  a  film  of  5000  A  thick  by  the  generation  of  edge  dislo¬ 
cations. 

Even  though  misfit  strain  is  relaxed,  we  still  have  to  contend  with  stress  caused  by  difference  in 
the  thermal  expansion  coefficients  of  GaAs  and  Si.  For  the  PL  experiments,  GaAs  films  were  grown  on 
Si  substrates  with  patterned  SiN  masks  by  molecular  beam  epitaxy.  The  pattern  consists  of  bare  Si 
stripes  of  width  ranging  from  10  pm  to  100  pm  surrounded  by  SiN  on  both  sides  and  a  reference  area 
of  bare  Si.  The  PL  experiments  were  carried  out  at  77K  to  yield  information  about  the  internal  stress 
and  the  optical  quality  of  the  hetero-epitaxial  film.  Under  stress,  the  degeneracy  of  the  valence  band  is 
removed,  and  the  heavy-hole  and  light-hole  gap-energies  shift  with  stress.  From  the  PL  peak  shift  and 
separation,  the  stress  can  be  determined  and  is  found  to  be  3.5xl09  dynes/cm2  at  77K  corresponding  to 
a  thermal  stress  of  2xl09  dynes/cm:  at  300K.  For  1.5  dm  and  3  pm  thick  films,  PL  intensity  from  the 
10  pm  stripes  showed  a  140%  and  75%  increase  over  the  reference  area,  respectively.  This  remarkable 
increase  in  the  PL  intensity  is  believed  due  to  a  reduction  of  dislocations  inside  the  window  area. 


TEM  pictures  of  selective-area  grown  samples  show  that  stacking  faults  are  confined  to  the  boun¬ 
dary  of  single  crystal  and  polycrystalline  regions.  TEM  pictures  of  the  rapid-thermal  annealed  (RTA) 
samples  show  a  significant  reduction  of  threading  dislocations  and  an  almost  elimination  of  stacking 
faults  in  the  single  crystal  (window)  region.  Work  is  in  progress  to  obtain  TEM  pictures  for  selective- 
area  grown  samples  after  RTA.  Preliminary  experiments  on  hydrogenated  samples  show  noticeable 
improvements  in  photo-response. 

In  summary,  we  have  tried  a  number  of  approaches  to  reduce  the  defect  density  in  and  to  improve 
the  optical  quality  of  hetero-epitaxial  GaAs/Si  films.  Selective-area  epitaxy  appears  effective  in  moving 
the  defects  to  the  area  near  the  single  crystal  and  polycrystalline  boundary.  Rapid  thermal  annealing 
also  seems  effective  in  defect  reduction.  Once  the  density  is  reduced  to  below  a  certain  level,  other 
techniques  such  as  doped  superlattices  may  become  more  effective  in  further  reducing  the  defect  den¬ 
sity.  We  believe  that  a  combined  use  of  selective -area  growth,  incorporation  of  doped  superiattices, 
post-growth  RTA,  and  post-growth  hydrogenation  can  bring  the  defect  density  to  a  level  suitable  for 
minority-carrier  devices. 

Publications 

Y.  H.  Huang,  P.  Y.  Yu,  M.-N.  Charasse,  Y.  H.  Lo,  and  S.  Wang,  “A  Raman  Study  of  Epitaxial  GaAs 
Layer  on  Si  (100)  Substrate,”  Appl.  Phys.  Lett.,  Vol.  51,  p.  192,  1987. 

Y.  H.  Lo,  M.-N.  Charasse,  H.  Lee,  D.  Vakhshoori,  Y.  Huang,  P.  Y.  Yu,  Z.  Liliental-Weber,  M.  Wemer, 
and  S.  Wang,  “Raman  Scattering  Spectrum  Along  a  Bevel-Etched  GaAs  on  Si,  TEM  Study,  and 
GaAs  PIN  Photodetector  on  SI,”  Materials  Research  Society  Symposium  Proceedings,  Vol.  91, 
Heteroepitaxy  on  Silicon  Technology,  p.  149,  1987. 

H.  P.  Lee,  S.  Wang,  Y.-H.  Huang,  and  P.  Yu,  “Photoluminescence  Studies  of  Selective-Area  Molecular 
Beam  Epitaxy  of  GaAs  Film  on  Si  Substrate,”  Appl.  Phys.  Lett.,  Vol.  52,  p.  215,  1988. 

Y.-H.  Huang,  P.  Y.  Yu,  H.  P.  Lee,  and  S.  Wang,  “Characterization  of  GaAs  Film  Grown  on  Si  Sub¬ 
strate  by  Photoluminescence  at  77K,”  Appl.  Phys.  Lett.,  Vol.  52,  p.  579,  1988. 


HFD.3.  Basic  Techniques  for  Electromagnetic  Scattering  and  Radiation: 

Transmission  Line  Systems  for  Millimeter/Submillimeter  Propagation 
D.  J.  Angelakos  and  K.  K.  Mei 

Copianar  transmission  lines  are  being  used  in  monolithic  submillimeter  devices.  Because  of  the 
geometry,  a  number  of  problems  arise  related  to  coupling  between  lines  and  sections  of  lines.  In  addi¬ 
tion,  some  topological  arrangements  between  balanced  and  unbalanced  lines  lead  to  problems  of  match¬ 
ing  and  radiation.  The  time  domain  finite  difference  method  described  above  is  very  suitable  for  such 
problem  as  well.  Some  preliminary  results  are  already  obtained.  A  full  investigation  for  the  coupling 
problem  will  require  the  access  of  a  Cray  2  computer. 

Pulse  Radiating  Antennas 
Jiayuan  Fang 
(Professor  K.  K.  Mei) 

There  are  many  types  of  broadband  antennas,  but  most  of  them  are  phase  dispersive.  So,  not  all 
broadband  antennas  can  radiate  short  pulses.  This  project  intends  to  identify  and  investigate  antennas 
which  have  minimal  phase  dispersion  so  that  they  are  high  fidelity  both  in  phase  and  amplitude. 

Computations  of  High-Frequency  Circuit  Components  in  Superconductivity  Devices 
Guo  Chun  Liang 

(Professors  K.  K.  Mei  and  T.  Van  Duzer) 

Superconducting  devices  involve  numerous  circuit  components  that  need  to  be  analyzed  and 
measured.  Both  computation  and  direct  measurement  are  important,  because  neither  approach  along  can 
provide  reliable  data.  Time-domain  techniques  will  be  used  to  calculate  passive  components  used  in 
superconducting  devices.  These  include  microstrips,  copianar  waveguides,  and  slot  lines.  Recently,  we 
also  find  that  it  may  be  easier  to  use  slow  wave  structure  in  the  measurement  and  computations. 
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HFD.4.  Submillimeter- Wave  Mixing  with  Superconductive  Tunnel  Junctions 
Michael  A.  Bruns 
(Professor  T.  Van  Duzer) 

The  goal  of  this  work  is  to  study  die  application  of  small  (=1  pm2)  superconductor-insulator- 
superconductor  (SIS)  tunnel  junctions  in  an  integrated  antenna-mixer  combination  for  submillimeter 
waves  at  frequencies  up  to  600  GHz.  Recent  research  at  UC  Berkeley  and  elsewhere  has  already 
demonstrated  that  single  SIS  junctions  and  junction  arrays  can  operate  as  efficient  heterodyne  mixers  for 
millimeter  waves  with  input  noise  approaching  the  quantum  limit  The  mixing  properties  of  the  SIS  tun¬ 
nel  junction  result  from  the  extremely  nonlinear  dc  current-voltage  characteristic,  the  major  feature  of 
which  is  a  sharp  rise  in  tunneling  current  at  the  gap  voltage.  Quantum  mixer  theory  predicts  very  low 
noise  and  even  conversion  gain  (an  impossibility  in  classical  mixer  theory)  if  the  current  rise  of  the  dc 
characteristic  is  sharp  on  a  voltage  range  of  hf  te ,  where  /  is  the  signal  frequency,  h  is  Planck’s  con¬ 
stant,  and  e  is  the  electron  charge. 

Another  important  factor  is  the  junction  shunt  capacitance,  which  should  be  large  enough  to  short 
out  undesired  higher  harmonics  of  the  local-oscillator  frequency,  but  not  so  large  that  the  wasted  dis¬ 
placement  current  constitutes  a  very  large  fraction  of  the  total  junction  current  As  the  local  oscillator 
frequency  increases,  the  junction  capacitance  must  decrease  to  maintain  good  mixer  performance.  This 
constraint  is  not  as  severe  when  using  series  arrays  of  tunnel  junctions,  because  the  overall  capacitance 
is  reduced  with  the  junction  capacitances  in  series.  Inductive  tuning  structures  may  also  be  used  to  can¬ 
cel  out  junction  capacitance. 

The  600-GHz  mixer  experiment  also  differs  from  similar  experiments  at  lower  frequencies  in  that 
the  signal  and  local  oscillator  are  not  coupled  and  guided  by  classical  waveguide  techniques.  Instead,  a 
quasi-optical  approach  is  used,  coupling  and  guiding  the  signal  with  mirrors  and  lenses  to  a  focus  on 
the  SIS  junction  itself. 

The  first  phase  of  the  quasi-optical  receiver  study  employes  90-GHz  RF  signals,  for  which  the 
design  and  testing  is  easier  and  more  accurate  than  for  the  600-GHz  system.  A  complete  test  system  has 
been  developed  in  the  Physics  Department  and  studies  are  in  progress  using  junctions  from  another 


source. 


The  first  part  of  the  research  on  this  project  has  been  directed  towards  the  fabrication  of  SIS  tun¬ 
nel  junctions  with  low  capacitance  and  dc  I-V  characteristics  with  a  sharp  current  rise  at  the  gap  vol¬ 
tage.  The  junctions  will  subsequently  be  incorporated  in  a  thin-film  antenna-mixer  with  integrated 
impedance -matching  structures.  It  is  expected  that  experiment  will  deviate  from  theory,  especially  for 
the  600-GHz  structures,  and  iterated  design  adjustments  will  be  required  for  the  thin-film  devices. 

We  have  devised  a  fabrication  procedure  for  making,  with  accurate  control,  1  pm2  Nb-A10x-Nb 
tunnel  junctions  using  a  crossed-strip  technique.  The  various  steps  are  being  developed. 


SPL.l.  Electro-Optic  Probing  of  Semiconductor  Materials  and  Devices 
S.  Wang 

The  work  on  electro-optic  probing  has  been  carried  out  in  two  fronts:  (1)  exploring  the  possibility 
for  three-dimensional  mapping  of  internal  field  and  (2)  demonstrating  the  capability  of  standing-wave 
measurements  in  coplanar  microwave  waveguides  (see  publications  1  and  2). 

The  electro-optic  probing  technique  is  for  the  first  time  proposed  to  detea  three-dimensional  field 
distribution  in  linear  electro-optic  material  like  GaAs.  By  changing  the  incident  angles  and  positions  of 
the  probing  beam,  sufficient  information  of  the  field  disribution  should  be  included  in  the  phase  retarda¬ 
tion  of  the  probing  beam.  We  have  shown  that  if  certain  conditions  on  the  probing  beam  are  saricficH  a 
very  simple  linear  relation  between  phase  retardation  and  each  field  component  can  be  found  A  three- 
dimensional  computer  simulation  was  undertaken  to  illustrate  the  relation  between  a  given  field  distribu¬ 
tion  and  the  predicted  electro-optic  signal.  A  simple  experimental  demonstrating  the  sensitivity  of  the 
electro-optic  signal  with  the  incidence  angle  was  also  performed.  These  results  show  that  three- 
dimensional  field  mapping  by  electro-optic  probing  is  not  only  mathematical  solvable  but  also  experi¬ 
mentally  feasible. 

We  have  also  successfully  measured,  for  the  first  time,  standing  waves  in  a  GaAs  coplanar 
waveguides  at  frequencies  of  8.2107  and  12.310  GHz  by  using  harmonic-mixing  electro-optic  probing 
technique.  The  technique  consists  of  mixing  the  microwave  signal  at/*  with  the  n-th  harmonic  nf0  of 
the  mode-locked  laser  pulse  of  repetition  frequency  f0  to  produce  an  electro-optic  probing  signal  at 
fm  ±  «/» •  This  technique  imporives  greatly  the  signal  to  noise  ratio,  making  possible  SWR  measure¬ 
ments  at  high  microwave  frequencies.  The  experiment  has  since  been  extended  to  20  GHz.  We  believe 
that  the  technique  can  provide  valuable  information  about  impedance  matching  in  the  design  of  mon¬ 
wave  circuits. 

In  addition  to  the  work  on  electro-optic  probing,  we  have  also  carried  out  exploratory  work  in 
several  areas:  (1)  resonant  tunneling  diodes,  (2)  second- harmonic  generation  in  GaAs  waveguides,  and 
(3)  lasing  action  in  GaAs/Si  films.  These  studies  are  undertaken  to  explore  potential  high-speed  elec¬ 
tronic  and  opto-electronic  devices  for  GaAs/Si  hetero-epitaxy.  Our  theoretical  study  shows  that  a  high 
peak  to  valley  ratio  can  be  achieved  in  resonant  tunneling  diodes  if  the  barrier  is  made  high  and  the 


barrier  layer  is  made  thin.  The  study  provides  guidance  for  an  optimal  design  of  the  device,  (see  publi¬ 
cation  3).  The  second-harmonic  experiment  originally  was  intended  to  detect  stress-induced  changes  in 
electro-optic  effect  During  the  investigation,  we  observed  a  coherent,  second-harmonic  emission  normal 
to  the  surface  (see  publication  4)  of  a  GaAs/AlGaAs  waveguide.  The  observation  of  SH  signal  may 
open  up  novel  possibilities  in  a  number  of  applications,  such  as  optical  interconnect,  auto-correlator  and 
spectrum  analyzer.  Finally,  we  have  fabricated  gain-guided  laser  structures  on  GaAs/Si  films  and 
observed  lasing  action  at  77K.  The  polarization  is  definitely  TM.  The  emission  spectrum  will  be 
analyzed,  for  the  first  time,  to  yield  information  on  the  internal  stress  in  GaAs/Si  films. 

In  summary,  we  have  demonstrated  the  feasibility  of  electro-optic  probing  for  three-dimensional 
field  mapping  and  for  SWR  measurements.  Earlier,  we  have  used  electro-optic  probing  to  ascertain  the 
cause  for  the  back-gating  effect  (see  Y.  H.  Lo,  Z.  H.  Zhu,  C.  L.  Pan,  S.  Y.  Wang,  and  S.  Wang,  “New 
Technique  to  Detect  the  GaAs  Semi-Insulating  Surface  Property  —  CW  Electro-Optic  Probing,”  Appl. 
Phys.  Lett.,  vol.  50,  p.  1125,  April  27,  1987).  One  fundamental  limitation  of  probing  an  electronic  cir¬ 
cuit  by  an  optical  beam  is  the  spatial  resolution  due  to  the  finite  spot  size.  The  spatial  resolution  is  not 
a  problem  for  measuring  SWR  and  for  studying  back-gating  effect,  but  is  a  problem  for  probing 
integrated  circuits.  Ways  to  overcome  this  problem  is  under  study.  In  the  mean  time,  exploratory  studies 
will  be  initiated  in  areas  related  to  GaAs/Si  hetero-epitaxy  to  enhance  the  chip  functionality  including 
opto-electronic  integration. 
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THEME  II  -  New  Architectures  for  Parallel  Computation 

MPC.l.  CAD  and  Architectural  Problems  in  Realizing  Massively  Concurrent  Computation: 

Acceleration  Techniques  for  Analog  Circuit  Simulation 
G.  Casino vi  and  A.  Sangiovanni-Vincentelli 

We  have  been  investigating  techniques  to  accelerate  circuit  simulation  of  large  scale  integrated 
circuits.  Our  interest  has  focused  on  two  ares:  aggregation  methods,  which  can  be  used  to  speed  up  the 
solution  of  large  systems  of  algebraic  equations,  and  macromodeling,  which  reduces  the  computational 
complexity  of  the  circuit  to  be  simulated. 

Aggregation  methods  are  used  to  improve  the  efficiency  of  relaxation-based  circuit  simulation, 
which  has  proved  very  effective  for  accurately  simulating  large-scale  digital  circuits.  Programs  such  as 
RELAX  and  SPLICE  are  based  on  relaxation  algorithms  that  have  guaranteed  convergence  properties. 
These  programs  can  analyze  large  circuits  with  the  same  accuracy  and  using  the  same  models  as 
SPICE2,  however  they  are  over  one  order  of  magnitude  faster  in  terms  of  CPU  time  used.  For  circuits 
with  tight  feedback  loops,  the  convergence  of  relaxation  algorithms  can  be  slow.  Aggregation  methods 
have  been  proposed  to  speed  up  relaxation  algorithms  used  in  the  solution  of  large  systems  of  algebraic 
equations. 

A  survey  of  aggregation  methods  was  compiled  by  F.  Chatelin  and  W.  L.  Miranker  [1],  who  also 
showed  that  all  the  existing  aggregation  methods  are  nothing  else  but  different  implementations  of  the 
same  mathematical  idea:  more  precisely,  aggregating  a  system  corresponds  to  performing  a  projection 
on  a  suitable  linear  subspace:  different  aggregation  methods  correspond  to  different  choices  of  sub¬ 
spaces. 

In  that  report,  the  authors  show  that  aggregation  methods  can  be  very  effective  in  improving  the 
performance  of  relaxation  processes.  At  the  same  time,  they  prove  that  the  choice  of  the  subspace 
plays  a  major  role  in  the  gain  that  can  be  obtained  from  an  aggregation  step;  it  follows  that  one  should 
choose  carefully  the  subspace  which  to  perform  the  projection  onto. 

We  developed  a  new  algorithm  that  uses  information  obtained  from  the  relaxation  process  to 
build  the  aggregated  system;  theoretical  analysis  and  numerical  experiments  performed  on  several 


matrices  show  that  this  algorithm  is  more  effective  than  the  other  aggregation  methods  proposed  so  far. 
The  results  of  this  research  were  presented  at  the  1985  International  Symposium  on  Circuits  and  Sys¬ 
tems  [2].  We  also  showed  that  our  method  can  be  applied  successfully  to  relaxation  based  circuit  simu¬ 
lation.  Our  conclusions  are  summarized  in  a  paper  that  will  appear  to  the  IEEE  Transactions  on 
Computer-Aided  Design  of  Integrated  Circuits  and  Systems  [3]. 

The  CPU  time  needed  for  cicruit  simulation  can  also  be  decreased  by  using  models  for  whole 
functional  blocks  rather  than  for  individual  devices.  Most  circuits  can  be  logically  divided  into  a  certain 
number  of  interconnected  functional  blocks,  and  in  many  instances  the  designer  is  not  interested  in  the 
detailed  waveforms  of  the  internal  nodes.  For  instance,  modularly  designed  systems  contain  several 
building  blocks  which  are  used  repeatedly.  As  a  consequence,  a  certain  amount  of  CPU  time  can  be 
saved  by  representing  each  block  with  a  macromodel  that  retains  approximately  the  same  input-output 
behavior  of  the  original  block,  while  discarding  all  redundant  information. 

A  natural  way  of  obtaining  a  macromodel  is  to  reduce  the  number  of  differential  equations 
describing  the  circuit,  trying  at  the  same  time  to  maintain  the  dynamics  of  the  original  system.  This 
approach  requires  us  to  be  able  to  quantify  the  intuitive  notion  that  "two  dynamical  systems  have 
approximately  the  same  input-output  behavior".  In  addition  to  that,  we  also  need  an  algorithm  that  can 
tell  us  how  to  modify  a  given  macromodel  in  order  to  make  its  dynamical  behavior  as  close  as  possible 
to  the  one  of  the  original  system. 

The  first  problem  can  be  solved  by  taking  the  difference  of  the  outputs  of  the  original  system  and 
the  macromodel  for  the  same  input,  and  computing  its  L2  norm  over  a  fixed  time  interval.  The 
supremum  of  that  norm,  as  the  input  varies  over  a  given  class  of  "admissible  inputs",  can  be  taken  as  a 
numerical  indication  of  how  much  the  two  systems  differ  in  their  dynamical  behavior.  From  a 
mathematical  point  of  view,  this  is  a  classical  optimal  control  problem,  and  in  the  literature  there  are 
numerous  algorithms  that  can  be  used  for  its  numerical  solution. 

In  sharp  contrast  to  this,  only  few  and  limited  results  have  been  published  answering  the  question 
of  how  the  dynamics  of  a  system  are  affected  by  changes  in  the  set  of  differential  equations  that 
describe  it.  In  fact,  we  can  show  that  this  problem  can  be  dealt  with  in  a  very  general  setting  by  a  suit- 


able  generalization  of  the  classical  optimal  control  problem.  Moreover,  we  already  have  a  mathematical 
tool  available  that  can  be  used  to  solve  it,  both  from  a  theoretical  and  computational  point  of  view.  This 
tool  is  the  theory  of  hamiltonian  systems,  which  has  been  studied  extensively  because  of  its  applications 
to  many  problems  of  classical  and  reladvistic  physics. 

Applying  those  theoretical  results  to  our  case,  we  have  developed  an  algorithm  that  can  be  used 
to  modify  a  macromodel  and  make  its  dynamics  fit  those  of  the  original  system  as  closely  as  possible. 
We  are  now  implementing  the  algorithm  in  a  form  suitable  for  modeling  actual  circuit  The  program 
MODEL  will  use  NUTMEG  as  input  and  output  interface,  giving  it  full  compatibility  with  SPICE  3  and 
other  simulators  using  the  same  language. 


References 

[1]  F.  ChateUn  and  Wl.  Miranken  Acceleration  by  aggregation  of  successive  approximation 
methods,  IBM  Thomas  J.  Watson  Research  Center  Report  No.  RC  8508,  1980. 

[2]  G.  Casinovi  and  A.  Sangiovanni-Vincentelli:  Aggregation  Methods  for  the  Solution  of  Large 
Scale  Systems  of  Linear  Algebraic  Equations.  Proceedings  of  the  1985  International  Symposium 
on  Circuits  and  Systems.  Kyoto,  Japan,  June  1985. 

[3]  G.  Casinovi  and  A.  Sangiovanni-Vincentelli:  A  New  Aggregation  Technique  for  the  Solution  of 
Large  Systems  of  Algebraic  Equations.  To  appear  on  IEEE  Transactions  on  Computer-Aided 
Design  of  Integrated  Circuits  and  Systems. 


Our  research  has  focussed  on  three  projects:  Massive  Parallel  Computation  with  particular  emphasis  on 
neural  networks,  and  macro-modeling  for  large-scale  analog  circuit  simulation  and  placement  of 
modules  for  large-scale  chips. 

Scalability  in  Feedforward  Neural  Networks 
A.  Kramer  and  A.  Sangiovanni-Vincentelli 

Neural  networks  have  received  a  great  deal  of  attention  recently  as  a  promising  new  model  of 
parallel  computation.  While  there  is  still  a  great  deal  of  variation  in  definitions  of  what  comprises  a 
neural  network,  in  general  these  all  involve  highly  interconnected  networks  of  simple  analog  computing 
nodes.  The  output  of  each  of  these  nodes  is  a  non-linear  function  of  the  weighted  sum  of  its  inputs. 
Typically  these  nonlinear  functions  are  either  step  functions  or  sigmoids.  The  computational  function 
of  a  neural  network  is  fully  described  by  the  nonlinear  function  used  by  its  nodes,  the  topology  of  the 
interconnections  between  nodes,  and  the  "weights'*  on  these  interconnections. 

While,  in  general,  neural  networks  may  have  feedback  paths  (a  directed  cycle  on  the  graph 
representing  its  topology),  this  results  in  a  dynamical  system  which  is  very  complex  to  analyze.  Feed¬ 
forward  networks  are  a  restricted  class  of  network  architectures  in  which  nodes  are  organized  into 
sequential  "layers"  and  connections  are  made  only  from  a  node  to  nodes  on  higher  layers.  The  first 
layer  is  the  "input"  layer,  the  last  is  the  "output"  layer,  and  any  layer  in  between  are  "hidden"  layers 
(typically  only  1).  This  class  of  networks  is  interesting  because  the  dynamics  of  the  network  do  not 
affect  the  function  it  is  computing  (and  hence  dynamics  may  be  ignored),  and  the  class  of  functions 
computable  by  feedforward  networks  appears  to  be  quite  rich. 

The  usefulness  of  FFNN’s  really  involves  three  separate,  though  related  questions: 

1)  What  functions  are  representable  by  FFNN’s 

2)  How  easy  is  it  to  find  these  representations  (teaming) 

3)  How  efficient  are  these  representations 

While  there  are  many  groups  investigating  representation  and  learning  in  FFNN’s,  the  question  of 
efficiency  is  usually  ignored.  This  is  unfortunate,  because  it  is  unlikely  that  the  technology  of  neural 
networks  will  be  developed  until  the  computations  being  performed  in  them  are  known  to  be  more 
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efficient  than  solutions  that  use  existing  computational  technology  (von  Neuman,  systolic  array,  etc...). 

There  are  many  ways  to  pose  the  question  of  efficiency  in  FFNN’s,  including  power  consumption, 
physical  size,  accuracy,  robustness,  etc...  More  important  than  any  of  these,  however,  is  the  question  of 
scalability:  how  does  the  network  representation  of  a  function  "scale"  with  function  "size"  (dimensional¬ 
ity).  This  metric  is  known  for  conventional  models  of  computation,  and  it  is  a  relevant  and  well-defined 
measure  by  which  to  compare  the  efficiency  of  computation  in  FFNN’s. 

For  most  interesting  classes  of  functions,  the  question  of  scalability  cannot  be  answered  analyti¬ 
cally,  and  for  this  reason,  any  investigation  of  scalability  must  address  issue  of  leamability.  For  these 
functions,  the  only  way  to  determine  whether  a  given  network  is  capable  of  representing  a  particular 
function  is  to  try  to  find  such  a  representation  via  a  "learning  algorithm,"  typically  gradient  descent  or 
backpropogation.  These  algorithms  search  the  "weightspace"  of  a  network  for  a  point  (a  set  of  weight 
values)  at  which  the  "error"  between  the  target  function  and  the  network  function  is  equal  to  zero. 

A  question  related  to  that  of  scalability,  which  is  of  interest  for  the  reasons  just  given,  is  that  of 
the  relationship  between  scalability  and  leamability:  for  a  function  of  a  given  size,  does  the  network 
representation  of  the  function  become  easier  to  find  as  the  network  size  is  increased?  (Does  learning 
become  easier  as  the  size  of  the  hidden  layer  is  increased?).  Clearly  any  measure  of  "leamability" 
depends  on  the  learning  algorithm  being  used.  Despite  this,  because  most  iterative  improvement  algo¬ 
rithms  have  problems  searching  spaces  where  "bad"  (unacceptable)  local  minima  abound,  we  have 
chosen  "abundance/badness  of  local  minima"  as  a  useful  and  semi-universal  measure  of  the  "leamabil¬ 
ity."  "Abundance"  we  measure  probabilistically,  while  "badness"  is  a  description  of  a  local  minima 
which  includes  depth,  basin  of  attraction,  etc... 

It  should  be  clear  that  our  investigation  involves  a  great  deal  of  simulation  of  FFNN’s.  In  particu¬ 
lar,  while  most  existing  research  has  concentrated  on  simulation  of  small  "toy"  networks,  we  are 
interested  in  scalability  and  thus  need  to  be  able  to  simulate  large  networks  in  reasonable  times.  Prohi¬ 
bitively  long  simulation  times  are  the  major  reason  past  work  has  not  addressed  scalability.  On  a  con¬ 
ventional  sequential  computer,  FFNN  simulation  is  0( Anodes  squared)  and  this  is  unacceptable  for  large 
simulations.  For  this  reason,  we  have  developed  a  massively  parallel  FFNN  simulator  which  runs  on  the 
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Connection  Machine.  This  SIMD  computer  has  64,000  1-bit  processors  connected  in  a  hypercube.  Our 
simulator  on  the  CM  runs  in  0(log(#nodes))  and  can  simulate  networks  of  up  to  250  nodes  per  layer. 
Currently,  the  simulator  is  in  the  final  stages  of  testing. 


We  also  need  function  classes  which  are  interesting  and  scalable.  Currently,  we  have  several 
classes  in  mind: 

1)  random  boolean  functions 

2)  structured  boolean  functions  (pin,  etc...) 

3)  real-valued  functions 

4)  random  network  functions  (a  random  set  of  weight  values) 

5)  decoding  of  error-correcting  codes 

6)  real-world  problems  (ala  expert  systems) 

7)  recognition  of  geometric  patterns 

For  the  moment,  we  are  concentrating  on  boolean  functions. 


State  -  Analysis  of  Simple  Neural  Nets 
CJL  Sdquin  and  C.  Ch inrun grueng 

1.  Motivation 


Neural  networks  —  here  understood  to  mean  highly  interconnected  networks  of  some  kind  of  sum¬ 
ming,  thresholding,  or  switching  elements  —  provide  an  intriguing  alternative  to  other  computing  circuitry 
such  as  the  classical  von  Neumann  computer  or  systolic  signal  processing  pipelines.  The  true  potential  of 
neural  networks  and  their  most  fruitful  areas  of  application  are  still  a  topic  of  hot  debate.  This  debate  will 
continue  undl  one  has  developed  a  quantitative  understanding  of  their  properties  such  as  knowledge  storage 
capacity,  information  retrieval  speed,  and  learning  performance. 

One  goal  of  our  research  is  to  fully  map  out  the  error  surfaces  of  a  few  very  simple  networks  and  to 
make  them  somehow  visible  to  the  reader.  These  error  surfaces  are  given  as  some  measure  of  the  overall 
deviation  of  the  actual  outputs  of  the  network  from  the  desired  outputs  (summed  over  all  input  patterns)  as 
a  function  of  all  possible  states  of  the  network  as  given  by  the  values  of  all  its  (synaptic)  weights.  It  is  our 
belief  that  from  these  simple  examples  we  can  gain  some  insight  into  the  structure  of  the  error  surfaces  of 
more  complicated  neural  networks.  A  parallel  goal  is  to  develop  new  visualization  techniques  for  these 
high-dimensional  error  surfaces. 

2.  The  Networks 

The  problem  we  have  concentrated  on  so  far  is  the  Boolean  exclusive-or  (XOR)  function  of  two 
inputs.  This  is  one  of  the  simplest  problems  that  cannot  be  solved  with  a  simple  Perceptron  with  a  single 
threshold  element  (Fig.  la).  If  we  disallow  feedback  within  the  same  layer,  then  we  need  a  feed-forward 
network  with  at  least  two  layers  of  switching  elements  to  construct  a  network  that  can  perform  the  XOR 
function. 

The  particular  networks  considered  in  this  study  are  shown  in  Fig.  1.  The  elements  in  the  input  layer 
simply  provide  an  impedance  match  to  the  input  signals,  limit  the  signals  to  the  amplitude  range  considered 
at  all  the  other  switching  elements,  and  provide  unlimited  fanout  The  switching  elements,  or  ‘neurons’  in 
the  middle  layer,  also  called  the  ‘hidden’  elements,  are  connected  to  some  or  all  of  the  input  elements 
through  adjustable  weights.  Their  outputs  are  not  directly  observable,  but  provide  —  again  through 


adjustable  weights  —  inputs  for  the  neurons  in  the  output  layer  (just  one  in  our  case).  To  feed  an  input 
directly  to  an  output  neuron,  one  of  the  hidden  neurons  can  be  replaced  with  a  dummy  element  with  a  sin¬ 
gle  input  with  a  fixed  weight  of  one  (see  h3  in  Fig.  1c). 


Figure  1.  Three  networks  used  in  our  analysis  and  their  labeling  of  the  weights: 

(a)  Perceptron,  (b)  regular  3 -layer  net,  (c)  net  with  missing  hidden  units. 

The  neurons  themselves  can  take  several  different  implementations.  Common  to  all  of  them  is  the 
fact  that  they  form  a  weighted  sum  of  all  their  inputs  multiplied  with  the  corresponding  weights  on  those 
connections.  In  addition,  a  bias  term  that  does  not  depend  on  any  input  is  added  to  the  sum;  in  effect  this 
permits  one  to  adjust  the  switching  threshold.  This  sum  is  subjected  to  a  ‘switching’  function  that  typically 

takes  the  form  of  a  sharp  thresholding  function  or,  in  our  case,  a  smooth  monotonic  sigmoid  function.  We 

have  assumed  these  functions  to  be  symmetrical  with  respect  to  the  origin  in  order  not  to  camouflage  any 
inherent  symmetries  that  arise  from  the  symmetry  of  the  network  or  of  the  specified  problem.  We  also  nor¬ 
malize  the  amplitudes,  so  our  network  inputs  and  outputs  from  any  neuron  will  be  in  the  range  [-1  . .  +1]. 

3.  Error  Surfaces 

Let  us  consider  the  network  shown  in  Fig.  lb.  The  state  of  this  network  is  fully  defined  by  the  nine 
adjustable  weights  wl  through  w9;  thus  the  state  space  of  this  network  is  9-dimensional.  If  these  weights 
are  known,  we  can  compute  the  network  output  for  each  pattern  of  inputs  assuming  a  smooth  sigmoid 
switching  function  swf  inside  the  neurons: 

h  l(i  1  j'2)  -  swf  (  w4  +  w5  x  il  +  w6  x  i2 ) 
h 2 (i  1 4 2)  =  swf  (  w  7  +  w  8  x  i  1  +  w  9  x  i  2 ) 

0  l(i  l,i2)  =  swf  (wl  +  w2xhl  +  w3  x  h2 ). 

If  we  further  assume  that  for  each  pair  of  inputs  the  goal  function  goal(il  j2)  is  known,  then  we  can  also 
compute  an  error  for  each  input  pattern: 
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err  (i  lj2)  =  goal  (»'  1 , 1 2)  -  o  1  (i  1  ,i  2). 

Finally  we  can  combine  all  these  error  terms  into  a  single  term  that  reflects  how  well  the  network  overall 
performs  its  designated  function.  We  are  really  interested  in  the  maximum  error  made  for  any  of  the 
specified  input  patterns,  so  a  function  like  MaxErr  might  be  most  appropriate: 

MaxErr  =  max  (  err  (i  1,  i  2)  )  for  all  legal  input  patterns  (i  1,  i  2). 

However,  computationally  it  is  often  preferable  to  have  a  smooth  function  to  work  with  during  the  analysis. 

Thus  often  a  total  square  error,  TsqErr ,  is  computed  by  summing  the  squares  of  the  individual  errors: 

TsqErr  =  T  {err  (i  1 J.  2f. 
iU2 

For  a  given  network  and  a  given  problem  statement,  such  an  error  surface  Err(w  1  •  •  •  w9),  defined  over 
the  domain  of  the  state  space,  fully  describes  the  possible  states  of  the  net  and  how  well  it  performs  for  any 
combination  of  values  of  its  weights.  If  we  could  globally  examine  this  surface,  we  could  readily  answer 
the  questions  how  well  a  network  could  perform  a  given  task.  Unfortunately  these  surfaces  can  be  of  very 
high  dimension  —  for  this  simple  example  it  is  already  a  9-dimensional  hyper-surface  embedded  in  10- 
dimensional  space. 

For  the  simplest  cases  we  can  actually  model  the  error  surface.  This  is  illustrated  with  the  trivial  net¬ 
work  given  in  Fig.  la.  While  the  system  basically  has  three  weights  and  thus  three  degrees  of  freedom,  the 
symmetry  of  the  problem,  i.e.  the  equivalence  of  the  two  inputs  allows  us  to  assume  that  w2  and  w3  are 
identical  and  then  restrict  us  to  exploring  the  2-dimensional  manifold  produced  as  a  function  of  wl  and  of 
w2=w3. 


Figure  2.  Error  surfaces  for  XOR  function  on  net  la. 

Figure  2  shows  two  surfaces  for  this  network  if  the  desired  function  is  the  Boolean  XOR  operation 
between  the  two  inputs.  Fig.  2a  shows  the  normalized  total  error  E  (x  ,y  )='UJsqErr  /4)  summed  over  all 
four  different  input  patterns.  The  range  of  the  two  weight  variables  is  from  -3.5  to  +3.5  each.  In  Fig.  2b  we 
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have  plotted  instead  the  maximum  error  for  any  of  the  four  patterns.  Notice  that  no  part  of  this  surface 
comes  down  near  the  ground  plane  as  would  be  indicative  of  a  possible  solution.  This  is  just  another  way  of 
realizing  that  a  single  neuron  (a  Perceptron)  cannot  perform  the  XOR  function. 


'(a)  TsqErr 


(b)  MaxErr 


Figure  3.  Error  surfaces  for  XOR  function  on  net  lc. 

If  we  add  a  second  neuron  in  the  middle  (hidden)  level  (Fig.  lc)  then  we  can  find  a  set  of  weights 
that  can  produce  a  solution  to  the  XOR  problem.  In  Figure  3  the  error  function  is  displayed  as  a  function  of 
the  threshold  wl  and  the  weight  w2.  The  diagonal  valley  goes  indeed  deep  enough  to  constitute  a  solution 
if  the  weight  range  is  increased  to  ±  6.0.  The  other  weights  are  set  to:  w4=3.0  and  w5=w6=3.2. 

The  network  in  Fig.  lb  is  already  complicated  enough  that  a  single  surface  in  3-dimensional  space  is 
inadequate  to  represent  its  solution  space.  In  this  case  we  need  to  rely  on  more  mathematical  techniques  to 
obtain  information  about  its  error  surface  and  possible  solution  states. 


4.  Symmetries 

Often  the  state  space  of  a  given  network  for  a  given  problem  is  rather  symmetric.  This  can  reduce  the 
scope  of  the  visualization  task.  The  net  in  Figure  lb  has  two  equivalent  hidden  neurons  that  are  connected 
in  the  same  manner  to  the  rest  of  the  network.  Thus  the  role  of  these  two  neurons  could  be  interchanged 
without  affecting  the  overall  i/o  behavior  of  the  net  as  a  whole.  Specifically  that  means  the  function  of  the 
network  will  be  unchanged  if  we  make  the  following  four  substitutions: 

opA:  w2  < — >  w3,  w4  *  >  w7,  w5  < — >  w8,  and  w6  < — >  w9. 

Similarly,  we  have  identified  five  more  symmetry  operations: 
opB:  w5  4 — >  w6  and  w8  < — *  w9. 
opC:  w2  — » -w2,  w4  — >  -w4,  w5  — >  -w5,  and  w6  — >  -w6. 
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opD:  w3  — »  -w3,  w7  — >  -w7,  w8  -»  -w8,  and  w9  ->  -w9. 

opE:  w5  -»  -w5,  w8  -*  -w8;  and  wl  ->  -wl,  w2  ->  -w2,  w3  — >  -w3. 

opF:  w6  -»  -w6,  w9  -» -w9;  and  wl  -»  -wl,  w2  ->  -w2,  w3  -»  -w3. 

These  six  operations  result  in  a  64-fold  symmetry  for  the  state  space  of  net  lb  computing  the  XOR  func¬ 
tion. 

Such  symmetries  are  very  important  to  discover  because  they  reduce  the  fraction  of  the  state  space 
that  needs  to  be  explored  in  order  to  gain  a  full  understanding  of  the  behavior  of  the  neural  net.  We  thus 
only  need  to  sample  a  suitably  defined  region  in  the  9-dimensional  space  that  amounts  to  l/64th  of  the  total 
domain. 


5.  Sampling  Results 

Within  one  representative  domain,  we  now  analyze  the  function  by  sampling  it  at  random  with 
several  hundred  arbitrary  starting  points,  from  where  we  start  a  steepest  descent  algorithm  to  find  the 
corresponding  local  minimum.  It  turns  out  that,  for  a  representative  domain  of  l/64th  of  that  total  state 
space,  we  find  only  three  separate  local  minima.  The  error  TsqErr  and  the  corresponding  values  of  the 
synaptic  weights  are  listed  in  Table  1. 


Table  1.  Minima  in  XOR  Function  on  Net  3b. 

Mult. 

Error 

wl 

w2 

w3 

w4 

w5 

w6 

w7 

w8 

w9 

4 

0.080 

2.86 

-3.00 

-3.00 

3.00 

-3.00 

3.00 

3.00 

3.00 

-3.00 

2 

2.097 

0.00 

-3.00 

3.00 

3.00 

-1.06 

3.00 

3.00 

1.06 

3.00 

2 

2.884 

0.82 

-0.88 

-0.74 

3.00 

-3.00 

3.00 

3.00 

-3.00 

3.00 

There  is  only  one  'good*  minimum  that  has  an  error  value  that  is  low  enough  so  that  the  state  of  this 
minimum  can  be  considered  a  solution  of  the  network  for  the  XOR  problem. 

Because  of  the  64-fold  symmetry  of  the  st?  space  for  this  problem,  each  such  minimum  will  be 
replicated  in  64  locations.  However  some  of  these  replicas  may  fall  on  top  of  one  another.  The  first  column 
in  the  above  table  indicates  the  multiplicity  of  the  observed  minima.  The  good  minimum  (line  1)  has  a  mul- 


WW'.VTN 


wffn. 


tiplicity  of  4,  and  thus  there  will  be  a  total  of  16  such  locations  in  the  whole  state  space. 

6.  Ongoing  Work 

We  are  now  in  the  process  of  determining  the  shape  of  the  various  local  minima.  In  particular  we  are 
trying  to  determine: 

•  The  maximum  extension  of  the  collection  zone,  i.e.,  what  is  the  farthest  known  source  point  from  the 
final  local  minima  point  where  it  comes  to  rest  ? 

•  The  maximum  radius  of  complete  collection,  i.e„  what  is  the  largest  9-dimensional  sphere  around  a 
local  minima,  inside  which  all  starting  points  will  converge  to  this  minimum  ? 

•  The  minimum  barrier  height  for  escape  from  any  minimum,  Le.,  what  are  the  lowest  barrier  heights 
between  any  pair  of  adjacent  minima  ? 

Once  this  static  analysis  of  the  error  surface  is  completed,  we  will  then  study  the  dynamic  aspects  of 
the  learning  behavior  of  these  networks.  This  will  involve  an  analysis  of  the  geodesic  trajectories  that  the 
system  takes  when  it  converges  to  one  of  the  local  minima  under  the  gradient  descent  or  the  back- 
propagation  algorithm. 


PROUD:  An  Efficient  Sea-of-Gates  Placement  Algorithm  for  Dense  Chip 
Ren-Song  Tsay  and  E.  S.  Kuh 

As  Application  Specific  Integrated  Circuits  (ASIC)  become  more  popular  and  the  gate  count  on  a 
chip  gets  increasingly  larger,  the  quality  and  speed  of  automatic  layout  algorithms  need  to  be  read¬ 
dressed.  It  is  well-known  that  the  problem  of  module  placement  is  especially  crucial  to  the  final  out¬ 
come  of  the  design.  For  a  sea-of-gates  chip  with  100K  or  more  gates,  conventional  deterministic 
methods  are  no  longer  suitable  because  sparsity  of  the  connectivity  specification  ought  to  be  considered. 
Random  algorithms  such  as  simulated  annealing  could  take  a  forbidden  amount  of  time  in  reaching  an 
acceptable  solution. 

We  introduce  a  hierarchical  method  of  placement  which  takes  full  advantage  of  the  sparsity 
inherent  in  the  placement  specification.  The  method  starts  from  a  quadratic  placement  formulation  and 
takes  the  I/O  pad  specification  as  input  and  solves  successively  linear  sparse  equations.  The  slot  con¬ 
straints  are  considered  gradually  by  a  simple  partitioning  on  the  global  placement  results  solved  from 
linear  equations.  After  partitioning,  we  use  an  interesting  global  optimization  technique  to  handle  the 
inter-influence  among  separated  regions.  From  the  result  of  several  real  chips,  the  quality  of  the  place¬ 
ment  is  excellent  It  is  superior  to  that  obtained  by  TimberWolf3.2  and  is  comparable  with  that  of  Tim- 
berWolf4.2.  The  run  time  complexity  of  our  method  is  0  (n  log2*)  and  the  memory  space  complexity 
is  linear.  For  a  100K  sea-of-gates  chip  with  26,000  modules,  the  run  time  is  about  50  minutes  on  a 
VAX  8650  (6  MIPs  machine)  and  the  memory  requirement  is  about  11  megabytes. 

The  algorithm  has  been  extended  to  do  placement  on  designs  with  mixed  big  macros  and  standard 
cells.  We  also  find  a  solution  for  I/O  pad  assignment  from  the  general  quadratic  formulation  but  using 
eigenvector  approach  instead  of  solving  linear  equations. 

Parallel  implementation  of  the  PROUD  placement  algorithm  is  now  under  study. 
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MPC.2.  Systems  Problems  in  Parallel  Distributed  Computation 
Shobana  Venkatraman 
(Professor  Shankar  Sastry) 

The  convergence  and  adaptation  properties  of  neural  networks,  their  memory  capacities  and  pro¬ 
grammability  have  not  been  studied  completely  thus  far.  We  are  investigating  these  properties  for  a 
general  class  of  neural  networks  of  which  the  Hopfield  and  Grossberg  networks  are  special  cases.  We 
have  found  that  in  such  networks  it  may  be  possible  to  establish  very  clearly  defined  basins  of  attraction 
around  each  of  the  equilibria  (the  nodes  we  wish  to  “remember”).  We  have  also  been  studying  the 
effect  of  the  hamming  distance  between  equilibria  on  the  location  of  the  spurious  equilibria,  if  any,  that 
are  generated.  We  are  working  on  making  the  algorithm  for  the  synthesis  of  this  class  of  neural  net¬ 
works  more  programmable,  in  the  sense  of  being  able  to  update  a  network  built  to  remember  *m’  nodes 
to  store  ‘m+r  equilibria  without  having  to  restructure  the  entire  network.  We  are  studying  the  charac¬ 
teristics  of  these  artificial  neural  networks  by  simulating  the  nonlinear  differential  equations  that 
describe  them. 
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MPC  J.  Interconnection  of  Dense  Computational  and  Memory  Elements: 


Implementation  of  Learning-Capable  Neural  Nets  in  an  Analog  MOS  VLSI  Environment 

Victor  Hu  and  Alan  Kramer 

(Professors  Ping  K.  Ko  and  Chenming  Hu) 

The  goal  of  this  project  is  to  investigate  the  feasibility  of  implementing  leam-capable  neural  nets 
using  analog  MOS  VLSI  technology.  While  there  is  a  great  deal  of  variety  in  approaches  to  neural-net 
computation,  most  proposed  architectures  are  characterized  by  the  weighted  interconnection  of  simple 
nonlinear  threshold  computing  elements.  These  connection  weights,  which  define  the  computational  con¬ 
tents  and  functions  of  the  network,  are  considered  one  of  the  most  critical  elements  in  its  practical 
implementation. 

The  computational  capabilities  of  these  networks  have  been  demonstrated  using  pre-programmed 
resistors  or  digitally  weighted  resistors  in  the  interconnections.  However,  the  real  potential  of  neural 
computation  lies  with  a  fully  analog  and  leam-capable  network.  We  are  exploring  the  possibility  of 
using  field  alterable  resistive  elements  in  this  network. 

In  the  past  few  months,  we  have  focused  on  the  use  of  the  floating-gate  EEPROM  devices.  We 
have  evaluated  the  analog  programming  characteristics  of  commercial  EEPROM  devices  as  well  as  test 
devices  we  obtained  through  our  industrial  contacts.  The  results  so  far  have  been  encouraging  and  sug¬ 
gest  that  analog  programming  of  a  EEPROM  device  is  indeed  feasible  and  controllable. 

In  order  to  build  a  MOS  neural  net  chip,  however,  a  integrated  MOS/EEPROM  technology  is 
necessary.  We  are  currently  designing  an  inhouse  technology  with  that  capability.  A  NMOS/EEPROM 
test  run  will  be  started  once  the  test  mask  design  is  completed.  Several  designs  of  a  simple  demonstra¬ 
tion  neural-net  chip  using  the  back-propagation  algorithm  are  being  investigated. 


High-Speed  MOSFET/BJT  Technology 
James  Moon 

(Professors  Ping  K.  Ko  and  Chenming  Hu) 

We  are  developing  a  high-speed  BJT/CMOS  technology  for  the  implementation  of  analog  and 
digital  circuits. 

For  the  MOS  transistors,  this  technology  features  Elevated  Drain  and  Source  (EDS)  regions  that 
are  grown  epitaxially  from  the  substrate,  and  the  Contact-Over-Oxide  (COO)  structure  developed  earlier 
in  our  laboratory.  The  COO  contact  scheme  occupies  a  much  smaller  area  than  conventional  junction 
contacts,  has  significantly  lower  junction-to-substrate  capacitances,  and  does  not  require  a  diffusion  bar¬ 
rier  layer  for  shallow  junctions.  In  addition,  this  structure  also  promise  greater  radiation  hardness  to 
SEU  (Single  Event  Upset)  than  conventional  techniques.  With  the  EDS  scheme,  an  LDD  (Lightly 
Doped  Drain)  structure  for  hot-electron  protection  can  be  realized  simply  by  controlling  the  implant 
energies  of  the  n-  and  n+  regions,  instead  of  using  the  more  difficult  oxide-spacer  technology.  One  of 
the  goals  in  our  device  design  is  to  make  the  devices  operational  even  at  channel  lengths  down  to  0.1  to 
02  pm.  For  the  bipolar  transistor,  we  are  investigating  the  designs  of  both  lateral  and  vertical  injection 
structures.  Simplifying  the  device  structure  and  making  it  MOS  technology  compatible  are  our  main 
goals. 

Since  the  project  started  six  months  ago,  we’ve  developed  a  deep-submicron  lithographic  process 
which  is  based  on  conventional  optical  lithographic  in  conjunction  with  a  novel  photoresist  ashing  tech¬ 
nique.  Using  this  process,  we  are  able  to  fabricate  FETs  with  controlled  channel  lengths  in  the  1/4  pm 
regime  routinely.  Some  of  the  important  achievements  so  far  are:  (1)  An  n-channel  NMOS  transistor 
with  channel  length  of  0.15  pm  and  transconductance  of  650  mS/mm  at  room  temperature  -  the  highest 
reported  for  the  MOS  transistor.  (2)  A  propagation  delay  of  21  ps  per  stage  has  been  achieved  for  a 
NMOS  depletion-load  logic  gate  chain  at  room  temperature  -  the  fastest  reported  for  a  MOS  logic  cir¬ 
cuit 

We  have  fabricated  n-channel  MOSFETs  with  large  variations  in  device  dimensions  and  process 
parameters.  A  comprehensive  study  was  performed  on  the  characteristics  of  these  devices.  Important 
scaling  guidelines  for  device  performance  and  reliability  have  been  established. 


-43- 


At  present,  we  are  focusing  on  various  technology  integration  issues.  We  have  also  begun 
designing  test  devices  and  circuits  for  evaluating  the  capabilities  of  the  technology. 


Process  Simulation  for  Extraction  of  Topography  Dependent  Electrical  Characteristics: 

Application  to  VLSI  Neural  Networks 
Edward  Scheduler 
(Professor  A.  R.  Neureuther) 

The  implementation  of  highly  interumnected  computational  structures,  such  as  neural  networks,  is 
critically  dependent  on  minimizing  performance  degradation  from  parasitic  loading  effects.  The  layout 
and  topographical  features  in  processing  determine  the  magnitude  of  parasitic  loading.  Efforts  to  link 
profiles  generated  by  the  SIMPL  and  SAMPLE  process  simulators  to  programs  for  device  and  RC- 
parasitic  modeling  are  being  pursued  in  order  to  investigate  questions  of  VLSI  neural  network  perfor¬ 
mance. 

Simulations  of  prototype  networks  have  yielded  basic  information  on  parasitic  effects.  Some  gen¬ 
eral  trends  regarding  speed  limitations  and  scaling  effects  with  larger,  higher  density  networks  have 
been  identified  for  Hopfield  networks  based  on  the  implementation  developed  at  AT&T  Bell  Labs. 
Investigation  of  other  architectures,  including  prototypes  being  developed  here  by  Professor  P.  K.  Ko,  is 
being  pursued.  General  trends  in  neural  network  research  are  being  followed  to  identify  which  areas 
are  most  significantly  affected  by  VLSI  issues. 
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Abstract — It  Is  shown  that  the  substrate  current  characterization 
method  and  modeling  approach  used  for  n-MOSFET’s  is  also  applicable 
to  p-MOSFETs.  The  impact  ionization  rale  extracted  for  holes  is  found 
to  be  I  x  10*  exp  ( -  3.7  x  10 */£),  where  £  is  the  electric  field.  Based  on 
oar  measurement  and  modeling  result,  roughly  twice  the  channel  electric 
field  la  required  for  p-MOSFET's  to  generate  the  same  amount  of 
substrate  current  as  n-MOSFET's.  The  hol-canicr*iaduced  breakdown 
voltage  is  therefore  also  about  two  times  larger. 


I.  Introduction 

IT  HAS  LONG  been  recognized  that  hot-carrier  effects 
impose  a  limit  upon  future  device  scaling.  Hot-carrier 
effects  in  n-MOSFET's,  which  can  be  correlated  with  the 
substrate  current  [I],  are  more  serious  than  in  p-MOSFET's. 
Models  of  the  substrate  current  in  n-MOSFET’s  have  been 
well  developed  12],  [3].  While  there  have  been  studies  of  hot- 
carrier  effects  in  p-MOSFET’s  [4J-[6],  it  is  not  clear  that  the 
substrate  current  model  for  n-MOSFET's  is  also  applicable  to 
p-MOSFET's.  In  this  paper,  we  show  that  the  channel  electric 
fields  and  the  substrate  currents  in  n-MOSFET's  and  p- 
MOSFET's  can  be  modeled  by  the  same  physical  approach, 
although  the  values  of  the  semiempirical  parameters  are 
different.  The  test  devices  used  in  this  study  are  n*  polysili¬ 
con-gate  MOSFET's  without  threshold  adjust  implant.  Three 
kinds  of  gate  oxide  thicknesses  were  fabricated  for  both  n-  and 
p-MOSFET's:  432,  152,  and  90  A .  The  substrate  concentra¬ 
tion  is  2  X  10“  cm  "3  for  all  the  n-MOSFET's  and  4  x  10“ 
cm*3  for  all  the  p-MOSFET’s.  The  sourcc/drain  junctions  are 
0.3-pm  arsenic  junctions  in  the  n-MOSFET's  and  0.4 -pm 
boron  junctions  in  the  p-MOSFET’s. 

II.  Channel  Field  and  Impact  Ionization  Coefficients 

It  has  been  shown  that  for  conventional  n-MOSFET’s  the 
maximum  lateral  channel  Held  at  the  drain  end  Em  can  be 
expressed  in  terms  of  the  drain  voltage  VD  as  [7],  [8] 


*,  \(  Kd”  Kdsat  V  .  !-•>  ~  V° SAT 

En^\ - 7 - )  +E™ - 7 - 


/  =  0.22r°J3A’?5 


(1) 


(2) 


saturation  velocity.  The  corresponding  drain  voltage  is  VD 
T<u  is  the  gate  oxide  thickness  and  Xj  is  the  source/drain 
junction  depth.  £Sat  is  about  4  x  104  V/cm  for  n-MOSFET’s. 
Using  this  relationship,  the  substrate  current  /SUB  can  be 
expressed  as 


'SUB1 


IoAiEJ 

I  '■  ■  ■  — ■ 

B, 


exp  (-Bj/Em) 


(3) 


where  Esat  is  the  channel  field  at  Ahicn  the  earners  reach 
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where  A,  and  B,  are  the  constants  in  the  expression  for  the 
impact  ionkaiion  rate,  i.e.,  a  *  A(  exp  (-  Bi/E),  and  Iq  is 
the  drain  current.  Equations  (1M3)  were  derived  and  have 
been  verified  by  experiments  and  simulations  for  n-MOS¬ 
FET’s  [8],  However,  there  is  nothing  in  the  derivations  to 
suggest  that  they  are  not  equally  applicable  to  p-MOSFET’s. 

VD  jut  can  be  determined  experimentally  for  n-MOSFET’s 
[9],  The  same  method  is  applied  to  p-  and  n-MOSFET’s  as 
shown  in  Fig.  1.  For  p-MOSFET’s,  parallel  shift  of  the 
constant  hud  Id  locus  is  observed  and  Vo  sat  can  be  experi¬ 
mentally  determined  as  well.  From  (l)-(3),  a  plot  of  In  (/su*/ 
(ID( Vd  -  VD sat)))  versus  l/(VD  -  Vdsat)  will  yield  one 
single  straight  line  for  all  the  devices  and  bias  voltages.  The 
vertical  intercept  and  slope  of  this  line  give  At/Bi  and  Bf, 
respectively.  In  Fig.  2,  such  a  plot  is  shown  for  both  the  n- 
MOSFET’s  and  the  p-MOSFET’s.  Using  E^t  -  1.2  x  10j 
V/cm,  a  single  straight  line  also  is  obtained  for  different  sizes 
of  p-MOSFET's  at  different  biases.  This  proves  that  (1M3) 
are  applicable  to  p-MOSFET’s,  too.  Using  (2)  for  /,  A,  and  B, 
deduced  from  Fig.  2  are 

n-MOSFET:  fi,- 1.7x10*  V  •  cm*',  A,»2x  10‘ cm'1 

p-MOSFET:  3.7x10*  V  •  cm"',  /L  =  8x  10‘  cm"1. 

A,  of  the  p-MOSFET  is  within  the  range  of  published  data 
[10HI2K-4  x  107cm~'  in  [10],  2.25  x  10‘cnr1  in  [11], 
and  2  x  I06  cm  "  ’  in  {12]).  The  values  of  2?,  arc  close  to  those 
obtained  from  p-n  junction  experiments  (*3.45  x  10‘ 
V-cm*1  in  [10]  and  3.26  x  ^‘V^cm-1  in  [11]),  suggesting 
that  (2)  is  applicable  to  p-MOSFET’s  as  well. 

After  A,  and  B,  for  holes  are  extracted  from  measurements, 
/sub  of  the  p-MOSFET's  can  be  calculated.  A  physical 
expression  for  VD  WT  [3] 


D  SAT  = 


\Va-V,\EiaL« 

\Va-yl\-hESATLtn 


is  found  to  fit  the  experimental  V0  saT  very  well  and  is  used  in 
the  calculation  of  Em.  V,  is  the  extrapolated  threshold  voltage. 
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Fig.  1.  Constant  Im/lo  loci  and  VBur  loci  are  superimposed  on  the  Io- 
V0  characteristics  of  an  n-  and  a  p-MOSFET.  FD1AT  loci  are  determined 
with  the  method  described  in  [9],  T„  -  152  A. 


Fig.  2.  Plot  of  log  [f!v,V(/o(F0-/D5AT)l  versus  1/(FD-F0JAT)  to  show 
the  different  impact  ionization  rates  between  n-  and  p-MOSFETs.  Circle: 
|  V0\  -  6  V;  square  \V0\  -  4  V.  T„  -  152  A. 


I  a  in  the  saturation  region  is 

/©  =*  Vjat  WC«(  Vc~  y,~  Va  sat)  (5) 

where  i/SAT  (=*  5.6  x  106  cm/s)  is  the  saturation  velocity  of 
holes  obtained  by  fitting  the  measured  Id  with  (5)  for  different 
channel  lengths.  W  is  the  device  width  and  C„  is  the  gate 
capacitance.  The  effect  of  drain  voltage  on  V,  can  be 
incorporated  into  (5)  to  improve  the  accuracy  of  ID •  /sub 
calculated  by  this  simple  model  is  shown  in  Fig.  3.  For  buried- 
channel  p-MOSFET’s.  modification  of  this  model  is  expected 
since  Em  also  depends  on  the  depth  and  the  doping  concentra¬ 
tion  of  the  buried  layer. 

In  (3).  Bt  may  be  interpreted  as  the  ratio  of  the  impact 
ionization  energy  ♦,  to  the  carrier  mean  free  path  X.  is 
expected  to  be  very  roughly  equal  for  electrons  and  holes  (if 


one  considers  the  different  effective  masses  of  holes  and 
electrons,  holes  can  have  a  larger  ♦,  than  electrons  [13],  [14]). 
Given  this  assumption  and  X  =  92  A  for  electrons  [15],  we 
find  X  =  43  A  for  holes.  Notice  that  B,l  of  the  p-MOSFET’s 
is  higher  than  that  of  the  n-MOSFET’s  by  a  factor  of  2.2.  This 
means  the  p-MOSFET’s  can  take  twice  the  Vp-  VD  SAT  to 
generate  the  same  /Sub  (see  (1)— (3))  as  n-MOSFET’s.  As  a 
result,  the  hot-carrier-induced  breakdown  voltage  of  the  p- 
MOSFET  is  also  about  two  times  larger  as  shown  in  Fig.  4. 
The  hot-carrier- induced  breakdown  voltage  VBD  is  defined  as 
the  drain  voltage  at  which  the  drain  current  is  20  percent  larger 
than  in  the  saturation  region  as  illustrated  in  Fig.  1.  In  Fig.  4, 
notice  that  for  the  152-  and  the  90- A  devices,  the  VBD-Va 
curves  of  both  the  n-  and  p-MOSFET’s  show  the  C-shape 
characteristics  of  hot-carrier-induced  breakdown  [16],  while 
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for  the  432-  A  devices  the  C-shape  is  observed  only  in  the  n- 
MOSFET’s.  The  breakdown  voltage  of  the  432- A  p-MOS- 
FET's  is  about  20  V,  limited  by  the  p*-n  junction  breakdown. 

m.  Conclusions 

The  dependence  of  the  maximum  channel  field  on  terminal 
voltages  is  found  to  be  the  same  in  conventional  n-  and  p- 
MOSFET’s.  The  substrate  currents  of  p-  and  n-MOSFET's 
thus  can  be  characterized  with  the  same  method.  Substrate 
current  measurement  and  modeling  reveal  that  p-MOSFET’s 
can  withstand  twice  the  channel  field  of  n-MOSFET’s.  For 
small-geometry  p-MOSFET’s,  the  breakdown  voltage  also  is 
dominated  by  hot-carrier-induced  breakdown  as  in  n-MOS- 
FET's,  but  the  value  is  about  twice  as  large. 
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Subbreakdown  Drain  Leakage  Current  in  MOSFET 
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Abstract — Significant  drain  leakage  current  can  be  delected  al  drain 
voltages  much  lower  than  the  breakdown  voltage.  This  subbreakdown 
leakage  can  dominate  the  drain  leakage  current  at  zero  Vc  in  thin-oxide 
MOSFET's.  The  mechanism  is  shown  to  be  band-to-band  tunneling  in  Si 
in  the  drain/gate  overlap  region.  In  order  to  limit  the  leakage  current  to 
0.1  pA/pm.  the  oxide  field  in  the  gate-lo-drain  overlap  region  must  be 
limited  to  2.2  MV/cm.  This  may  set  another  constraint  for  oxide 
thickness  or  power  supply  voltage. 

I.  Introduction 

N  A  MOSFET  with  a  thin  gate  oxide  and  with  the  gate 
grounded,  significant  drain  leakage  current  can  be  observed 
at  drain  voltages  much  be!c»*’  what  is  usually  considered  the 
“breakdown  voltage"  [1].  We  shall  refer  to  this  as  the 
subbreakdown  leakage  current.  This  subbreakdown  leakage 
current  is  also  observable  in  the  subthreshold  1q-  Vc  character¬ 
istics  of  thin-oxide  devices.  Fig.  1  shows  the  typical 
subthreshold  characteristics  for  a  88- A  gate  oxide  n-MOS- 
FET.  The  large  leakage  current  near  zero  Vc  is  another 
manifestation  of  the  subbreakdown  leakage.  It  increases  with 
increasing  VD  and  decreasing  Vc.  A  good  understanding  and 
control  of  this  MOSFET  subbreakdown  leakage  current  is 
very  important  to  thin-oxide  VLSI  MOSFET's  and  some  high- 
voltage  devices  such  as  nonvolatile  memory  cells  or  program¬ 
ming  circuits. 

Deep-depletion  breakdown  in  MOS  capacitors  has  been 
studied  in  detail  [2].  In  [1],  Feng  et  al.  have  studied  the  effect 
of  scaling  oxide  thickness  and  channel  doping  concentration  on 
the  drain  avalanche  breakdown  voltage.  This  paper  explores 
the  mechanisms  responsible  for  the  MOSFET  drain  sub- 
breakdown  leakage  current. 

n.  Experiment 

Measurements  were  performed  on  n-channe!  MOSFET’s 
fabricated  with  conventional  n-type  polysilicon  gate  technol¬ 
ogy  The  p-bulk  was  doped  at  10 16  cm"3  and  the  source/drain 
diffusions  were  0.4-jim  As  junctions.  The  channel  was  20  pm 
long  and  20  nm  wide.  Gate  oxide  thicknesses  ranging  from  55 
to  350  A  were  examined. 

Fig.  2  shows  the  measured  drain  current  at  Vc  =  0  V  for 
several  oxide  thicknesses.  These  currents  flow  from  the  drain 
to  the  substrate.  For  the  155- A  gate  oxide  device,  data  are 
presented  for  three  different  temperatures,  25.  100.  and 
150*C.  As  temperature  increased,  a  ^-independent  leakage 
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Fig.  1_.  Subthreshold  characteristics  for  two  n-MOSFET's  with  T„,  « 
88  A,  W  -  50  jam.  and  L  «  4.}  and  0.6  jim.  Significant  drain  leakage 
current  car.  be  observed  when  Vqc  is  high 


Fig.  2.  Drain  current  characteristics  for  n-MOSFET's  at  different  tempera¬ 
tures  with  gate  grounded.  A  kVtndependent  leakage  current,  believed  to  be 
the  thermal  generation  current,  has  been  subtracted  for  clanry  The  drain 
current  characteristics  for  devices  with  T„  ■  55  and  350  A  at  room 
temperature  are  also  included. 

current  (2.28  x  10-’  at  100*C  and  1.17  x  1 0 ~ 7  at  150*0. 
believed  to  be  the  thermal  generation  currents  collected  by  the 
drain  junction,  appeared  as  expected.  The  data  shown  in  Fig.  2 
have  this  current  component  subtracted  for  clarity. 

ID.  Model  and  Discussion 

One  may  attempt  to  attribute  the  subbreakdown  leakage 
cun-ent  to  thermal  current  amplified  by  impact  ionization.  The 
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Fig  3  A  deep-depletion  region  is  formed  in  the  gate-to-drain  overlap 
region.  The  energy-band  diagram  illustrates  the  band-to-band  tunneling 
process  and  the  flow  of  earners.  Valence-band  electrons  tunnel  into  the 
conduction  band  and  are  collected  at  the  drain.  The  holes  created  flow  to  the 
substrate. 


data  at  different  temperatures  show  that  this  leakage  current  is 
insensitive  to  temperature  while  the  thermal  leakage  current 
varies  by  more  than  105  times  over  the  temperature  range. 
This  indicates  that  the  subbreakdown  drain  current  is  not 
amplified  thermal  leakage  current.  Another  possible  origin  of 
this  subbreakdown  current  is  the  amplification  of  the  electron 
tunneling  current  from  the  gate.  However,  the  drain  current  is 
more  than  three  orders  of  magnitude  larger  than  the  gate 
current  which  rules  out  this  possibility. 

We  interpret  the  subbreakdown  current  as  due  to  the  band- 
to-band  tunneling  process  in  the  gate-to-drain  overlap  region 
as  illustrated  in  Fig.  3.  Drain  current  is  due  to  the  tunneling  of 
valence-band  electrons  into  the  conduction  band.  The  holes 
created  by  the  tunneling  of  electrons  flow  to  the  substrate. 
Note  that  tunneling  is  only  possible  in  the  presence  of  a  high 
electric  field.  The  field  in  silicon  at  the  Si-SiO;  interface 
depends  on  the  doping  concentration  in  the  diffusion  region 
and  the  difference  between  VD  and  Va,  i.e.,  Vqc- 

Band-to-band  tunneling  current  density  is  the  highest  where 
the  electric  field  is  the  largest  and  the  band  bending  is  larger 
than  the  energy  bandgap  Er  A  simple  expression  for  the 
surface  electric  field  at  the  tunneling  point  in  the  gate-to-drain 
overlap  region  can  be  obtained  as  follows: 


where  E,  is  the  vertical  electric  field  at  silicon  surface.  3  is  the 
ratio  of  silicon  permittivity  to  oxide  permittivity .  and  T0I  is  the 
oxide  thickness  in  the  overlap  region.  A  band  bending  of  5.2 
eV  is  the  minimum  necessary  for  band-to-band  tunneling  to 
occur.  The  theory  of  tunneling  current  predicts  [3] 


T  irm*  >,2E}/1 1 

Id  =  AE,  exp - — -  -AE,  exp  (  -  B/E,)  (2) 

L  2\2qhE,  J 


where  A  is  a  preexponential  constant  and  B  =  21.3  MV/cm 
with  m*  =  0.2mo  [4],  According  to  (1),  E,  is  proportional  to 
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Fig.  4.  Plot  of  In  (/o/(F'oc  ~  1 .2))  versus  1  /{Vx  -  1.2).  where  Foe  is  the 
voltage  difference  between  the  drain  and  gate.  All  measured  data  with  the 
same  oxide  thickness  fall  on  a  straight  line,  which  is  in  agreement  with  the 
band-to-band  tunneling  model. 


Vqc  -  1.2.  Fig.  4  shows  that  measured  In  (/p/f^oc  ~  1.2)) 
plotted  against  1  /(k^c  ~  •  -2)  is  linear  in  agreement  with 
(2).  Measured  data  yield  an  experimental  value  of  B  =  18 
MV/cm.  This  is  in  reasonable  agreement  with  the  theoretical 
value  of  B  =  21.3  MV/cm.  which  has  never  been  verified  to 
this  accuracy.  Fig.  4  also  verifies  that  the  tunneling  current 
only  depends  on  the  vertical  field,  namely,  the  difference 
between  drain  and  gate  voltages  as  in  [1].  These  facts  strongly 
support  our  interpretation  that  the  currents  are  due  to  band-to- 
band  tunneling.  Using  measured  A  and  B,  a  drain  leakage 
current  of  ID  »  0.1  pA//im  can  be  created  when  E,  is  0.75 
MV/cm  and  the  oxide  field  is  about  2.2  MV/cm.  Empirically, 
we  have  found  ID  —  0.1  pA/jim  when  Voc  s  1.2  +  r,„  x 
£ml  with  £m,  ranging  from  2.2  to  3  MV/cm  for  a  wide  range 
of  MOSFET  technologies,  where  Ttoz  is  the  measured  gate 
(channel)  oxide  thickness.  The  fact  that  Ecm  is  larger  than  2.2 
MV/cm  suggests  that  the  oxide  near  the  gate  edge  is  thicker 
than  T,OI  due  to  polysilicon  gate  reoxidation.  This  sets  another 
limit  to  the  minimum  oxide  thickness  or  the  maximum  power 
supply  voltage. 

As  seen  in  Fig.  2,  there  is  a  noticeable  kink  in  the  drain  /-  V 
characteristics.  So  far  we  have  only  addressed  Region  1.  to  the 
left  of  the  kink.  In  Region  I.  lD  increases  slightly  with 
increasing  temperature  probably  due  to  a  decrease  in  Et  in  (2). 
In  Region  II.  lD  rises  above  the  predicted  tunneling  current. 
We  attribute  the  rise  to  the  amplification  of  the  tunneling 
current  by  avalanche  impact  ionization.  Id  in  Region  II  is 
characterized  by  a  negative  temperature  coefficient  due  to  the 
decrease  of  the  impact  ionization  coefficient  at  higher  temper¬ 
ature  and  by  the  sensitivity  to  the  substrate  bias  (not  presented 
here),  of  which  ID  in  Region  I  is  independent. 

IV.  Summary 

Subbreakdown  leakage  is  the  significant  drain  leakage 
current  observable  at  beiow-breakdown  voltage.  The  domi- 
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nant  mechanism  for  this  current  is  band-to-band  tunneling  at 
the  Si-SiO;  interface  in  the  gate-to-drain  overlap  region.  At 
higher  drain  voltages,  impact  ionization  amplifies  this  current 
and  results  in  even  larger  leakage  currents  until  breakdown 
occurs. 

In  thin  oxide  devices,  the  OFF-state  drain  current  often 
increases  with  decreasing  Vc  measurably  due  to  the  same 
mechanism  (Fig.  1).  The  voitage  required  to  cause  band-to- 
band  tunneling  will  decrease  when  oxide  thickness  decreases. 
This  may  be  another  barrier  to  VLSI  MOSFET  scaling. 
Possible  solutions  to  the  problem  are  the  use  of  lightly-doped 
drain  devices  which  have  lower  or  graded  gate  oxide 
devices  [5]  in  which  the  oxide  is  thicker  in  the  gate-to-drain 
overlap  region  than  that  in  the  channel. 
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Abstract 

Significant  gate-induced  drain  leakage  current  can  be 
detected  in  thin  gate  oxide  MOSFETs  at  drain  voltages  much 
lower  than  the  juncrion  breakdown  voltage.  This  current  is 
found  to  be  due  to  the  band-to-band  tunneling  occurring  in  the 
deep-depletion  layer  in  the  gate-to-drain  overlap  region.  In 
order  to  limit  the  leakage  current  to  O.lpA/pm,  the  oxide  field 
in  the  gate-to-drain  overlap  region  must  be  limited  to 
1.9MV/cm.  This  may  set  another  constraint  for  the  power  sup¬ 
ply  voltage  and/or  oxide  thickness  in  VLSI  MOSFET  scaling. 
Device  design  considerations  for  minimizing  the  gate-induced 
drain  leakage  current  are  discussed. 


I.  Introduction 

As  MOS  technology  is  scaled  down  to  the  sub-half 
micron  regime,  it  is  anticipated  that  the  power  supply  voltage 
will  need  to  be  reduced  to  below  5  volts:  Hot-electron  reliabil¬ 
ity,  oxide  reliability,  and  punch-through  are  widely  discussed 
reasons!  1],  This  paper  repons  a  new  mechanism  -  Gate- 
Induced  Drain  Leakage  Current  -  that  imposes  an  additional 
constraint  on  the  acceptable  supply  voltage  and/or  oxide  thick¬ 
ness.  Significant  drain  leakage  current  can  be  detected  at  drain 
voltage  much  lower  than  the  breakdown  voltage.  This  drain 
leakage  current  is  caused  by  the  gate-induced  high  electric 
field  in  the  gate-to-drain  overlap  region  and  can  dominate  the 
drain  leakage  current  at  zero  gate  bias  in  thin  oxide  MOS¬ 
FETs.  A  model  based  on  band-to-band  tunneling  is  proposed 
and  can  satisfactorily  explain  the  experimental  results. 

Its  dependence  on  fabrication  process  and  device  param¬ 
eters  has  been  carefully  studied.  The  exact  oxide  thickness  and 
doping  profiles  in  the  gate-to-drain  overlap  region  are  found  to 
play  important  roles  in  the  gate-induced  drain  leakage  current. 
The  guidelines  for  designing  thin  oxide  VLSI  MOSFETs  and 
high  voltage  devices  such  as  those  used  in  non-volatile 
memory  circuits  are  discussed. 

IL  Gate-Induced  Drain  Leakage  Current  Characteristics 

A  typical  subthreshold  Iq-Vg  characteristics  for  88-/? 
gate  oxide  n-MOSFETs  wtth  two  different  channel  lengths  are 
shown  in  Fig.  1.  Significant  amount  of  drain  leakage  current 
near  zero  VG  is  observed.  This  drain  leakage  current  increases 


with  increasing  Vq  and  decreasing  VG,  and  shows  no  apparent 
dependence  on  channel  length.  Fig.  2  shows  the  off-state  (VG 
a  0V)  drain  leakage  currents  for  several  oxide  thicknesses. 
These  currents  flow  from  the  drain  to  the  substrate.  For  the 
155-:?  gate  oxide  device,  data  are  presented  for  three  different 
temperatures  -  25,  100.  and  150°C.  A  Vo-independent  leakage 
current,  believed  to  be  the  thermal  generation  current,  has 
been  subtracted  from  the  high  temperature  curves  for  clarity. 


IIL  Leakage  Mechanism 

One  may  attempt  to  attribute  the  off-state  drain  leakage 
current  to  thermal  current  amplified  by  impact  ionization.  The 
insensitivity  of  the  leakage  current  to  temperature  shown  in 
Fig.  2  indicates  that  the  drain  leakage  current  is  not  amplified 
thermal  leakage  current.  Another  possible  origin  of  this  drain 
leakage  current  is  the  amplification  of  the  electron  tunneling 
current  from  the  gate.  However,  the  drain  current  is  more 
than  three  orders  of  magnitude  higher  than  the  gate  currant 
which  rules  out  this  possibility. 

We  attribute  the  gate-induced  drain  leakage  current  to 
the  band-to-band  tunneling  process  in  silicon  in  the  gate-to- 
drain  overlap  region  as  illustrated  in  Fig.  3.  The  cross-section 
shown  in  Fig.  3  is  simply  a  gated-diode  configuration.  When 
high  voltage  is  applied  to  the  drain  with  the  gate  grounded,  a 
deep-depletion  region  is  formed  underneath  the  gate-to-drain 
overlap  region.  Electron-hole  pairs  are  generated  by  the  tun¬ 
neling  of  valence  band  electrons  into  the  conduction  band  and 
collected  by  the  drain  and  substrate,  separately.  Since  all  the 
minority  carriers  generated  thermally  or  by  band-to-band  tun¬ 
neling  in  the  drain  region  (holes  in  n-MOSFETs  and  electrons 
in  p-MOSFETs)  flow  to  the  substrate  due  to  the  lateral  field, 
the  deep-depletion  region  is  always  present  and  the  band-to- 
band  tunneling  process  can  continue  without  creating  an 
inversion  layer. 

Band-to-band  tunneling  is  only  possible  in  the  presence 
of  a  high  electric  field  and  when  the  band  bending  is  larger 
than  the  energy  band  gap,  E,.  The  field  in  silicon  at  the  Si- 
SiOj  interface  also  depends  on  the  doping  concentration  in  the 
diffusion  region  and  the  difference  between  VD  and  VG,  i.e. 
Vjxj.  A  simple  expression  for  the  surface  field  at  the  dom¬ 
inant  tunneling  point  can  be  expressed  as 


*  Present  address  :  Intel  Corporation,  3065  Bowers  Ave¬ 
nue,  Santa  Clara,  CA.  95051. 


where  E,  is  the  vertical  electrical  field  at  silicon  surface.  3  is 
the  ratio  of  silicon  permittivity  to  oxide  permittivity,  and  T0, 
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is  the  oxide  thickness  in  the  overlap  region.  A  band  bending  of 
1.2  V  is  the  minimum  necessary  for  tunneling  process  to 
occur.  The  theory  of  tunneling  current  predicts  [2] 

Id  =  AE,exp(-B/E»)  (2) 

where  A  is  a  pre-exponential  constant  and  B  *  21.3  MV/cm 
[3J.  Fig.  4  shows  the  measured  ln(Iiy'(VDG-1.2)]  plotted 
against  1/(Vdc_1-2).  All  measured  data  with  the  same  oxide 
thickness  fall  on  straight  lines,  in  agreement  with  the  tunnel¬ 
ing  theory  (Eq.  (2)).  Fig.  4  also  verifies  that  the  tunneling 
current  only  depends  on  the  vertical  field,  namely,  the  voltage 
difference  between  the  drain  and  gate.  The  slope  of  the  lines 
in  Fig.  4  is  3T„B  (Eqs.  (1)  and  (2))  and  is  plotted  in  Fig.  5. 
An  experimental  value  of  B  »  18  MV/cm  was  obtained,  in  rea¬ 
sonable  agreement  with  the  theoretical  value  of  B  =  21.3 
MV/cm  [4],  These  facts  strongly  support  our  interpretation 
that  the  currents  are  due  to  the  band-to-band  tunneling. 

At  high  Vdg.  Id  was  found  to  rise  above  the  tunneling 
current  predicted  by  Eq.  (2)  as  shown  in  Fig.  6.  We  attribute 
the  rise  to  the  amplification  of  the  tunneling  current  by 
avalanche  impact  ionization  which  accounts  for  the  negative 
temperature  coefficient  due  to  the  decrease  of  the  impact  ioni¬ 
zation  coefficient  at  higher  temperature  (Region  II  as  indi¬ 
cated  in  Fig.  2)  and  the  sensitivity  to  the  substrate  bias  (Fig. 
7),  of  which  Id  at  low  Vdg  (Region  I  as  indicated  in  Figs.  2 
and  7)  is  independent.  The  kink  observed  in  n-MOSFETs 
shown  in  Fig.  6  is  due  to  hole  trapping  in  the  oxide  which 
affects  the  increase  of  the  surface  field. 

IV.  Impact  on  MOSFET  Scaling  and  Design 

To  limit  the  undesirable  gate-induced  drain  leakage 
current  to  O.lpA  per  pm  channel  width,  the  oxide  field  in  the 
gate-to-drain  overlap  region  must  be  limited  to  1.9MV/cm  as 
can  be  obtained  in  Fig.  8.  This  sets  another  limit  to  the 
minimum  oxide  thickness  or  the  maximum  power  supply  vol¬ 
tage  in  MOSFET  scaling, 

Vcc  =1.2  +  T„,x  1 ,9MV/cm  (3) 

With  this  constraint,  a  100-i£  device  can  only  tolerate  a  max¬ 
imum  supply  voltage  of  3.1V.  The  measured  Vdg  for  Id  » 
0.  Ip  A/pm  is  usually  larger  than  predicted  by  Eq.  (3)  depend¬ 
ing  on  the  process  technologies  as  demonstrated  in  Fig.  8. 
This  is  due  to  the  larger  oxide  thickness  near  the  gate  edge 
resulted  from  source-drain  re-oxidation.  Since  the  gate- 
induced  drain  leakage  current  depends  on  the  oxide  field,  the 
doping  type  of  the  polvsilicon  gate  used  in  thin  oxide  MOS¬ 
FETs  strongly  affects  the  acceptable  supply  voltage.  The 
acceptable  Vcc  for  O.lpA/pm  is  expected  to  be  smaller  by 
about  IV  if  the  poly-gate  and  the  dram  region  have  opposite 
doping  types,  e.g.  p-MOSFET  with  n-type  poly-gate,  because 
of  the  work  function  difference.  This  gives  another  reason  to 
use  p*  poly-gate  in  sub-half  micron  p-MOSFETs.  in  addition 
to  the  concern  for  making  surface-channel  devices 

Possible  means  to  reduce  this  gate -induced  drain  leakage 
current  are  the  use  of  lightly-doped  <  LDD  dev  ices  or  graded 
gate  oxide  devices  [5],  However,  the  effectiveness  of  using 
LDD  device  strongly  depends  on  the  d'r "<  .■•.■■ncentraoon  in 
the  n'  region.  Fig.  9  shows  the  caic-.ated  • needed  to 
create  an  oxide  field  of  l.9MV/cm  :.n  rt  gi  t  dram  overlap 
region  and  induce  1.2V  band  bending  ni'.n  i-  'f  the  dop¬ 
ing  concentration.  A  minimum  is  otv.a  -ed  l-  -.-J  '. O-’cm-3. 
This  suggests  that  LDD  MOSFETs  ~-m  *a-r  i  v  doping 
lower  than  10l8cm-3  to  raise  the  accept- -t  -  ,  vet  that  of 


non-LDD  MOSFETs.  LDD  MOSFETs  usually  have  a 
moderately  doped  n*  region  in  order  to  reduce  the  series  resis¬ 
tance  and  minimize  the  hot-electron  degradation.  If  the  doping 
concentration  is  higher  than  1018cm"J  at  the  gate  edge,  a  point 
with  doping  concentration  of  1018cm"3  always  exists  in  the 
gate/drain  overlap  region  ,  which  makes  LDD  MOSFETs  no 
better  than  conventional  MOSFETs  for  the  gate-induced  drain 
leakage  current.  Buried  LDD  MOSFETs  with  the  peak  n~  dop- 
in-  concentration  several  hundred  angstroms  underneath  the 
S1-S1O2  interface  seems  to  be  the  more  favorable  device  struc¬ 
ture  for  both  hot-electron  reliability  [6]  and  gate-induced  drain 
leakage  current. 

Since  the  T0»  to  be  used  in  Eq.  (3)  is  that  present  over  the 
point  where  the  drain  doping  concentration  is  10l8cm"3, 
graded  gate  oxide  device  structure  could  effectively  reduce 
the  leakage  current.  However,  the  weakened  gate  control  near 
the  gate  edge  might  become  a  serious  concern  in  device  relia¬ 
bility  [7],  Judicious  design  of  the  gate  edge  profile  to 
compromise  device  reliability  and  gate-induced  drain  leakage 
current  is  necessary  if  the  graded  gate  oxide  device  structure  is 
adopted.  The  gate-induced  drain  leakage  current  in  n- 
MOSFETs  can  also  be  significantly  enhanced  after  hot- 
electron  stressing  due  to  electron  trapping  in  the  oxide  and  at 
the  interface  as  shown  in  Fig.  10,  while  p-MOSFETs  show 
continuous  decrease  in  the  drain  leakage  current  because  the 
dominant  trapped  charges  are  still  electrons.  This  may  pose  a 
long-term  device  drift  issue  in  n-channel  devices. 

V.  Conclusion 

The  gate-induced  drain  leakage  current  will  inevitably 
become  an  important  issue  in  device  design  and  fabrication 
when  thinner  gate  oxides  are  used  in  sub-micron  and  deep 
sub-micron  MOSFETs.  Tne  voltage  required  to  cause  band- 
to-band  tunneling  decreases  when  oxide  thickness  decreases. 
Since  the  band-to-band  tunneling  process  is  caused  by  the 
high  oxide  field  in  the  gaie-to-drain  overlap  region,  the  exact 
oxide  thickness  and  doping  profiles  in  the  gate/drain  overlap 
region  require  great  care  in  current  and  future  device  design 
and  fabrication. 

VI.  Acknowledgement 

This  work  is  supported  by  AMD,  Hughes  and  Rockwell 
International  under  State  of  California  MICRO  program  and 
JSEP  under  contract  F49620-84-C-0057. 

References 

[1]  Masakazu  Kakumu.  Masaaki  Kinugawa,  Kazuhiko 
Hashimoto  and  Jun'ichi  Matsunaga,  1986  IEDM  Techni¬ 
cal  Digest,  pp.  399-402. 

[2]  J.L.  Moll,  "Physics  of  Semiconductor,"  McGraw  Hill,  P. 
253,  1964. 

[3]  S.M.  Sze,  "Physcis  of  Semiconductor  Devices,"  2nd  edi¬ 
tion,  Wiley  New  York,  p.  525, 1981. 

[4]  J.  Chen,  T.Y.  Chan.  I.C.  Chen,  P.K.  Ko.  and  C.  Hu.  IEEE 
Electron  Device  Lett..  November  1987. 

[5]  P.K.  Ko.  S.  Tam  and  C.  Hu,  1984  IEDM  Technical  Dig¬ 
est,  p.  88. 

[6]  Hugh  R.  Grinolds,  Masaaki  Kinugawa,  and  Masakazu 
Kakumu.  1985  IEDM  Technical  Digest,  p.  246. 

[7]  T.Y.  Chan.  A.T.  Wu.  PX  Ko,  and  C.  Hu,  IEEE  Electron 
Device  Lett.,  vol.  EDL-8,  p.  326,  July  1987. 


31.3 


IEDM  87-719 


1 


Fig.  1  Subthreshold  characteristics  of  thin  oxide  n-MOSFETs 
with  T„  =*  88/? ,  and  Ltff  *  0.6  and  4.5  fim.  Significant  off- 
state  (Vc  -  0V)  drain  leakage  current  can  be  observed  when 
VD  is  high. 

VG  *  electro" 

•  hole 


Fig.  2  Drain  current  characteristics  for  n-MOSFETs  with  gate 
grounded.  The  drain  leakage  currents  show  strong  dependence 
on  the  oxide  thickness,  but  are  insensitive  to  temperature.  A 
Vo -independent  leakage  current,  believed  to  be  the  thermal 
generation  current,  has  been  subtracted  for  clarity. 
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Fig.  3  (a)  A  deep-depletion  region  is  formed  in  the  gate/drain  overlap  region,  (b)  The  en¬ 
ergy  band  diagram  illustrates  the  band-to-band  tunneling  process  in  silicon  in  the 
gate/drain  overlap  region.  Electron-hole  pairs  are  generated  by  the  tunneling  of  valence 
band  electrons  into  the  conduction  band  and  collected  by  the  drain  and  substrate, 
separately. 
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Fig.  4  Plot  of  ln(fo/(Vo -1.2)|  versus  :  2'.  where 

Vqq  is  the  voltage  difference  between  •.re  t/j  -  i.-d  gate  All 
measured  data  with  the  same  oxide  thick-esi  ‘i..  -n  straight 
lines,  which  is  in  agreement  with  the  bond  >  and  unnehng 
theory. 
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Fig.  5  The  slope,  3 T„  B.  of  the  lines  in  Fig.  4  as  a  funcdon  of 
the  oxide  thickness  for  both  n-  and  p-MOSFETs.  An  experi¬ 
mental  value  of  B  *  18  MV/cm  was  obtained,  in  agreement 
with  the  theoretical  value  of  B  *  21.3  MV/cm. 
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Fig.  6  lo~vDC  characteristics  from  measurement  and 
predicted  by  the  simple  model.  The  rise  of  measured  Ip  above 
the  predicted  tunneling  current  is  attributed  to  the 
amplification  of  the  tunneling  current  by  avalanche  impact 
ionization.  The  kink  observed  in  n-MOSFETs  is  due  to  hole 
trapping  in  the  oxide. 


(a) 


(b) 


Fig.  8  Measured  Vqc  for  O.lpA/pim  drain  leakage  current  as 
a  function  of  oxide  thickness  for  both  n-  and  p-MOSFETs.  It 
is  found  that  Ip  =  O.lpA/nm  when  VpQ  *  1.2  +  Tm  x 
1.9MV/cm.  The  higher  Vdg  for  some  process  technologies  is 
due  to  the  larger  oxide  thickness  near  the  gate  edge  resulted 
from  source-drain  re-oxidation. 


Fig.  9  VDC  needed  to  create  an  oxide  field  of  1.9MV/cm  in 
the  gate/drain  overlap  region  and  to  induce  1.2V  band  bending 
as  a  function  of  the  doping  concentration.  A  minimum  is  ob¬ 
served  around  1018c/n-3. 


Fig.  7  Iq-Vqq  characteristics  at  different  substrate  bias  for 
(a)  n-MOSFET  and  (b)  p-MOSFET.  The  senstuvity  of  the 
drain  leakage  current  to  the  substrate  bias  at  high  Vdg  indi¬ 
cates  the  amplification  of  the  tunneling  current  by  avalanche 
impact  ionization. 


Fig.  10  The  gate-induced  drain  leakage  current  for  a  n- 
MOSFET  measured  at  Vc  *  0V  and  VD  »  5V  as  a  function  of 
the  stressing  time.  The  leakage  current  is  significantly 
enhanced  after  hot-electron  stressing  due  to  election  trapping 
in  the  oxide. 
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Ultra-Thin  Silicon-Dioxide  Breakdown 
Characteristics  of  MOS  Devices  with 
n+  and  p+  Polysilicon  Gates 
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Abstract— Is  this  work  we  investigate  the  effect  of  the  gate  material  on 
the  breakdown  characteristics  of  ultra-thin  silicon  dioxide  films  at  low 
voltages  (<4  V).  When  MOS  capacitors  are  stressed  with  a  positive  gate 
voltage,  the  charge  to  breakdown  and  time  to  breakdown  at  a  fixed  oxide- 
voilage  drop  are  significantly  smaller  in  p*  polysilicon-gale  capacitors 
than  in  n'  polysilicon-gale  capacitors.  The  results  are  interpreted  in  terms 
of  a  simple  model  of  hole  tunneling  resulting  from  hot-hole  generation  in 
the  anode  by  hot  electrons  entering  from  the  silicon  dioxide.  Extrapola¬ 
tion  of  high-voltage-breakdown  lifetime  measurements  for  relatively 
thick-oxide  devices  to  low  voltages  may  be  complicated  by  this  mecha¬ 
nism. 

I.  Introduction 

REAKDOWN  in  thin  films  of  silicon  dioxide  has  received 
considerable  attention  in  the  recent  technical  literature 
[l]-[5].  In  these  previous  studies  the  voltage  applied  during 
the  breakdown  measurement  is  usually  greater  than  10  V, 
which  is  much  larger  than  the  actual  power-supply  voltage. 
Hence  the  study  of  silicon-dioxide  breakdown  in  the  low- 
voltage  regime  ( *  5  V)  is  of  much  interest. 

The  investigation  of  low-voltage  breakdown  in  a  conven¬ 
iently  short  time  requires  the  use  of  utlra-thin  oxides  (thinner 
than  »  5  nm).  For  oxide  thicknesses  in  this  range  it  has  been 
reported  that  the  charge  to  breakdown  QBD  (the  total  charge 
injected  into  the  oxide  necessary  to  cause  breakdown)  and  time 
to  breakdown  tBD  are  strongly  dependent  on  the  stressing 
electric  field  [3),  [6],  more  so  than  in  thicker  oxides.  In 
addition,  ultra-thin  silicon  dioxide  films  are  being  considered 
for  subhalf-micrometer  MOS  transistors  [7],  [8]  and  therefore 
the  study  of  breakdown  in  these  films  is  of  practical 
importance. 

II.  Sample  Preparation 

Poly  silicon-gate  MOS  capacitors  were  fabricated  for  this 
study.  The  starting  substrates  were  4-in-diameter  6-12-Q'cm 
n-type  and  15-30-0  cm  p-type.  Gate  oxides  were  grown  at 
800"C  after  patterning  and  etching  of  a  steam-grown  field 
oxide.  Before  gate  oxidation  the  furnace  tube  was  cleaned  for 
8  h  at  1050*C  in  an  (tytnchlorethane  (TCA)  ambient. 
Polysilicon  was  deposited  undoped  at  650*C  immediately  after 
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gate  oxidation.  For  each  wafer  the  polysilicon  was  doped  p* 
on  one  half  and  n*  on  the  other  half  by  using  a  CVD  oxide  as  a 
diffusion  mask.  The  n*  poiysiiicon  was  doped  from  a  POCly 
source  at  950*C  and  the  p*  polysilicon  was  doped  from  a  solid 
source  containing  B;Oj  at  900*C. 

III.  Experimental  Results  and  Discussion 

Fig.  1  shows  the  measured  tBD  and  QBD  versus  the  oxide- 
voltage  drop  yax  during  stress  for  MOS  capacitors  of  area  1  x 
10‘4  cm2  and  thickness  3.3  nm.  Data  are  shown  for  both  p* 
and  n'  polysilicon-gate  devices.  The  solid  curve  is  a  theoreti¬ 
cal  plot  that  will  be  discussed  later.  The  gate  voltage  during 
stress  was  held  constant  (positive  with  respect  to  the  sub¬ 
strate).  The  tunneling  current  during  stress  remained  essen¬ 
tially  constant  until  breakdown  occurred.  In  some  cases  the 
tunneling  current  was  observed  to  increase  slightly  during  the 
measurement  indicating  the  possibility  of  positive  charge 
trapping.  Breakdown  was  detected  as  a  sudden  large  increase 
in  the  current  through  the  device.  V„  is  the  actual  voltage  drop 
across  the  silicon  dioxide,  i.e..  the  applied  voltage  minus  the 
work-function  difference  between  the  polysilicon  and  n-type 
substrate  minus  any  band  bending.  For  the  p*  polysilicon-gate 
devices  a  work-function  difference  plus  band  bending  of  1 . 1  V 
was  assumed.  We  found  that  for  the  n*  polysilicon-gate 
devices  a  work-function  difference  plus  band  bending  of  0.4  V 
was  necessary  in  order  to  have  equal  tunneling  currents  on  a 
gate  current  versus  V„  plot  (equal  current  implies  equal  Vox). 
Similar  experiments  on  magnesium-  and  aluminum-gate  MOS 
capacitors  yielded  QBd  values  on  the  order  of  0.1  C/cm:  [9], 
and  hence  metal-gate  devices  are  not  suitable  for  this  type  of 
study. 

As  can  be  seen  in  Fig.  1 .  there  is  a  reduction  in  tBD  and  QB0 
for  the  p*  polysilicon-gate  devices  of  approximately  three 
orders  of  magnitude  when  compared  to  the  n*  devices.  For 
this  experiment  the  cathode  (the  negatively  biased  electrode) 
during  stress  is  the  n-type  substrate  and  hence  is  the  same  for 
both  types  of  devices.  The  anode  during  stress  (the  polysilicon 
gate)  is  not  the  same,  however. 

Fig.  2  shows  the  results  of  a  similar  experiment  except  that 
in  this  case  the  gate  is  biased  negative  with  respect  to  the 
substrate.  The  substrate  is  p-type  and  the  capacitor  area  is  1  x 
10" 4  cm2.  The  solid  curve  is  a  theoretical  plot.  A  constant  gate 
voltage  was  applied  and  as  before  the  tunneling  current  during 
stress  was  essentially  constant.  For  this  case,  the  tBD  at  a  fixed 
current  density  (hence  fixed  oxide-voltage  drop)  is  indepen- 
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Fig  I.  (a)  Time  10  breakdown  and  (b)  charge  to  breakdown  Q,D  versus 
oxide-volugc  drop  comparing  fl*  and  p*  poly  silicon  gate  MOS  capacitors 
(oxide  thickness  »  3.3  nm,  capacitor  area  »  1  x  10'*  cm:).  The  gate 
voltage  was  constant  and  positive  with  respect  to  the  substrate.  All  devices 
are  from  the  same  wafer.  The  solid  curve  is  the  theoretical  model. 

dent  of  gate  type.  In  this  experiment  the  anode  during  stress 
(the  p-type  substrate)  is  the  same  for  both  types  of  devices. 

Since  the  breakdown  characteristics  of  Figs.  1  and  2  are 
compared  at  equal  oxide-voltage  drops,  the  energy  gained  by 
the  tunneling  electrons  traversing  the  silicon  dioxide  should  be 
the  same  for  both  types  of  devices.  Hence  it  is  unlikely  that  the 
breakdown  is  a  result  of  hot -electron  effects  in  the  silicon 
dioxide  itself.  Also,  the  maximum  potential  energy  that  can  be 
gained  by  the  tunneling  electrons  in  the  silicon  dioxide  is  q(  Vox 
-  0,),  where  <t>„  is  the  barrier  height  between  the  silicon  and 
silicon  dioxide.  For  this  set  of  experiments  the  maximum 
energy  gained  in  the  silicon  dioxide  is  approximately  1-2  V, 
which  is  probably  too  small  to  cause  impact  ionization  in  the 
silicon  dioxide. 

However,  the  electrons  will  enter  the  anode  electrode  with 
q  Kw  electronvolts  of  energy  relative  to  the  conduction  band  in 
the  anode  if  scattering  in  the  silicon  dioxide  is  negligible  as 
shown  in  Fig.  3.  Since  this  energy  can  be  substantially  greater 
than  the  bandgap  of  silicon,  it  is  possible  to  generate  hole- 
electron  pairs  in  the  anode  electrode  [10].  It  is  thus  possible  to 
generate  hot  holes,  a  fraction  of  which  can  tunnel  into  the 
silicon  dioxide  as  originally  proposed  by  Weinberg  et  al.  ( 1 1 J. 
These  holes  will  tend  to  increase  the  cathode  electric  field 
which  in  turn  could  lead  to  breakdown.  In  the  n  *  -gate  case  the 


J  mA/cm* 

Fig.  2.  Time  lo  breakdown  l,D  versus  stressing  current  density  J  comparing 
n*  and  p*  poly  silicon  gate  MOS  capacitors  (oxide  thickness  ■  4  nm. 
capacitor  area  »  I  x  10'*  cm:).  The  gate  voltage  was  constant  and 
negative  with  respect  to  the  substrate.  All  devices  are  from  the  same  wafer 
The  solid  curve  is  the  theoretical  model. 

valence-band  electrons  must  be  excited  to  the  bottom  of  the 
conduction  band.  However,  in  the  p*  case  the  valence-band 
electrons  only  need  to  excited  to  the  top  of  the  valence  band 
since  there  are  available  states  there.  Thus,  for  the  p*  case  the 
barrier  height  for  hole  tunneling  will  be  1.1  eV  (the  silicon 
bandgap)  lower  than  in  the  n*  case  as  shown  in  Fig.  3. 

If  the  fraction  of  tunneling  electrons  which  give  up  all  of 
their  energy  by  creating  hole-electron  pairs  is  aH.  then  the 
hole -tunneling  current  arising  from  hole  generation  in  the 
anode  will  be  given  by  Jp  =  aHJ„  &H  where  J„  is  the  electron 
tunneling  current  from  the  cathode  and  9W  is  the  hole¬ 
tunneling  probability.  The  charge  associated  with  the  trapped 
holes  at  time  t  is  thus  Qp(t )  =  r\Jpt  where  r>  is  the  hole 
trapping  efficiency.  Assuming  that  breakdown  occurs  when  a 
critical  number  of  holes  have  tunneled  into  the  SiO;  and  have 
been  trapped  [5],  then  QBD  will  be  given  by 

Qbd  =  JntBD  =  CQ  - 1  =  C  exp  ^  [0p{  V0X)]  j  (1 ) 

<t>p(Vox)  =  f’g.s.O;-  Q<t>n-  p*  polysilicon 

^p(voi)~  EttSxOl-q0„-  qVox  n*  polysilicon 

where  C  is  a  constant.  Xox  is  the  oxide  thickness.  B'  =  Bi 
<t>3n,:  is  the  Fowler-Nordheim  slope  parameter.  £,s,0.  and  £,.s, 
are  the  bandgap  energies  of  silicon  dioxide  and  silicon, 
respectively,  and  <t>„  and  <t>p(  Vox)  are  the  barrier  heights  for 
electrons  and  holes,  respectively.  (Equation  (1)  is  derived 
assuming  Fowler-Nordheim  tunneling  and  equal  electron  and 


<bl 

Fig  3  Energy-band  diagram  depicting  the  mechanism  of  hole  tunneling 
due  10  hole  generation  in  the  polysilicon  gate  (the  anode  for  the  case 
showni.  The  gate  is  biased  positive  with  respect  to  the  n-type  substrate.  «, 
is  the  barrier  height  for  hole  tunneling,  (a)  p"  polysilicon  gate,  <bi  n ' 
polysilicon  gate. 

hole  effective  masses.  In  addition,  the  field  dependence  in  the 
pre-exponential  factor  of  the  Fowler-Nordheim  expression 
has  been  neglected  for  simplicity,  i.e..  7„  »  A  exp  [  -  BXot/ 
V„\  where  A  is  a  material  constant.)  The  dependence  of  tBo 
on  Vox  can  be  determined  by  dividing  (1)  by  7,. 

The  solid  curves  shown  in  Figs.  I  and  2  are  the  theoretical 
predictions.  For  Fig.  1  the  pre-exponential  factor  C  was 
arbitrarily  fit  at  one  point  for  the  n*  data  and  the  same  pre- 
exponenttai  factor  was  used  to  generate  the  theoretical  curve 
for  the  p  ‘  data.  The  model  parameters  used  were  9  eV  for  the 
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Fig.  4.  Time  to  breakdown  t,D  versus  the  oxide  voltage  drop  for  oxide 
thicknesses  4.3.  6.3.  and  9  3  nm.  The  solid  curve  is  the  theoretical  plot  for 
the  4.3-nm  devices.  The  gate  voltage  was  constant  and  positive  with  respect 
to  the  substrate. 

silicon-dioxide  bandgap,  3.2  eV  for  the  barrier  height  for 
electron  tunneling,  and  =  250  MV/cm  [12].  As  is 
evident  from  Fig.  I,  both  the  dependence  on  Vox  and  the 
difference  in  magnitude  resulting  from  p*  and  n*  gate  doping 
are  well  represented  by  the  simple  model.  In  Fig.  2  V„  was 
determined  from  the  measured  Jn~Vox  characteristics. 

Fig.  4  compares  the  effect  of  both  the  gate  type  and  oxide 
thickness  on  the  time  to  breakdown.  Theoretical  plots  are 
shown  for  the  4.3-nm  oxide.  However,  the  6.3-  and  9.5-nm 
oxides  will  have  significant  electron  scattering  in  the  oxide  and 
hence  the  simple  model  presented  does  not  apply  directly.  In 
addition,  for  the  thicker  oxides  the  difference  in  iBD  due  to  gate 
type  is  not  as  large  as  in  the  4.3-nm  case,  and  the  difference 
becomes  smaller  at  higher  voltages  suggesting  that  a  different 
mechanism  of  hole  generation  may  be  operable  at  relatively 
high  voltages.  If  so,  a  simple  extrapolation  of  high-voltage  tB0 
measurements  may  not  give  a  valid  estimation  of  the  device 
lifetime  at  low  voltages. 

IV.  Summary 

In  conclusion,  we  have  presented  a  simple  model  for  the 
breakdown  of  ultra-thin  silicon-dioxide  films  at  low  voltages. 
The  difference  between  p~  and  n*  polysilicon-gate  tBD  and 
Qbd  can  be  explained  if  one  assumes  that  at  low  voltages  holes 
are  generated  in  the  anode  and  that  some  fraction  of  these 
holes  then  tunnel  into  the  oxide.  In  this  paper  the  discussion 
has  been  in  terms  of  V„.  When  one  compares  the  results  at  the 
same  applied  voltage  the  differences  due  to  gate  material  are 
not  as  significant.  Hence  we  do  not  believe  that  this  effect 
should  determine  the  choice  of  n*  and  p*  gate  technology. 
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ABSTRACT 

A  study  of  the  performance  and  hot-electron 
reliability  of  submicron  n-channel  MOSFET’s  is 
presented.  Well-established  hot-electron-based 
physical  models  are  adequate  in  explaining  the 
general  behaviors  of  the  drain,  substrate,  and  gate 
currents  of  these  devices.  These  results  suggest 
that  the  basic  physics  is  rather  well-understood  and 
the  design  criteria  developed  for  micron-size  dev¬ 
ices  can  be  extended  to  cover  their  deep-submicron 
counterparts.  Hot-electron  studies  reveal  a 
channel-length  dependence  in  device  degradation. 
This  phenomenon  together  with  gate-induced  drain 
leakage  current  (1)  will  impose  an  upper  limit  on 
the  supply  voltage  and  a  lower  limit  on  the  gate 
oxide  thickness.  Based  on  device  degradation 
results  alone,  the  power  supply  voltage  for  a 
quarter-micron  device  with  oxide  thickness  of  86 
A  should  be  limited  to  2.5  V  if  no  degradation- 
resistant  structure  is  used. 


INTRODUCTION 

Performance  and  hot-electron  reliability  of  sub¬ 
micron  MOSFET’s  have  been  the  two  major  con¬ 
cerns  for  their  feasibility  in  circuit  applications. 
Recent  years  have  seen  increased  activities  in  the 
studies  of  deep-submicron  MOSFET’s  due  to  their 
potential  applications  in  high-performance  VLSI 
circuits  and  systems  (2I-(4J.  However,  the 
emphases  of  most  previous  reports  have  been  to 
demonstrate  the  feasibility  of  fabricating  these  dev¬ 
ices.  This  paper  attempts  to  construct  a  unified 
understanding  of  their  characteristics,  thereby 
establishing  some  basic  device  design  guidelines. 

The  threshold  voltage,  drain-induced-barrier¬ 
lowering  (DIBL),  transconductance,  current  driving 
capability,  substrate  current,  and  degradation  of 
submicron  devices  were  characterized  for  various 
channel  dimensions,  doping  concentrations,  and 
gate  oxide  thicknesses.  Discussion  of  the  experi¬ 
mental  results  are  presented  in  later  sections. 


DEVICE  FABRICATION 
The  devices  used  in  this  study  were  fabricated 
using  an  NMOS  technology  for  which  the  tempera¬ 
ture  cycle  after  source/drain  implantation  (As, 
3xl015cm~2,  50  Kev,  O’  inclination)  was  limited  to 
900  °C  or  below  to  reduce  the  junction  depth. 
Conventional  g-line  optical  lithography  was  used  to 
define  all  levels  including  that  for  the  gate.  The 
ultra-fine  gates  were  obtained  by  calibrated  thin¬ 
ning  of  the  optically  defined  photoresist  features  in 
oxygen  plasma.  An  SEM  photograph  of  a  transis¬ 
tor  cross-section  with  Leff  =  0.22  pm  is  shown  in 
Fig.  1.  With  this  technique,  we  have  been  able  to 
fabricate  deep-submicron  MOSFET’s  routinely 
with  excellent  device  characteristics  usu»s  standard 
semiconductor  processing  equipment.  Transistors 
of  four  oxide  thicknesses  (  36  A,  56  A,  86  A,  and 
155  A  )  with  substrate  dopings  ranging  from 
6xl016cm-3  to  8xIOncm-3  were  fabricated.  The 
junction  depth  is  about  0.15  pm  for  all  devices.  In 
Fig.  2(a),  the  strong-inversion  drain  current  of  a 
device  with  L^t  =  0.15  pm  and  Tox  ®  36  A  is 
shown  which  exhibits  very  good  characteristics. 
This  device  also  exhibits  very  well-behaved 
subthreshold  characteristics  (Fig.  2(b))  with 
insignificant  punchthrough  problem.  The  practical 
channel  length  appears  to  be  limited  by  the  tech¬ 
nology  controllability  than  by  any  inherent  physical 
mechanism. 


DEVICE  CHARACTERISTICS 
AND  DISCUSSION 

(A)  Drain  current  and  short-channel  effects 
A  good  indicator  of  the  extent  of  short-channel 
effects  for  a  given  process  technology  is  the  thres¬ 
hold  voltage  shift  due  to  charge-sharing  and  DIBL 
effects.  In  Fig.  3(a),  the  linear-region  threshold 
voltage  is  plotted  versus  effective  channel  length 
for  four  oxide  thicknesses.  The  substrate  doping  for 
each  oxide  thickness  is  chosen  so  that  the  threshold 
voltages  for  long-channel  devices  are  in  the  range 
of  0.5  -  0.7  V.  The  threshold  voltage  shift  at  VDS 
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at  3V  versus  effective  channel  length  is  shown  in 
Fig.  3(b).  Both  figures  show  that  although  thicker 
oxide  devices  present  worse  threshold  shift  as 
expected,  the  threshold  voltages  do  not  suffer 
severe  change  until  the  effective  channel  length  is 
shorter  than  0.25  jam  for  all  oxide  thicknesses. 

The  transconductance  and  the  saturation  drain 
current  of  an  array  of  devices  are  plotted  in  Fig.  4. 
Very-well-behaved  trends  agree  with  a  simple  phy¬ 
sical  model  [5]-(6],  indicating  that  the  basic  phy¬ 
sics  is  well-understood.  The  model  equations  are 
summarized  below. 

JdSaT=  WeiT  vsai C ox  (  ^GS  "  Vih  ~  VDSAT  )  (0 


where 
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Es  is  the  effective  vertical  field  in  the  channel  and 
can  be  approximated  by  Es  =  [2QB  + 
Co*(vgs-vJ]/2£s.-  Qb  ^  the  depletion  bulk 
charge,  Eq  *  0.67  MV/cm,  n  =  1.6,  pg  *  670 
cm:/V  sec,  and  v^sgxlO6  cm  /sec.  The  same 
model  parameters  were  used  for  all  device  dimen¬ 
sions  and  oxide  thicknesses.  The  data  shown  in 
Fig.  4  have  been  corrected  for  the  source/drain  and 
contact  resistance  (=30Q/side).  The  inversion- 
layer  capacitance  which  is  more  important  for 
thin-oxide  devices  [7]  at  low  gate  bias,  has  not 
been  included  in  these  equations.  The  device  with 
36  A  gate  oxide  and  Leff  =  0.15  pm  exhibits  a 
maximum  transconductance  of  640  mS/mm,  the 
highest  ever  reported  for  a  MOSFET  at  room 
temperature.  While  simplistic  in  formulation,  these 
model  equations  do  provide  good  physical  as  well 
as  quantitative  understanding  of  the  current  perfor¬ 
mance  of  MOS  devices  down  to  submicron  regime. 
They  would  also  serve  well  as  first  order  calibra¬ 
tors  for  process  control  and  diagnosis. 


(B)  Hot-electron  effects 

The  substrate  current  peak  is  higher  for  a  dev¬ 
ice  with  shorter  channel  length  and  thinner  oxide 
as  expected.  The  peak  substrate  current  at  VDS  =  3 
V  versus  channel  length  for  different  oxide 
thicknesses  is  shown  in  Fig.  5.  More  informative 


plots  of  Isub^ds  versus  1/(vds-Vdsat>  [8]  for 
devices  with  Leff  =  0.3  pm  and  different  gate  oxide 
thicknesses  are  shown  in  Fig.  6.  The  straight  lines 
suggest  that  the  same  basic  physical  mechanism 
prevail  independent  of  channel  length.  However, 
the  slopes  for  the  lines  seem  to  deviate  from  the 
TgX1/3  dependence  predicted  by  the  model. 

In  the  degradation  studies,  the  device  lifetimes 
show  a  dependence  on  the  channel  length.  Results 
for  devices  with  86  A  gate  oxide  are  presented  in 
Fig.  7.  Similar  results  are  observed  for  other  oxide 
thicknesses.  This  channel-length  dependent  degra¬ 
dation  can  partly  be  explained  by  the  increase  of 
the  ratio  of  damaged  interface  area  to  the  total 
channel  area  as  the  channel  length  decreases.  A 
useful  variation  of  Fig.  7  which  provides  direct 
device  design  guidelines  is  shown  in  Fig.  8,  where 
the  maximum  supply  voltage  to  ensure  a  10-year 
device  lifetime  (corresponding  to  a  3%  change  in 
drain  current  [9])  is  plotted  against  channel  length. 
This  maximum  supply  voltage  appears  to  decrease 
as  the  channel  length  decreases.  For  a  quarter- 
micron  device,  the  maximum  supply  voltage  is  2.5 
V  suggesting  that  some  kinds  of  hot-electron- 
resistant  structures  are  still  desirable  even  if  the 
power  supply  is  lowered  to  3.3  V. 

CONCLUSION 

Performance  and  hot-electron  reliability  of 
deep-submicron  MOSFET’s  are  characterized  for 
various  oxide  thicknesses  and  substrate  dopings.  It 
is  found  that  the  basic  physics  involved  in  deep- 
submicron  devices  is  similar  to  that  of  micron-size 
devices.  Existing  physical  models  are  applicable  to 
these  devices  with  the  need  for  only  minor 
modifications.  The  device  degradation  results  sug¬ 
gest  that  the  maximum  supply  voltage  for  a  con- 
ventioal  quarter-micron  device  with  86  A  oxide 
thickness  should  be  limited  to  2.5  V  based  on  a 
10-year  lifetime  extrapolation. 
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Fig.  2  Drain  current  characteristic',  t  i  desp- 
submicron  MOSFET  fabricated  *  tn  the 
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Fig.  1  SEM  photograph  of  a  transistor  cross-section 
with  Ljfj  =  0.22  pm.  The  junction  depth  is  about 
0.15  pm. 
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Fig.  3  (a)  linear-region  threshold  voltage,  (b)  thres¬ 
hold  voltage  shift  at  VDS  =  3  V,  as  a  function  of 
effective  channel  length  for  various  oxide 
thicknesses. 
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Fig.  4  Saturation  characteristics  as  a  function  of 
effective  channel  length.  Symbols  represents  meas¬ 
ured  data.  Lines  are  calculated  from  a  simple  phy¬ 
sical  model,  (a)  transconductance,  (b)  drain  satura¬ 
tion  current. 
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Fig.  7  Device  lifetime  versus  subsiraif  :-rrent  per 
unit  channel  width  for  T0I  =  86  A  ard  various 
channel  lengths.  Lifetime  is  defined  n  '•  "t  .  j-rent 
degradation. 
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Fig.  5  Peak  substrate  current  at  VDS  =  3  V  as  a 
function  of  effective  channel  length  for  various 
oxide  thicknesses. 


Fig.  6  Normalized  substrate  currert  versus 
1/(VD$  —  VDSAT)  for  transistors  with  L^  *  0.3  Jim 
and  various  oxide  thicknesses. 
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Fig.  8  Maximum  operating  voltage  versus  effective 
channel  length  for  86  A  devices  based  on  a  10- 
year  lifetime  extrapolation. 
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A  bevel  has  been  etched  in  a  GaAs  epitaxial  film  grown  on  a  Si  substrate,  so  that  the  Raman 
spectrum  of  the  GaAs  layers  can  be  measured  as  a  function  of  distance  from  the  GaAs/Si 
interface.  The  amount  of  strain  and  disorder  in  the  GaAs  film  has  been  estimated  from  the 
GaAs  longitudinal  optical  phonon  line  shape  and  frequency.  Both  the  strain  and  the  amount  of 
disorder  were  found  to  decrease  with  increase  in  the  distance  from  the  interface. 


Recently  there  has  been  much  interest  in  GaAs  epitaxial 
films  grown  on  Si  substrates.  This  is  due  to  the  many  poten¬ 
tial  applications  which  can  result  from  combining  the  fast 
speed  and  excellent  electro-optic  properties  of  GaAs  with 
the  highly  developed  Si  technology. 1-14  Two  major  problems 
have  been  encountered  in  the  growth  of  high  quality  GaAs 
films  on  Si  substrates.  Firstly  there  is  a  4. 1  %  lattice  mis¬ 
match  between  the  two  materials,  and  secondly  the  random 
occupation  of  the  Ga  and  As  sublattices  leads  to  the  presence 
of  antiphase  domains.  The  first  effect  results  in  high  densities 
of  dislocations  and  possibly  also  strain  in  the  GaAs  films. 
The  antiphase  domains  contain  large  amounts  of  antisite  de¬ 
fects  at  the  domain  boundaries.  Raman  scattering  has  been 
useful  in  the  study  of  defects  in  GaAs  caused  by  ion  implan¬ 
tation,  polishing,  and  high-energy  neutron  and  electron  irra¬ 
diations. 5-10  We  report  here  the  first  Raman  study  of  epitax¬ 
ial  films  of  GaAs  grown  on  Si.  From  this  study  we 
determined  the  crystallinity  and  strain  in  the  GaAs  film  as  a 
function  of  its  distance  from  the  interface. 

Our  sample  was  a  2.3-/atn-thick  GaAs  film  grown  by 
molecular  beam  epitaxy  (MBE)  on  a  Si  substrate  oriented  4* 
off  the  [  100]  direction  towards  the  [  1 10]  axis.  To  avoid  the 
formation  of  islands,  the  first  few  hundred  angstroms  of 
GaAs  was  grown  at  T<  350  *C.  This  thin  layer  of  GaAs  was 
amorphous.  After  the  Si  surface  was  completely  covered  by 
GaAs,  the  substrate  temperature  was  increased  to  600  *C. 
Two  layers  of  strained  InGaAs-GaAs  superlattices,  each  of 
five  periods,  are  included  in  the  GaAs  film  as  shown  in  the 
inset  of  Fig.  1.  The  strained  superlattices  bend  the  disloca¬ 
tions  and  drastically  reduce  the  density  of  the  dislocations 
which  can  propagate  through  the  superlattice.  A  small  angle 
bevel  was  etched  by  gradually  dipping  the  sample  into 
NH«0H:H20::H20  (1:1 :20)  to  expose  the  GaAs  layers  near 
the  interface  as  shown  in  the  inset  of  Fig.  1.  To  maintain  a 
uniform  etching  rate,  the  sample  passes  over  de-ionized  wa¬ 
ter  before  touching  the  etchant. 

The  Raman  experiment  was  performed  at  room  tem¬ 
perature  in  thex(yr)3E  backscattering  geometry,  where x.^.z 
refer  to  the  [  100],  [010],  and  [001  ]  crystallographic  3xes, 
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respectively.  A  SPEX  double  spectrometer  and  a  multichan¬ 
nel  detector,  which  are  similar  to  those  described  in  Ref.  1 1 , 
were  used  to  analyze  the  Raman  signal.  Typically  300  m  W  of 
Ar+  laser  power  (5145  A)  was  focused  with  a  cylindrical 
lens  to  a  spot  of  about  10  X  40  (/am  )J  on  the  sample  with  the 
longer  side  parallel  to  the  edge  of  the  bevel.  The  GaAs  lattice 
temperature  at  the  focus  was  estimated  to  be  360  K. 

The  Raman  spectra  were  obtained  successively  by  trans¬ 
lating  the  sample  in  a  direction  perpendicular  to  the  bevel 
edge  with  a  micrometer.  From  the  known  bevel  angle  and 
distance  of  laser  focus  to  the  bevel  edge,  the  thickness  of  the 
GaAs  film  probed  by  the  laser  can  be  calculated.  The  accura¬ 
cy  of  this  thickness  measurement  is  about  300  A.  When  the 
focus  was  scanned  towards  the  bevel  edge  from  the  GaAs 
side,  the  GaAs  Raman  signal  decreased  while  the  Si  optical 
phonon  (at  519.5  cm-1)  intensity  increased  (see  Fig.  1). 
Since  the  laser  beam  has  to  penetrate  the  GaAs  film  to  reach 
the  Si  substrate,  we  can  determine  the  absorption  coeflficient 
of  the  GaAs  film  from  the  variation  of  the  Si  phonon  intensi¬ 
ty  provided  that  this  absorption  coefficient  is  not  varying 


FIG.  I.  Variation  of  the  Si  optical  phonon  and  GaAs  LO  phonon  intensities 
with  the  thickness  of  the  GaAs  film.  The  solid  curve  is  a  fit  to  the  Si  phonon 
intensity  calculated  by  assuming  that  the  GaAs  film  has  an  absorption 
length  ofl  580  A.  The  broken  curve  guides  the  eyes.  The  GaAs  LO  phonon 
intensities  have  been  multiplied  by  a  factor  of  12  in  order  to  be  plotted  on  the 
same  scale  as  the  Si  phonons.  The  inset  shows  the  sample  geometry. 
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with  the  GaAs  film  thickness.  The  solid  curve  in  Fig.  1 
shows  a  fit  to  the  experimental  points  with  a  constant  GaAs 
penetration  depth  of  1580  ±  20  A.  This  value  is  consistent 
with  that  of  bulk  GaAs.12 

The  GaAs  longitudinal  optical  (LO)  phonon  intensity 
in  Fig.  I  was  actually  the  sum  of  the  Raman  scattering  from 
GaAs  layers  at  different  depths  from  the  interface  up  to  the 
surface.  To  determine  the  Raman  spectra  of  the  GaAs  layers 
as  a  function  of  its  distance  to  the  interface,  we  have  adopted 
the  following  procedure.  First  the  Raman  spectrum  of  GaAs 
was  measured  with  the  laser  focus  at  a  spot  as  close  as  possi¬ 
ble  to  the  edge  of  the  bevel.  This  is  taken  to  be  the  Raman 
spectrum  of  the  GaAs  layer  at  the  interface.  The  laser  focus 
is  then  translated  by  a  few  hundreds  of  angstroms  away  from 
the  interface.  The  resultant  Raman  spectrum  is  then  the  sum 
of  the  Raman  spectrum  from  the  GaAs  nearest  to  the  inter¬ 
face  plus  that  from  an  additional  overlayer.  By  subtracting 
one  Raman  spectrum  from  the  other,  we  obtain  the  Raman 
spectrum  of  any  additional  layer.  The  attenuation  of  light 
due  to  absorption  in  the  GaAs  film  has  to  be  corrected  for  in 
this  subtraction  process.  This  process  can  be  continued  until 
the  spot  reaches  the  surface  of  the  GaAs  film.  Figure  2  shows 
the  normalized  GaAs  LO  phonon  line  shape  of  GaAs  layers 
at  different  distances  to  the  interface  obtained  in  this  way. 
When  the  layers  are  far  away  from  the  interface,  their  LO 
phonon  spectra  are  narrow,  symmetric,  and  Lorentzian  in 
shape.  After  correction  for  the  instrumental  response,  the 
intrinsic  GaAs  LO  phonon  linewidth  is  about  2  cm  1  for 
layers  farther  than  half-micro.i  from  the  interface.  The 
GaAs  LO  phonon  peaks  become  broader  and  more  asymme¬ 
tric  when  the  layers  are  closer  to  the  interface.  Figure  3(a) 
shows  the  asymmetry  and  broadening  of  the  GaAs  LO 
phonon  as  a  function  of  the  layer  distance  to  the  interface. 

In  a  perfect  crystal  only  the  zone  center  (q  =  0)  phonon 
modes  are  observed  in  one-phonon  Raman  scattering  due  to 
quasi-momentum  conservation.  The  phonon  line  shapes  are 


RAMAN  FREQUENCY  (cm'1) 


FIG.  2.  Raman  spectra  of  GaAs  layers  ai  different  distances  to  the  interface 
obtained  by  taking  the  difference  of  two  Raman  spectra  measured  from  two 
adjacent  spotJ  along  the  bevel.  At  0-500  A:  B:  500-800  A.  C  800-1000  A, 
D:  1000-1500  A;  E:  1500-2100  A;  F:  2100-5200  A. 
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FIG.  3.  (a)  Linewidth,  frequency  shift,  and  asymmetry  of  the  GaAs  LO 
phonon  as  a  function  of  distance  to  the  interface.  These  results  have  been 
obtained  from  Raman  spectra  of  the  GaAs  LO  phonon  measured  with  bet¬ 
ter  resolution  than  those  shown  in  Fig.  2.  The  inset  shows  quantitatively  the 
definition  of  the  parameters  1~,  and  f ,  and  the  asymmetry  parameter 
r,/ry  (b)  Strain  and  spatial  correlation  lengths  of  the  GaAs  film  as  a 
function  of  distance  to  the  interface  deduced  from  (a). 


usually  Lorentzian  with  the  phonon  linewidths  determined 
by  the  phonon  lifetime.  In  a  disordered  crystal,  however 
phonons  can  be  localized  to  regions  of  extent  L.  Asa  result 
quasi-momentum  conservation  is  relaxed,7'9  and  the  LO 
phonon  Raman  peak  becomes  asymmetric  and  shifted 
towards  lower  frequency.  On  the  other  hand,  compressive 
stress  can  shift  the  LO  phonon  towards  higher  frequen- 
Cy  S.U.U.IS  Thus  tj,e  LO  phonon  line  shape  and  peak  posi¬ 
tions  contain  information  on  both  the  crystal  quality  and  on 
the  amount  of  strain  in  the  GaAs  film.  We  adopt  this  spatial 
correlation  model  to  interpret  our  experimental  results.7  9  In 
this  model,  L  can  be  deduced  from  one  of  these  quantities: 
shift  in  LO  phonon  frequency,  LO  phonon  width,  or  the 
asymmetry  of  LO  phonon  line  shape.  Figure  3(b)  shows  the 
correlation  lengths  of  the  GaAs  film  as  a  function  of  distance 
to  the  interface  deduced  from  the  linewidth  and  the  asymme¬ 
try.  Both  show  the  same  increase  in  the  correlation  length  as 
the  GaAs  film  becomes  farther  from  the  GaAs/Si  interface. 

We  found  that  the  GaAs  LO  phonon  frequency  is  blue 
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shifted  for  GoAs  layers  nearer  to  the  interface  in  contrast  to 
the  prediction  of  the  spatial  correlation  theory.  We  attribute 
this  difference  to  the  presence  of  compressive  strain  in  the 
GaAs  layers  near  the  interface  due  to  the  lattice  mismatch 
between  GaAs  and  Si.  To  determine  the  magnitude  of  this 
strain  we  assume  a  uniform  and  isotropic  two-dimensional 
strain  for  simplicity.  We  also  assume  that  the  GaAs  film  free 
surface  is  unstrained  and  use  us  LO  phonon  frequency  as 
reference  co0.  From  the  phenomenological  theory  in  Refs. 
8,13-15  the  strain-induced  frequency  shift  A<u  is  given  by 
Aa)/a>0  =  —  1.3 A",  where  X  is  the  strain.  Figure  3(b)  shows 
the  variation  in  the  compressive  strain  as  a  function  of  dis¬ 
tance  from  the  interface.  The  strain  measured  in  our  sample 
is  much  smaller  than  the  4%  strain  expected  from  the  lattice 
mismatch.  This  suggests  that  the  strain  in  GaAs  films  on  Si  is 
relieved  by  dislocations.  For  GaAs  layers  more  than  ~  2000 
A  away  from  the  interface,  the  strain  is  below  our  measure¬ 
ment  sensitivity. 

A  very  small  TO  phonon  Raman  peak  was  also  observed 
in  our  sample.  Since  the  TO  phonon  is  forbidden  by  a  selec¬ 
tion  rule  in  our  scattering  geometry,  we  have  interpreted  this 
peak  as  due  to  disorder-induced  scattering.  The  intensity  of 
the  TO  phonon  relative  to  the  LO  phonon  increases  towards 
the  interface.  This  is  consistent  with  our  result  that  disorder 
increases  for  GaAs  layers  closer  to  the  interface. 

In  conclusion,  we  have  used  the  Raman  scattering  from 
a  bevel  etched  on  the  GaAs  film  grown  on  a  Si  substrate  by 
MBE  to  examine  the  GaAs  quality  as  a  function  of  its  dis¬ 
tance  to  the  interface.  The  lattice  spatial  correlation  and 
strain  were  deduced  from  the  LO  phonon  line  shapes  and 


frequency  shifts.  We  found  that  both  the  strain  and  the  disor¬ 
der  in  the  GaAs  film  decre  ised  in  the  first  2000  A  to  a  level 
below  the  sensitivity  of  the  Raman  measurement. 
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RAMAN  SCATTERING  SPECTRUM  ALONG  A  BEVEL 
ETCHED  GaAs  ON  Si,  TEM  STUDY  AND 
GaAs  P-I-N  PHOTODETECTOR  ON  Si 
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Abstract 

Raman  scattering  is  measured  along  a  bevel  etched  GaAs  epitaxial  film  grown  on  Si  by 
molecular  beam  epitaxial  (MBE).  From  the  correlation  length  profile  of  Raman  scattering, 
most  dislocation  lines  are  confined  in  the  200oA  regions  close  to  the  interface.  The  strain 
profile  calculated  from  the  Raman  peak  shift  shows  that  about  0.6%  compressive  strain  exists 
near  the  interface  because  of  lattice  mismatch.  However,  as  one  moves  away  Grom  the  inter¬ 
face,  the  compressive  strain  is  gradually  counterbalanced  by  thermal  expansion.  Transmission 
electron  microscope  (TEM)  studies  of  the  local  dislocation  image  and  properties  show  that  an 
ultra  clean  Si  surface  is  essential  for  dislocation  confinement.  From  high  resolution  TEM,  we 
find^that  the  distance  between  dislocations  at  the  interface  is  nonuniform,  varying  from  5oX  to 
125A  with  an  average  distance  at  81A.  Finally,  a  GaAs  p-i-n  photodetector  on  Si  substrate  is 
fabricated.  Even  though  a  normal  photoresponse  curve  is  obtained,  the  high  dark  cutrent 
(50nA)  and  relatively  low  responsivity  (0.01A/W)  show  that  the  material  quality  needs  to  be 
further  improved  to  make  a  minority  carrier  vertical  transition  device. 

Introduction 

Lately  the  research  on  epitaxial  growth  of  GaAs  on  Si  substrate  is  developing  quickly 
because  it  not  only  invokes  a  number  of  interesting  scientific  challenges  but  also  provides 
practical  benefits  in  applications  such  as  GaAs  VLSI  digital  circuit  and  opto-electronic  integra¬ 
tion.1'5  However,  several  fundamental  problems  exist  in  the  epitaxial  growth  of  GaAs  on  Si. 
First,  the  growth  of  a  polar  crystal  on  a  non-polar  crystal  tends  to  form  anti-phase  domains 
(APD)  which  induce  static  charge  sheets  in  the  GaAs  thin  film.  Second,  the  relatively  large 
lattice  mismatch  (=4%)  between  GaAs  and  Si  will  induce  many  dislocations  near  the  interface 
and  some  of  them  may  propagate  towards  the  crystal  surface.  Third,  the  large  difference  in 
thermal  expansion  coefficients  between  GaAs  and  Si  causes  a  very  high  stress  which  often 
makes  the  sample  crack  during  wet  etch  processing. 

To  investigate  the  strain  and  dislocabon  properties  as  mentioned  above,  we  perform  a 
Raman  scattering  experiment  which  can  measure  the  strain  and  dislocation  density  profile 
from  the  GaAs  and  Si  interface  to  the  GaAs  surface.  On  the  other  hand,  cross  section  TEM 
and  high  resolution  TEM  provide  us  with  information  about  dislocation  image  and  properties. 
A  GaAs  on  Si  p-i-n  photodetector  is  also  fabricated  for  the  first  time  and  the  results  are 
compared  with  a  normal  GaAs  p-i-n  photodetector. 
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Raman  Scattering 

GaAs  is  grown  on  a  Si  substrate  oriented  4  degrees  off  the  [100]  direction  towards  the 
[110]  axis  by  MBE.  As  shown  in  the  insert  of  Fig-1,  a  bevel  of  0.1  degree  has  been  etched 
in  the  GaAs  epitaxial  film,  so  that  the  Raman  spectrum  of  GaAs  layers  can  be  measured  as  a 

Q 

function  of  distance  from  the  GaAs/Si  interface.  The  514JA  line  from  an  A t*  laser  is  used  to 
trigger  the  Raman  signal.  The  Raman  experiment  is  performed  at  room  temperature  in  the 
x(yz)x  back-scattering  geometry,  where  x,y,z  refer  to  the  [100],  [010]  and  [001]  crystallo¬ 
graphic  axes  respectively.  In  this  scattering  geometry,  the  longitudinal  optical  (LO)  phonon  of 
GaAs  is  allowed  by  the  selection  rules,  while  the  transverse  optical  (TO)  phonon  is  forbidden. 
The  LO  phonon  Raman  scattering  intensity  variation  with  distance  from  the  interface  is  shown 
in  Fig-1.  We  can  determine  the  transmission  coefficient  of  the  GaAs  film  form  the  variation  of 
the  Si  phonon  intensity  provided  that  this  absorption  coefficient  is  not  varying  with  the  GaAs 
film  thickness.  The  solid  curve  in  Fig-1  shows  a  fit  to  the  experimental  points  with  a  constant 
GaAs  penetration  depth  of  1580±20A. 

The  normalized  Raman  spectrum  at  different  depths  from  the  interface  is  shown  in  Fig-2. 
The  GaAs  LO  phonon  peaks  become  broader  and  mote  asymmetric  when  the  layers  are  closer 
to  the  interface.  This  phenomenon  can  be  explained  as  follows. 

In  a  perfect  crystal,  only  zone  center  (q«0)  phonon  modes  are  excited  for  one-phonon 
Raman  scattering  because  of  quasi-momentum  conservation.  For  a  disordered  lattice,  how¬ 
ever,  phonons  can  be  localized  to  a  region  of  space  L,  and  as  a  result  quasi-momentum  con¬ 
servation  will  be  relaxed  so  that  off-center  phonons  (q*Q)  can  be  excited.6-®  The  participation 
of  off-center  phonon  scattering  causes  the  asymmetric  non-Loren tzian  lineshape  near  the  inter¬ 
face.  On  the  other  hand,  the  phonon  linewidth  near  the  interface  is  broader  than  the  phonon 
linewidth  for  layers  farther  than  half  micron  from  the  interface  because  the  phonon  lifetime  is 
reduced  by  disordered  material.  Both  the  LO  phonon  lineshape  and  linewidth.  therefore,  con¬ 
tain  infonnation  about  the  crystal  quality.  We  adopt  the  spatial  correlation  model  developed 
by  Richter6  and  Tiong  et  al.,*  to  calculate  the  spatial  correlation  length  L  from  either  the  pho¬ 
non  width  or  the  lineshape  asymmetry.  Fig-3  shows  the  correlation  lengths  of  the  GaAs  film 
as  a  function  of  distance  from  the  interface.  Those  calculated  correlation  lengths  are  related 
to  the  average  distance  between  two  dislocation  lines.  The  abrupt  increase  of  correlation 
length  around  2000A  from  the  interface  shows  that  most  dislocations  can  not  propagate 
further. 

In  addition,  the  strain  profile  can  also  be  calculated  from  the  Raman  peak  frequency  shift 
if  a  uniform  and  isotropic  two-dimensional  strain  is  assumed.7  As  shown  in  Fig-3,  a  compres¬ 
sive  strain  is  present  in  the  GaAs  layers  near  the  interface  is  present.  This  may  be  due  to  the 
lattice  mismatch  between  GaAs  and  Si.  This  compressive  strain  from  lattice  mismatch  is 
counterbalanced  by  the  tensile  strain  from  thermal  expansion,  so  the  GaAs  surface  becomes 
almost  strain  free  for  a  2.3um  film. 

© 

For  the  region  thicker  than  5000A,  the  dislocation  density  is  too  low  to  be  detected  by 
Raman  scattering  techniques.  The  discrete  feature  and  local  image  of  dislocation  lines  will  be 
studied  by  TEM. 


Figs-4  and  5  show  the  TEM  cross  sections  of  two  GaAs  on  Si  samples  grown  by  MBE. 
Fig-4  shows  the  growth  of  GaAs  on  an  oxygen  and  carbon  contaminated  Si  surface.  Very 
characterisdc  triangular  protrusions  are  observed  near  the  Si  surface  giving  rise  to  a  high 
dislocation  density.  Many  dislocations  and  stacking  faults  generated  at  the  interface  propagate 
to  the  surface  layer  of  GaAs.  The  density  of  dislocations  estimated  in  the  near  surface  area  is 
in  the  range  of  5  to  9xlO,cm~2.  If  the  Si  surface  is  clean  (Fig-S),  triangular  features  ate  not 
formed  at  the  interface  and  the  density  of  stacking  faults  formed  at  the  interface  is  lower  as 
well.  Many  of  the  misfit  dislocations  are  tangled  and  confined  in  the  near  Si/GaAs  area 
within  0.3pm  from  the  Si  substrate.  This  him  quality  is  quite  acceptable  even  without  any 
strain  layer  superlatdce  or  thermal  annealing.  This  observation  indicates  that  an  ultra  clean  Si 
substrate  is  necessary  for  high  quality  GaAs  on  Si  growth. 

As  was  found  before,3  two  types  of  dislocations  are  formed  at  the  interface.  Most  of  the 
dislocations  (=8 5%)  have  their  Burgers  vector  parallel  to  the  interface  (Fig-6)  and  only  about 
15%  have  their  Burgers  vector  inclined  to  the  interface.  This  small  fraction  of  dislocations 
with  inclined  Burgers  vector  is  probably  the  reason  for  the  much  lower  density  of  dislocations 
=5x  107cm'2)  in  the  near  surface  area  of  the  GaAs.  The  average  distance  between  the  dislo¬ 
cations  at  the  interface  is  one  dislocation  per  81 X.  However,  the  distribution  of  the  disloca¬ 
tions  at  the  interface  in  not  uniform.  There  are  areas  where  the  distance  between  the  disloca¬ 
tions  (estimated  from  high  resolution  TEM  pictures)  is  only  50A  apart  but  there  are  other 
areas  where  dislocations  are  separated  as  much  as  123A.  Higher  stress  (=2.5%)  is  built  up  in 
the  areas  where  dislocations  are  closer  to  each  other.  Stress  is  almost  completely  released 
(=0.1%)  in  the  areas  where  dislocations  are  far  from  each  other.  One  would  expect  that  dislo¬ 
cations  which  are  close  together  interact  differently  from  dislocations  which  are  further  apart 
This  nonuniformity  will  influence  the  final  dislocation  distribution  in  the  top  surface  layer  of 
GaAs,  where  one  can  find  both  areas  which  are  dislocation  free  and  areas  which  have  higher 
dislocation  densities. 

P-I-N  Photodetector 

In  terms  of  device  characteristics,  most  majority  carrier  devices  like  FETs  have  been 
reported  by  several  groups1,2  showing  comparable  performance  with  those  devices  made  on  a 
GaAs  substrate.  LEDs  and  lasers  have  also  been  recently  reported  even  though  reliability  and 
lifetimes  are  much  worse  than  those  on  a  GaAs  substrate.9,10  Phoiodetectors.  however,  are 
almost  not  reported  yet  besides  the  lateral  conductive  photoconductors  which  is  most  insensi¬ 
tive  to  dislocations  and  other  defects.  We  fabricated  a  p-i-n  homojunction  GaAs  photodiode 
having  the  device  structure  shown  in  Fig-7.  The  commonly  used  two-step  growth  were 
adopted  and  three  InGaAs/GaAs  strain  layer  superlattices,  ten  50&I0oX  periods  each  super¬ 
lattice.  are  grown  as  dislocation  filtering  layers.  The  high  dark  current  (50nA  at  -0.2volt)  res¬ 
tricts  this  device  from  working  at  higher  bias.  The  normalized  photoresponse  with  respect  to 
wavelength  is  shown  in  Fig-8.  The  peak  responsivity  at  room  temperature  occurs  at  0.78|rm 
(1.59ev)  and  cuts  off  at  0.88pm  (1.4 lev).  The  lack  of  abnormal  peaks  at  longer  wavelength 
(up  to  1.5pm)  means  that  there  is  no  significant  deep  level  absorption.  High  dark  current  can 
be  explained  by  higher  dislocation  density  (approximately  103  to  104  times)  in  the  device 
active  region  than  in  the  GaAs  substrate.  The  responsivity  without  an  anti-reflection  coating 
is  about  0.01A/W,  which  is  also  relatively  low. 
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Conclusion 

In  summary,  we  use  Raman  scattering  to  scan  across  a  beveled  edge  of  GaAs  on  Si  to 
iind  out  the  strain  as  well  as  the  average  dislocation  density  profiles.  Compressive  strain  is 
present  at  the  interface  because  of  the  GaAs  and  Si  lattice  mismatch  but  this  strain  is  gradu¬ 
ally  counterbalanced  by  the  tension  from  thermal  expansion.  The  dislocadon  image  and  its 
propagation  are  studied  by  TEM,  which  indicates  that  an  extremely  clean  Si  surface  is  essen¬ 
tial.  The  distance  between  two  dislocations  near  the  interface  is  nonuniform.  Finally,  the  first 
GaAs  on  Si  p-i-n  photodetector  is  demonstrated.  Despite  its  high  dark  current  and  relatively 
low  tesponsivity,  a  normal  photoresponse  curve  without  any  deep  level  peak  is  observed. 
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Figure  1.  The  LO  phonon 
Raman  scattering  intensity 
profile.  The  epitaxial  structure 
of  the  sample  with  bevel  angle 
0.1  degree  is  also  shown. 
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Figure  2.  Normalized  Raman  scatter¬ 
ing  a:  different  distances  from  the 
interface.  Both  LO  phonon  and  TO 
phonon  peaks  can  be  identified. 


***■»■ 
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Figure  3.  Calculated  strain  and  corre¬ 
lation  length  profiles  from  Raman 
scattering.  The  correlation  length 
obtained  from  either  Raman  linewidth 
or  lineshape  asymmetry  corresponds 
to  the  average  dislocation  distance. 
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Figure  4.  TEM  picture  of  GaAs 
grown  on  an  oxygen  and  carbon  con¬ 
taminated  Si  surface.  A  number  of 
dislocauons  and  stacking  faults  are 
observed  and  many  of  them  propagate 
to  the  GaAs  surface  laver. 


Figure  5.  A  TEM  picture  of  GaAs 
grown  on  a  clean  Si  surface.  In  con¬ 
trast  to  Fig-4,  most  of  the  dislocations 
are  confined  near  the  interface  or  are 
bent  in  the  middle  so  that  the  GaAs 
surface  has  very  low  distocadon  den- 
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Figure  6.  High  resolution  TEM  picture  of  GaAs  and  Si  interface.  In  this  picture. 
We  look  at  (110)  plane  where  dislocations  at  the  interface  are  found  on  the  (100) 
plane  with  Burgers  vector  parallel  to  the  interface.  All  of  these  dislocations  in 
this  picture  are  not  inclined. 
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Figure  7.  Schematic  diagram  of  the 
GaAs  on  Si  p-i-n  photodetector. 


Figure  8.  The  normalized  pho¬ 
toresponse  curve  of  the  GaAs  on  Si 
p-i-n  diode. 
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GaAs  films  have  been  grown  on  Si  substrates  patterned  with  SiN  by  molecular  beam  epitaxy. 

The  pattern  consists  of  bare  Si  stripe  of  width  ranging  from  10  to  IOQ^m  surrounded  by  SiN 
on  both  sides  and  a  reference  area  of  bare  Si.  77  K  photoluminescence  (PL)  spectrum  and 
intensity  are  measured  on  the  single  crystalline  GaAs  films  grown  on  these  stripes  and  the 
reference  area.  For  1.5-  and  3-pm-thick  films,  PL  intensity  from  the  10/jm  stripe  shows  140% 
and  75%  increase  over  the  reference  area,  respectively.  This  remarkable  increase  in  the  PL 
intensity  is  believed  due  to  the  reduction  of  dislocations  inside  the  window  area.  The 
improvement  in  the  optical  quality  makes  selective-area  molecular  beam  epitaxy  a  very 
attractive  technique  for  the  fabrication  of  optical  devices  on  Si  substrate. 


Heteroepitaxy  of  GaAs  thin  films  on  Si  substrate  offers  a 
unique  opportunity  for  monolithic  integration  of  GaAs 
based  optoelectronics  and  high-speed  devices  with  Si  inte¬ 
grated  circuits  (IC’s).  Since  1984,  GaAs  on  Si  metal-semi¬ 
conductor  field-effect  transistors  IC's, '  high  electron  mobil¬ 
ity  transistors,  heterojunction  bipolar  transistors, 
light-emitting  diodes  and  lasers,  have  been  successfully  dem¬ 
onstrated  by  various  groups2  using  molecular  beam  epitaxy 
(MBE)  and  metalorganic  chemical  deposition  (MOCVD). 
Despite  these  impressive  achievements,  epitaxial  GaAs  films 
on  Si  substrates  have  one  inherent  problem  independent  of 
the  epitaxy  techniques  used;  namely,  the  high  threading  dis¬ 
location  density  ( 107-108/cm2)  arising  from  the  4%  lattice 
mismatch  between  the  GaAs  and  Si  lattice  constants.  High 
dislocation  density  and  large  tensile  stress  (low  10’  dynes/ 
cm2)  are  two  main  obstacles  to  be  overcome  before  GaAs  on 
Si  technology  can  realize  its  full  potentials. 

To  reduce  dislocation  density  away  from  the  heteroin¬ 
terface,  strained-layer  superlattices  (SLS's)  such  as 
InGaAs/GaAsJ  and  GaP/GaAs/GaP  As*  have  been  incor¬ 
porated  by  many  groups  during  the  epitaxial  growth  to  bend 
upward  propagating  dislocations.  In  this  letter,  we  study  the 
possibility  of  laterally  terminating  dislocations  through  the 
use  of  selective-area  molecular  beam  epitaxy  (SAMBE)5 
with  SiN  as  the  on-wafer  insulating  mask  In  SAMBE 
growth,  single  crystalline  epitaxial  films  are  formed  only  on 
the  exposed  area  of  crystalline  substrate  \reas  cohered  by 
noncrystalline  films  (such  asSiO.  SiN  ■  -e-.ult  n  the  growth 
of  polycrystalline  (poly)  films.  Since  r-  -ms  contain 
many  grain  boundaries,  dislocations  enema  “tern  from  the 
single  crystalline  areas  are  effectively  jrvr*w  'he  poly 
film.  We  report  in  this  letter  successful  >aMHF  <r"»th  of 
GaAs  films  on  Si  substrates  and  their  K  Pt  -  arictens- 
tics.  We  found  significant  increase  in  PL  - •>  ;  •  •  r  both 


■'  Also  Center  for  Advanced  Material.  Lawrence  9r-»r  . -»-i  i-> 


1.5-  and  3-^m-thick  films  grown  by  SAMBE  over  films 
grown  by  planar  MBE.  This  finding  provides  strong  evi¬ 
dence  that  crystalline/poly  interface  is  responsible  for  the 
reductions  of  dislocation  density. 

Heavily  doped  /t7' -Si  (001 )  substrate  tilted  3.5*  towards 
( 1 10)  axis  is  used  in  this  work.  The  wafer  is  first  cleaned 
repeatedly  in  5:1  H2S04.H202  before  depositing  800  k  of 
SiN  by  thermal  CVD.  The  wafer  is  patterned  and  etched  to 
form  2-mm-long  stripes  of  exposed  Si  with  widths  ranging 
from  10  to  100 /rm  along  (T 10)  orientations.  Three  identical 
stripes  separated  by  20  p m  of  SiN  are  grouped  together  as 
shown  in  Fig.  1 .  Next  to  the  stripe  pattern  is  a  4  mm  by  4  mm 
area  of  bare  Si.  GaAs  film  grown  over  this  area  is  used  as  a 
reference  for  comparison  since  its  properties  are  identical  to 
GaAs  film  grown  over  a  patternless  substrate.  The  etching  of 
SiN  is  done  by  reactive  ion  etching  (RIE)  using  a  mixture  of 
CHFj  and  O-  and  a  20%  of  overetch  time.  Etching  of  Si 
substrate  is  estimated  to  be  100  A.  After  RIE,  the  wafer  is 
thoroughly  cleaned  and  processed  by  the  usual  pre-MBE 
preparation  procedures6  before  loading  into  the  MBE  cham¬ 
bers.  A  two-step  MBE  growth  procedure  is  used  where  the 
initial  2000  A  layer  grown  at  a  growth  rate  of  0.6/itn/h  and 
substrate  temperature  of  400  *C  is  followed  by  a  normal 
growth  rate  of  1  pmAi  at  substrate  temperature  of  570  *C. 
Two  undoped  samples  having  a  thickness  of  1.5  and  3  pm  are 
grown  in  this  way,  respectively.  The  films  grown  inside  the 
stripe  windows  show  identical  color  and  texture  compared 
with  the  film  grown  over  the  patternless  area  when  viewed 
under  the  microscope.  Significant  lateral  growth  in  ( 1 10) 
direction  is  observed  as  in  the  case  of  SAMBE  on  GaAs 
substrate.7 

The  two  samples  are  mounted  next  to  each  other  in  a 
liquid-nitrogen  Dewar  sitting  on  a  movable  X-Y  stage.  A  1- 
pm- thick  unintentionally  doped  GaAs  film  is  also  included 
as  a  standard  sample  for  comparison.  The  488  nm  line  of  an 
Ar*  laser  with  a  power  of  150  mW  focused  to  1.5  mm  by 
0.76  mm  spot  size  is  used  as  the  excitation  source  in  the 
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FIG.  1.  Layout  of  the  stripe  windows  (not  in  full  detail).  The  darkened 
region  is  area  covered  with  SiN  film.  The  rectangular  bo*  formed  by  the 
broken  line  represents  the  slit  opening  where  only  luminescence  inside  the 
bo*  enters  the  spectrometer. 


photolumincscence  (PL)  measurement.  Relative  low  power 
density  (28  W/cm3)  is  used  here  to  minimize  the  tempera¬ 
ture  variations  of  GaAs  films  inside  different  window  sizes. 
The  detector  consists  of  a  SPEX  double  grating  spectrom¬ 
eter  and  a  scanning  photocounting  GaAs  detector  cooled  to 
—  20  *C.  The  spectrometer  is  scanned  from  12  350  cm-1 
(8097  A)  to  11  750  cm-1  (8510  A)  during  the  measure¬ 
ment.  The  scanning  rate  of  the  spectrometer  is  set  at  2  cm  ~ '/ 
s  with  a  spectral  resolution  of  4  cm'"1.  The  luminescence 
signal  from  the  sample  is  first  imaged  onto  an  adjustable  slit 
before  entering  the  spectrometer.  By  adjusting  the  width  of 
the  slit,  only  luminescence  signal  from  areas  under  study  is 
detected.  The  large  size  of  the  laser  spot  covers  illumination 
over  a  large  area  of  the  sample  and  makes  accurate  compari¬ 
son  meaningful.  Before  each  measurement  the  wafer  posi¬ 
tion  is  carefully  adjusted  to  maximize  detector  counts  at  a 
wavelength  near  the  PL  peak.  This  ensures  that  all  the  emis¬ 
sion  from  a  particular  area  of  the  wafer  under  study  is  col¬ 
lected. 

The  PL  spectra  of  a  standard  GaAs  film,  a  GaAs  on  Si 
film  inside  10,um  and  100 /4m  stripes,  and  in  the  reference 
area  are  shown  in  Fig.  2  for  the  1.5-  and  3-/im-thick  samples, 
respectively.  All  GaAs  on  Si  films  show  two  clearly  resolved 
PL  peaks  (A  and  B)  shifted  towards  longer  wavelength  com¬ 
pared  to  standard  GaAs  sample.  This  has  been  observed  by 
many  groups4,8'9  and  is  generally  attributed  to  the  splitting  of 
valence  bands  by  the  biaxial  tensile  stress  inside  the  film.  The 
calculated  spectra  taking  into  account  of  a  uniform  biaxial 
tensile  stress,  exciton  transitions,  and  temperatures  of  excit¬ 
ed  carriers  are  also  plotted  as  solid  lines  in  Fig.  2.'°  The 
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FIG.  2.  PL  spectra  of  standard  GaAs  sample  and  GaAs  on  Si  samples  inside 
a  4  mm  by  4  mm  reference  area.  The  calculated  PL  spectra  ( shown  in  solid 
line)  using  a  biaxial  stress  model'0  are  superimposed  on  the  measured  data. 
Excellent  agreement  between  the  measured  and  calculated  spectra  is 
achieved  here. 


agreement  between  the  calculated  and  measured  spectrum  is 
excellent.  From  the  theoretical  calculation,  the  biaxial  ten¬ 
sile  stress  not  only  splits  the  J  —  3/2  and  /=  1/2  valence 
bands  but  also  causes  an  overall  shift  of  the  two  splitted 
bands-  The  tensile  stress  calculated  from  the  shift  of  peak  A 
and  B  and  their  splitting  agrees  well.  The  tensile  stress  is 
found  to  be  3.2  ±  0.4  kbar  for  the  1.5-/im  film  and  3.7  ±  0.5 
kbar  for  the  3-^m  film  respectively  at  77  K.  The  PL  peak 
wavelengths  from  the  film  inside  the  10-/xm  width,  the  100- 
p m  width,  and  the  reference  area  show  no  appreciable  differ¬ 
ence,  indicating  that  the  tensile  stress  of  the  films  is  not  af¬ 
fected  by  the  width  of  the  window  stripes  in  which  the  films 
are  grown.  The  linewidth  of  the  peak  B  is,  however,  broader 
for  the  lO^um-wide  film  compared  to  the  100-/um-wide  film 
probably  due  to  some  nonuniform  distribution  of  stress  in¬ 
side  the  \0-fim  film.  The  total  PL  intensity  is  calculated  by 
first  subtracting  off  the  dark  counts  of  the  detector  and  then 
integrating  PL  counts  over  all  frequency.  PL  intensity  from 
various  stripes  is  then  normalized  with  respect  to  their  win¬ 
dow  area,  the  slit  size,  and  the  PL  intensity  of  the  standard 
GaAs  sample.  The  results  are  shown  in  Fig.  3  for  both  sam¬ 
ples.  PL  intensity  on  the  4  mm  by  4  mm  reference  area  is 
taken  as  the  PL  intensity  for  infinite  stnpe  width.  The  im¬ 
provement  in  PL  intensity  for  the  film  in  the  10  /im  stnpe 
over  the  100  fim  stripe  is  over  140%  for  the  1.5-/4m-thick 
sample  and  75%  for  the  3-/rm-thick  sample.  Indirect  evi¬ 
dences  from  high-temperature  annealing  experiments" 
indicate  that  increase  in  PL  intensity  is  related  to  the  reduc¬ 
tion  of  dislocation.  Although  the  dislocation  reduction 
mechanism  is  not  yet  clear  at  this  stage,  we  believe  that  in¬ 
crease  in  PL  intensity  resulted  from  dislocation  density  re¬ 
duction  inside  the  select.ve-area  grown  film.  It  is  expected 
that  careful  transmission  electron  microscopy  (TEM)  stud¬ 
ies  will  reveal  the  mechanism  of  dislocation  termination  on 
the  poly/crystalline  boundary. 
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FIG.  3.  PL  intensity  normalized  with  respect  to  the  stripe  area  and  PL  in¬ 
tensity  of  standard  GaAs  sample  for  Sim  of  different  stripe  widths.  The  4 
mm  by  4  mm  reference  area  is  taken  as  having  inSnite  stripe  width. 


In  summary,  we  have  found  significant  increase  in  the 
PL  intensity  for  GaAs  films  grown  on  Si  substrates  by 
SAMBE.  We  attribute  the  increase  in  PL  intensity  to  the 
reduction  of  dislocation  density.  The  PL  spectra  of  GaAs  on 
Si  films  can  be  fully  explained  by  the  presence  of  biaxial 
tensile  stress  in  the  film.  It  is  also  shown  that  the  presence  of 
poly-GaAs  on  either  side  of  the  strip  does  not  change  the 
stress  of  the  crystalline  GaAs  film  grown  in  between.  The 


remarkable  improvement  in  the  film  quality  makes  SAMBE 
a  very  attractive  technique  for  the  fabrication  of  optical 
source  and  detector  on  Si  substrates. 
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Photoluminescence  of  GaAs  films  grown  on  Si  substrate  has  been  investigated  quantitatively  at 
77  K.  The  peak  shift  and  splitting  of  the  exciton  luminescence  are  shown  to  result  from  tensile 
stress  in  the  film.  Information  on  carrier  lifetime  has  been  deduced  from  the  line  shape  of  the 
photoluminescence. 


GaAs  film  grown  on  Si  substrate  by  molecular  beam 
epitaxy  has  been  of  much  interest  recently. 1-1  Many  different 
types  of  devices  have  been  fabricated  successfully  on  these 
thin  films.  But  there  are  still  serious  problems  remaining  to 
be  solved.  One  such  problem  is  the  high  density  of  disloca¬ 
tions  resulting  from  the  lattice  mismatch  between  GaAs  and 
Si.  This  high  dislocation  density  results  in  short  minority- 
carrier  lifetimes  making  it  difficult  to  fabricate  high-quality 
laser  diodes.  The  second  problem  is  thermal  stress  which 
results  from  difference  in  thermal  expansion  coefficients 
between  GaAs  and  Si.  Several  techniques  have  been  used  to 
measure  this  thermal  stress  in  GaAs  films  grown  on  Si.  Mea¬ 
surement  of  the  stress-induced  curvature  is  possible  for  large 
wafers  such  as  those  grown  by  metalorganic  chemical  vapor 
deposition.5  The  lattice  constant  of  the  strained  GaAs  films 
can  also  be  measured  by  x-ray  diffraction.6  But  both  tech¬ 
niques  measure  the  magnitude  of  the  stress  averaged  over 
large  areas.  The  line  shape  and  frequency  shirt  of  phonon 
Raman  scattering  have  been  shown  to  yield  information  on 
both  the  crystal  quality  and  stress.’  This  technique  suffers 
from  low  sensitivity  .  Determination  of  stress  in  GaAs  film 
on  Si  using  photoluminescence  technique  has  been  reported 
previously.'1  Zemon  et  al.  found  that  the  band-edge  lumines¬ 
cence  peak  was  shifted  and  split.  They  suggested  that  this 
shitt  and  splitting  were  caused  by  strain  and  could  be  used  to 
determine  the  thermal  stress  in  the  GaAs  film.9  In  this  letter 
we  report  a  quantitative  analysis  of  the  photoluminescence 
(PL)  of  several  GaAs  on  Si  samples  with  different  film 
thicknesses  and  widths.  We  show  that  from  an  analysis  of  the 
splitting  of  the  PL  peaks,  their  line  shapes,  and  relative  inten¬ 
sities.  one  can  deduce  besides  the  thermal  stress  other  useful 
information  on  sample  quality  such  as  minority-camer  life¬ 
time. 

Our  GaAs  films  were  grown  by  molecular  beam  epitaxy 
on  either  Si  or  GaAs  substrates.  The  film  thtcknesses  are 
listed  in  Table  I.  Sample  D  is  only  10  pi  wide  while  other 
samples  were  more  than  2  mm  wide.  The  samples  were  unin¬ 
tentionally /i-type  doped  to  less  than  I0l6em  "  The  5  14-5-  A 
line  of  an  Ar  -  laser  was  used  to  excite  the  PL  Since  the  line 
shape  of  the  PL  spectra  changed  at  high  excitation  power 
densities,  the  laser  power  density  was  kept  at  a  low  level  of  29 
W/cm:  which  we  will  show  to  cause  minimal  heating  of  the 
samples.  A  Spex  1403  double  monochromator  and  i  photon 


counting  system  were  used  to  measure  the  PL  spectra.  The 
experiment  was  performed  at  liquid-nitrogen  temperature 
with  all  the  samples  mounted  side  by  side  on  the  same  cold 
finger.10 

Sampled  was  assumed  to  be  free  of  thermal  stress  since 
its  substrate  was  also  GaAs.  Indeed  its  PL  spectra  showed 
one  sharp  peak  (see  Fig.  1 )  with  no  detectable  sign  of  split¬ 
ting.  The  frequency  of  this  peak  was  used  as  the  reference  for 
determining  the  stress-induced  shift  in  the  PL  of  other 
GaAs/Si  samples.  In  contrast  the  PL  spectra  of  all  the 
GaAs/Si  films  showed  two  peaks  similar  to  those  reported 
by  Zemon  et  al Although  these  spectra  were  qualitatively 
similar  to  each  ether,  the  widths  of  the  PL  were  quite  differ¬ 
ent  between  spectra  B,  C,  and  D. 

To  interpret  these  PL  spectra  quantitatively,  we  have 
used  the  following  model.  First  we  assume  that  each  peak  in 
the  PL  spectrum  contains  two  contributions:  transitions  due 
to  the  lowest  ( Is)  exciton  level  and  the  corresponding  band- 
to-band  transition.  These  two  contributions  were  not  fully 
resolved  in  our  experiment  because  of  the  smallness  of  the 
.  xciton  binding  energy  ( £„  )  in  GaAs  and  the  high  tempera¬ 
ture  ( >  77  K )  of  the  experiment.  In  the  absence  of  strain  the 
heavy  hole  (hh)  and  light  hole  (Ih)  l^ands  are  degenerate  at 
the  Brillouin  zone  center  so  that  these  bands  contribute  only 
one  luminescence  peak.  The  split-off  band  is  0.34  eV  higher 


TABLE  I.  Summary  of  characterization  of  GaAs  films  and  sample  param- 
eteis  deduced  by  fitting  the  experimental  photoluminescence  spectra. 


Sample 

A 

GaAs/GaAs 

B 

GaAs/Si 

c 

GaAs/Si 

O 

GaAs/ Si 

Thickness  (pi ) 

Exciton 

linewidtn  (cm  " 1 1 

1.5 

3 

1.5 

1.5 

( 10 a  wide) 

d,  (heavy  hole) 

44 

53 

85 

A,  i  light  hole) 
Carrier 

i ; 

52 

'0 

no 

temperature  1  K.) 

SI 

127 

187 

216 

Stress  ( Itbar) 

0 

3.'  -0.5 

3.2  -0.4 

3  6  -  0,5 

A 

Carrier 

1.8 

1.9 

3.5 

4  8 

lifetime  (ps) 

:ooo 

79 

28 

18 

579 
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show  easily  that  for  a  GaAs  film  grown  on  ( 100)  Si  substrate 
this  stress  is  equivalent  to  the  sum  of  a  hydrostatic  pressure 
and  a  uniaxial  stress  perpendicular  to  the  film.  The  effect  of 
such  uniaxial  stress  on  the  optical  transitions  of  GaAs  can  be 
described  by  the  strain  Hamiltonian  proposed  by  Bir  et  al.14 
For  GaAs  subjected  to  a  tensile  stress  in  the  [  100]  direction, 
this  Hamiltonian  gives  the  following  stress-induced  band- 
gap  shifts  to  the  lowest  order  in  stress: 


X  3  V  dP  J 


=  —  10.4  meVAbar, 


(4) 


FIG  1.  Experimental  photoluminescence  spectra  of  GaAs  films  at  77  K 
(solid  curves).  A:  GaAs  on  GaAs  substrate.  B:  GaAs  onSi  substrate  (3  n 
thick).  C  GaAs  on  Si  substrate  ( 1 .5 /a  thick).  D :  GaAs  on  Si  substrate  ( 1.5 
fi  thick  and  10^  wide)  The  broken  curves  are  theoretical  fits  to  the  experi¬ 
mental  curves. 


in  energy  so  at  77  K  its  population  is  negligible.  To  simplify 
the  calculation,  we  have  neglected  the  discrete  excited  states 
of  the  excitons  and  the  excitonic  enhancement  of  the  contin¬ 
uum  band-to-band  transition. 1 1  The  error  introduced  by  this 
simplification  was  not  serious  for  GaAs  because  of  its  very 
small  exciton  binding  energy.  The  exciton  luminescence  is 
represented  by  a  Lorentzian  L(E)  centered  at  the  exciton 
energy  {Eg,  —  £„, )  while  the  band-to-band  transition  is  de¬ 
scribed  by  the  function  F(E) : 


L,  (£)  =  ]M,}2m!/2 


“  (£gi  -  Et*>)/k7‘ 

(£-£,,  +£„,)=- f-  Af  * 


(i) 


and 


£,(£)  =  ia/,:vH(v£/£‘' 


\)e~E/kT  E>Eg, 


In  Eqs.  ( 1 )  and  ( 2),  /refer  to  either  hh  or  Ih  bands,  IA/(  j:  are 
their  optical  transition  matrix  elements,  Egl  and  /a,  are  thetr 
band  gaps  and  joint-density -of-states  masses,  £„,  and  A,  are 
the  exciton  binding  energies  and  linewidths,  7"  is  the  carrier 
temperature,  and  k  is  the  Boltzmann  constant.  The  function 
El  E )  describes  an  asymmetric  peak  whose  high-energy  side 
is  determined  by  the  temperature  T  of  the  earners.  The  total 
PL  intensity  is  obtained  by  summing  the  contributions  from 
both  the  heavy  hole  and  light  hole  bands: 


/pi. 1  T  [AL,  (£)  f  £,(£)].  (3) 

i  -  hh 
Ih 

A  is  an  adjustable  parameter,  which  represents  the  relative 
intensity  of  the  lowest  exciton  level  to  the  band-to-band  tran¬ 
sition. 

An  external  hydrostatic  stress  changes  the  lattice  con¬ 
stant  and  shifts  the  energy  bands,  while  a  uniaxial  stress  can 
lift  some  of  band  degeneracies.  The  stress-induced  band-gap 
shifts  and  splittings  in  GaAs  have  been  studied  extensively 
by  PL  and  many  other  experiments. i:u  We  assume  that  the 
stress  distribution  inside  the  GaAs/Si  films  is  biaxial  with 
the  stress  perpendicular  to  the  film  equal  to  zero,  while  the 
stress  in  the  plane  of  the  GaAs  films  is  isotropic  One  can 


=  —  4.9  meVAbar, 

where  XE,  and  A Eh  are,  respectively,  the  lh  and  hh  band- 
gap  shifts,  5,,  and  Sn  are  components  of  the  compliance 
tensor  of  GaAs,  X  is  the  magnitude  of  the  isotropic  stress  in 
the  plane  of  the  films,  ( dEg/dP )  is  the  hydrostatic  pressure 
coefficient  of  the  band  gap,  and  b  is  the  deformation  poten¬ 
tial  for  uniaxial  stress  along  the  [  100]  direction.  To  fit  the 
PL  spectra  for  GaAs  films  grown  on  Si  substrate,  we  use  Eq. 
( 4 )  to  calculate  the  band  gaps  Eg,  under  stress  and  substitute 
these  values  into  Eq.  (3).  Figure  1  shows  the  experimental 
PL  spectra  (solid  lines)  and  the  corresponding  theoretical 
fits  (broken  curves).  In  fitting  curve  A  (GaAs/GaAs  sam¬ 
ple)  the  adjustable  parameters  are  A,  E„,  and  T and  their 
values  are  listed  in  Table  I.  Within  our  experimental  accura¬ 
cy  we  found  that  £„  =  0  meV  gave  the  best  overall  fit  to  the 
experimental  spectrum.  If  we  use  Etx  =  5  meV,  which  is  the 
generally  accepted  exciton  binding  energy  in  GaAs,  the 
theoretical  curves  show  a  dip  between  the  exciton  peak  and 
the  continuum  band-to-band  contribution.  The  dip  would 
not  be  present  if  the  higher  discrete  exciton  states  plus  the 
exciton  enhancement  to  the  continuum  density  of  states 
were  included  in  our  model.  Thus  the  use  of  the  adjustable 
parameter  A  and  the  assumption  that  £„  =  0  compensate 
for  these  omissions.  We  also  noted  that  in  sample  A  the  car¬ 
rier  temperature  T  was  very  close  to  the  lattice  temperature 
suggesting  that  laser  heating  of  the  samples  was  negligible. 
In  fitting  the  spectra  for  the  GaAs/Si  samples  we  again  as¬ 
sumed  £„,  =  0  for  both  hh  and  lh  transitions.  The 
adjustable  parameters  in  these  cases  are  A",  A,,  A*,  T,  and  A. 
The  results  are  shown  in  Table  I  also.  From  Eq.  (4)  the 
stress  X  causes  a  shift  in  the  “center  of  gravity”  of  the  hh  and 
lh  transitions  in  addition  to  splitting  them.  The  fact  that  a 
single  value  of  X  can  fit  both  the  shift  and  splitting  of  the  hh 
and  Ih  transitions  is  a  strong  indication  that  the  biaxial  stress 
we  have  assumed  is  basically  correct.  Furthermore,  the  fact 
that  X  deduced  from  the  different  GaAs/Si  samples  are 
more  or  less  the  same  is  consistent  with  its  thermal  origin. 

Besides  the  splitting  of  the  PL  peaks,  there  are  also  dif¬ 
ferences  in  line  shapes  and  relative  peak  intensities  between 
all  the  samples.  For  example  the  linewidths  of  the  exciton 
peaks  in  the  GaAs/GaAs  sample  were  considerably  narrow¬ 
er  than  those  of  the  GaAs/Si  samples.  The  broadening  of  the 
exciton  peaks  in  the  GaAs/Si  samples  can  result  from  either 
shorter  exciton  lifetimes  in  the  GaAs/Si  samples  due  to  dis- 
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locations  or  from  inhomogeneous  broadening  due  to  vari¬ 
ation  of  stress  in  the  sample.  We  will  show  later  that  the 
exciton  lifetime  in  our  GaAs/Si  samples  is  longer  than  tens 
of  picoseconds  ( ps )  so  homogeneous  lifetime  broadening  in 
these  samples  is  negligible.  Thus  to  reconcile  this  inhomo¬ 
geneous  stress  broadening  with  the  biaxial  -tress  distribu¬ 
tion.  we  have  to  assume  that  the  stress  is  locally  homogen¬ 
eous  and  isotropic  but  there  are  variations  over  an  area  of  the 
size  of  our  focal  spots  which  were  typically  larger  than  300 
fim  across.  Presumably  this  stress  variation  is  largest  near 
the  edge  of  the  film.  This  will  explain  the  larger  exciton 
linewidth  in  the  narrower  sample  D.  The  stress  variation 
over  the  film  can  be  mapped  out  by  a  tighter  focusing  of  the 
laser  spot  and  this  will  be  the  subject  of  future  investigations. 

From  the  relative  intensities  of  the  hh  and  lh  exciton 
peaks  and  the  line  shapes  of  the  higher  energy  tail  of  the 
band-to-band  contribution,  we  found  that  the  carrier  tem¬ 
peratures  in  the  GaAs/Si  samples  are  considerably  higher 
than  the  lattice  temperature.  (See  Table  I  for  the  carrier 
temperatures.)  We  propose  to  explain  this  “hot-carrier”  ef¬ 
fect  in  the  GaAs/Si  samples  by  the  shorter  lifetime  of  pho- 
toexcited  carriers  in  these  samples  due  to  the  higher  density 
of  dislocations.  To  correlate  the  carrier  lifetime  with  its  tem¬ 
perature  we  assume  that  the  photoexcited  electron-hole  ( e - 
h )  pair  lifetime  is  r.  The  rate  at  which  energy  is  transferred 
to  each  e-ft  pair  from  the  laser  beam  is  given  by 
H  =  (faj  —  Et  )/r  where  fuu  is  the  laser  photon  energy.  The 
rate  at  which  each  e-h  pair  will  lose  energy  due  to  emission  of 
longitudinal  optical  (LO)  phonons  is13 


where  &uLO  is  the  LO  phonon  energy,  r0  is  the  time  for  an  e- 
h  pair  to  emit  a  LO  phonon  and  is  ~0. 1  ps,16  and  TL  is  the 
lattice  temperature  assumed  to  be  77.4  K.  At  steady  state  we 
expect  H  =  L  and  from  this  equation  we  can  calculate  r  from 
the  experimentally  deduced  carrier  temperature.  The  resul¬ 
tant  r  are  shown  in  Table  I  for  samples  A-D.  The  shorter 
carrier  lifetimes  in  samples  C  and  D  as  compared  with  sam¬ 
ple  B  are  presumably  caused  by  their  smaller  thickness  so 
there  are  more  dislocations  at  the  surface  of  these  films.  The 
carrier  lifetimes  estimated  in  this  way  from  these  films  are 
somewhat  shorter  than  the  lifetime  of  photoexcited  free-car- 


riers  in  GaAs  containing  nonradiative  recombination 
centers  such  as  Cr.17  They  are  consistent  with  the  relative 
quantum  efficiencies  of  these  films  determined  with  the  same 
optical  setup.  The  results  of  these  quantum  efficiency  mea¬ 
surements  will  be  published  elsewhere.18 

in  nonclusion,  we  have  shown  that  from  the  photolumi¬ 
nescence  spectra  of  GaAs  films  grown  on  Si,  one  can  deter¬ 
mine  the  thermal  stress,  carrier  lifetime,  and  stress  inhomo¬ 
geneity  in  these  films. 
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Proposal  for  three-dimensional  internal  field  mapping  by  cw  electro-optic 
probing 
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The  cw  electro-optic  probing  technique  is  for  the  first  time  proposed  to  detect  the  three- 
dimensional  internal  field  distribution  in  linear  electro-optic  material  like  GaAs.  By  changing 
the  incident  angles  and  positions  of  the  probing  beam,  sufficient  information  of  the  electric 
field  distribution  is  included  in  the  phase  retardation  of  the  probing  beam.  If  certain  conditions 
on  the  probing  beam  are  satisfied,  a  very  simple  linear  relation  between  phase  retardation  and 
each  field  component  can  be  found  and  the  whole  problem  becomes  not  only  mathematically 
solvable  but  experimentally  feasible.  Finally,  a  three-dimensional  computer  simulation  is 
undertaken  to  illustrate  the  relation  between  field  distribution  and  detected  electro-optic 
signal. 


To  understand  the  physics  of  more  and  more  complicat¬ 
ed  III- V  compound  devices  and  to  find  out  proper  physical 
models  to  describe  or  predict  the  device  performance,  it  is 
essential  to  know  the  internal  potential,  field,  and  charge 
distributions.  However,  there  is  no  technique  presently 
available  to  detect  such  crucial  information  in  either  bulk 
semiconductor  or  devices.  For  those  compound  materials 
having  linear  electro-optic  properties  such  as  GaAs,  the 
phase  retardation  of  light  is  affected  by  the  electric  field  in¬ 
side  the  material,  and  thus  intrinsically  contains  the  infor¬ 
mation  of  electric  field  in  the  area  swept  by  the  probing  light. 
This  attractive  property  has  been  widely  applied  to  detect 
the  logic  state  and  transient  response  of  a  very  fast  GaAs 
circuit,'"5  or  in  cw  operation,  to  detect  the  surface  potential 
variation  and  the  boundary  of  a  p-n  junction.'’  7  However, 
none  of  the  work  has  ever  explored  the  possibility  of  the  great 
importance  of  extracting  three-dimensional  internal  field 
distribution  from  the  electro-optic  signals.  In  this  letter,  we 
propose  for  the  first  time  a  nondestructive,  high-resolution 
three-dimensional  field  measurement  method  by  means  of 
the  cw  electro-optic  probing  (CWEOP). 

Because  GaAs  is  the  most  extensively  used  electronic 
and  optic  material  having  linear  electro-optic  effect.  *t  will 
consider  it  as  an  example  in  this  paper.  The  index  ellipsoid  of 
GaAs  under  the  electric  field  is  described  in  Eq.  ( 1 ). 

(jr  J- y:  -+-  zf)/n\  +  2yix(Ex  yz  +  Ey  zx  4-  E.  xy)  =  1  , 

(1) 

where  y4l  is  the  only  nontrivial  electro-optic  coefficient  for 
the  cubic  crystal  GaAs.  If  the  wave  vector  k  of  the  probing 
laser  beam  has  directional  cosines  a,0.  y  with  respect  to  the 
three  crystal  axes,  the  plane  normal  to  k  and  passing  through 
the  origin  is  described  as  Eq.  (2). 

ax  -  0y  +  yz  =  0  .  (2) 

Geometrically,  those  points  satisfying  both  Eqs.  (  1 )  and  (2) 
form  an  ellipse  C  normal  to  k  as  shown  in  Fig.  1.  The  direc- 
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tions  of  the  major  and  minor  axes  of  the  ellipses  represent 
two  eigenvectors  of  the  polarization  states  of  the  light,  and 
their  magnitudes  are  equal  to  the  refraction  indices  of  the 
slow  and  fast  optical  axes.  In  spite  of  the  clarity  in  geometry, 
the  algebra  to  find  those  two  axes  is  very  involved.  However, 
if  the  k  is  in  the  y-z  plane  and  the  angle  with  the  z  axis  is  less 
than  45‘,  the  problem  can  be  greatly  simplified.  These  condi¬ 
tions  are  satisfied  for  all  the  following  derivations  and  pro¬ 
posed  measurements.  The  simplest  mathematical  approach 
is  to  first  project  the  ellipse  C  normal  to  k  into  an  ellipse  C '  in 
the  x-y  plane  as  shown  in  Fig.  1. 

After  obtaining  the  directions  and  magnitudes  of  the 
axes  of  C '  in  the  way  of  planar  analytic  geometry,  we  can 
easily  find  the  axes  of  ellipse  Cby  Eq.  ( 2 )  since  each  point  in 
C '  possesses  the  same  x  and y  coordinates  as  the  correspond¬ 
ing  point  in  C.  The  results  are  summarized  in  the  following 
set  of  equations: 

cot(20)  =£  =  {2a  4-  c)/d  ,  (3) 


a  =  (aJ)/(r)  -  2nlla/y)y,xEy  ,  (4a) 

c  -2n\y4X(aEy  4-  0Ex)/y-  {a1  4-  02)/yz,  (4b) 

d  =  2[  4-  /»o  r*i  ( Y^z  ~  <*EX  -  PEy  )/y] 

(4c) 

The  angle  <9  in  Eq.  ( 3 )  is  the  angle  between  the  x  axis  and  the 


its  projection  C '  in  the  x-y  plane,  t)  is  the  acute  angle  between  the  axis  I  major 
or  minor )  of  C ’  and  the  x  axis. 
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axis  (major  or  minor)  of  ellipse  C'  such  that  0  <  9<(ir/  2). 
The  directional  cosines  aJ3,y  are  those  of  k.  The  two  ex¬ 
trema  of  the  ellipse  normal  to  k  have  the  coordinates  de¬ 
scribed  as  follows: 


x  t  =  n{)  cos  6  v  ( 1  +  c‘ )  ±d'  , 
y  £  =  n()  sin  d 4 ( 1  +  c')  ±  d '  , 
z±  =  -  (ax  ,  +  Py  ±  )/y , 


and 


(5a) 

(5b) 

(5c) 


nx  =  +■ 


(6) 


where  c'  =  c/2  and  d'  =  d /[2  sin(2#)  ]  . 

After  finding  the  field  dependence  of  the  eigenstates  for 
arbitrary  incident  light,  we  look  for  the  relation  between  the 
detected  variation  and  spatial  distribution  of  electric  field  by 
the  prescribed  electro-optic  effect.  For  simplicity  but  with¬ 
out  loss  of  generality,  let  us  assume  the  probing  light  initially 
linearly  polarized  along  the  x  axis,  and  then  its  polarization 
state  can  be  represented  by  the  Jones  vector  as 


To  find  the  field  distribution  from  a  finite  number  of  probing 
beams,  we  need  to  set  up  the  mesh  points  in  the  tested  area. 
For  each  small  unit  of  those  meshes,  the  transformation  of 
the  polarization  state  from  input  beam  to  output  beam  is 
described  as  follows’’: 


T,  =s  (W,/c)(n„  -nt)  ,  (10) 

A  =  (<y/,/2o)  (n„  -!-«*).  ( 1 1 ) 


where  R,  is  the  matrix  rotating  the  initial  polarization  state, 
along  x  axis  in  this  case,  to  the  slow  axis  of  the  index  ellipse. 
The  quantity  T,  is  the  phase  retardation,  and  and  n„  are 
slow  and  fast  refraction  indices  of  the  ith  mesh,  and  /,  is  the 
length  of  the  optical  path  in  mesh  i.  Because  the  exp(  —  (A, ) 
term  in  Eq.  ( 8 )  does  not  have  any  physical  significance  but  a 
phase  shift,  we  will  omit  it  in  the  following  derivations. 

If  they  component  of  the  polarization  state  is  detected 
from  the  emerging  beam,  the  final  Jyk  is  obtained  from  a 
series  of  multiplication  of  the  transformation  matrices  IV,  in 
Eq-  (8)  along  the  light  path.” 


=-y-^  T,  sin(2d>, )  .  (12) 

where  the  summation  is  over  all  the  mesh  points  passed  by 


In  general,  both  T,  and  d>,  in  Eq.  (12)  are  functions  of£„, 
£y,,and£.,  resulting  in  a  number  of  cross  terms  in  Eq.  (12). 
To  simplify  the  analysis,  we  must  choose  the  incident  beam 
to  be  either  iny-z  plane  or  in  x-z  plane  (i.e.,a  =  0or£  =  0). 
Assuming  k  in  y-z  plane  (i.e.,  a  =  0)  as  an  example,  we 
simply  set  4>,  in  Eqs.  (9)  and  ( 12)  equal  to  d ,  in  Eq.  (3)  \fd 
in  Eq.  (4c)  >0  or  equal  to  6,  +  jt/2  ifd<0.  The  ir/2  differ¬ 
ence  between  4>,  and  6,  as  d  <  0  is  due  to  the  previous  defini¬ 
tion  that  the  angle  <J>,  is  between  x  axis  and  the  slow  axis  but 
6,  is  confined  in  (0,ir/2)  to  ensure  sin  ( 29, )  in  Eq.  (6)  larger 
than  zero.  By  combining  Eqs.  (3),  (4),  (6),  (10),  and  (12) 
and  the  above  relation  between  9,  and  <t>,  with  the  condition 
a=0,Jyk  can  be  explicitly  represented  in  terms  of  the  field 
components  as  follows: 


K  =  -  i(a>lun30yi,)/(2c)  . 


(14) 


The  field-dependent  nonlinear  term  in  the  square  root  of 
Eq.  ( 1 3 )  is  negligible  as  long  as  the  incident  beam  has  a  small 
angle  with  the  z  axis,  in  other  words,  0  '/y3  <  1 .  Therefore, 
Eq.  (13)  can  be  approximated  as 


A.  =  | <■=> 

Equation  ( 15)  is  the  equation  setting  up  the  linear  relation 
between  the  field  components  and  the  detected  probing 
beam  from  different  angles  and  positions.  The  Ex,  in  Eq. 
( 15)  is  missing  because  of  the  condition  a  =  0.  To  find  £„, 
we  can  replace  0Ey,  with  aE„  in  Eq.  ( 1 5 )  assuming  that  the 
light  is  in  x-z  plane.  The  signal  from  a  photodetector  is  pro¬ 
portional  to  the  magnitude  of  Jyk  based  on  the  same  princi¬ 
ple  as  the  amplitude  modulation  of  electro-optic  material.9 
Therefore,  the  three-dimensional  electric  field  distribution 
measurement  becomes  not  only  theoretically  possible  but 
practically  feasible  by  means  of  the  recently  developed  cw 
electro-optic  probing  technique. ft-7  From  the  measured  field 
distribution,  the  real  values  of  potential  profile,  field  intensi¬ 
ty,  and  charge  density  can  be  substantially  obtained  since  the 
applied  voltage  supplies  the  necessary  boundary  conditions. 

Figure  2(c)  shows  the  computer  simulation  results  of 
the  electro-optic  signals  from  Eq.  (15)  under  some  given  Ey, 
E.  distributions  shown  in  Figs.  2(a)  and  2(b).  The  detected 
area  in  Fig.  2  is  100  X  100  /im;  in  the  y-z  plane  and  the  inci¬ 
dent  probing  beam  scans  from  —  1 8*  to  18°  with  1°  difference 
between  adjacent  beams.  Not  only  does  the  beam  incident 
angle  vary,  the  probing  beam  also  moves  along  the  y  axis 
with  0. 1 8  urn  per  step  to  supply  enough  linearly  independent 
equations  to  solve  the  field  components  at  each  mesh  point. 
A  larger  mesh  size  can  be  used  in  those  regions  of  slow- 
varying  field  to  reduce  the  number  of  probing  beams  and 
save  computer  memories  without  sacrificing  the  accuracy. 
Finally,  a  simple  experiment  is  done  to  justify  our  proposed 
method  in  practice.  As  shown  in  Fig.  3,  the  electro-optic 
signals  measured  under  normal  incidence  and  1 .4°  tilted  inci- 


the  probing  beam  with  wave  vector  k.  In  Eq.  ( 12),  we  have  dence  are  clearly  distinguishable.  Those  features  occurring 
approximated  exp(  —  tT, )  =  1  -  iT,  since  the  phase  retar-  in  the  tilted  incidence  curve  are  partly  due  to  the  different 
dation  in  each  mesh  is  always  much  smaller  than  one  radian.  optical  path  from  the  normal  incident  one  and  partly  due  to 


FIG  2.  Computer  simulation  of  the  electro-optic  signal  under  the  given 
field  profile  in y-2  plane.  Ev  and  E.  profiles  are  shown  in  I  a!  and  i  bl  but  E , 
is  not  shown  here  because  it  does  not  affect  the  electro-optic  signal  if  the 
light  has  the  wave  vector  in  the y-z  plane;  (c)  shows  the  signal  variation  as 
both  the  mctdent  angle  and  the  incident  position  along  the  v  axis  change 


FIG.  2.  Measured  electro-optic  signal  across  two  electrodes  on  GaAs  sub¬ 
strate.  The  dark  curve  is  measured  under  normal  incidence  and  the  light  one 
is  measured  under  incidence  at  an  angle.  Lock-in  technique  is  used  in  this 
experiment. 

the  influence  of  lateral  field.  The  more  complete  experiment 
and  the  detailed  analysis  of  the  experimental  data  will  be 
published  elsewhere. 

In  conclusion,  we  propose  a  very  new  method  to  nonde¬ 
structive^  detect  the  three-dimensional  field  distribution  in 
linearly  electro-optic  materials  such  as  GaAs,  and  math¬ 
ematically  formulate  it.  If  an  appropriate  incident  angle  is 
chosen,  the  electro-optic  signal  is  a  linear  combination  of  the 
field  components  in  the  light  path.  This  property  facilitates 
the  data  analysis  and  makes  the  experiment  more  feasible. 
To  obtain  more  accurate  results,  the  finite  beam  spot  size, 
multiple  reflection  inside  the  material,  reflection  at  the  inter¬ 
face,  as  well  as  the  light  shape  function  all  need  to  be  taken 
into  account. 

The  authors  are  grateful  to  J.  Y.  Fang  for  his  generous 
help  in  computer  graphics.  This  work  is  supported  in  part  by 
Bell  Communications  Research  and  by  Joint  Services  Elec¬ 
tronics  Program  AFOSR  F49620-84-C-0057. 
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We  have  successfully  measured,  for  the  first  time,  standing  waves  in  a  GaAs  coplanar 
waveguide  at  frequencies  of  8.2107  and  12.310  GHz  by  using  harmonic-mixing  electro-optic 
probing.  The  technique  is  nondestructive  and  has  great  potential  in  expanding  the  measuring 
frequency  band  to  millimeter  wave.  This  letter  describes  the  principle  of  the  technique,  the 
experimental  setup,  and  the  measurement  results. 


In  recent  years,  many  exciting  applications  of  the  linear 
electro-optic  effect  as  a  measurement  technique  have  been 
explored.  The  electro-optic  sampling  technique  has  been  ap¬ 
plied  to  detect  the  transients  of  ultrafast  electronic  and  op¬ 
toelectronic  devices.1’  More  recently,  cw  electro-optic  prob¬ 
ing  (CWEOP)  has  been  used  to  measure  the  potential 
profile,  field,  and  charge  distribution  in  GaAs  material  and 
devices.’-5 

In  GaAs  microwave  integrated  circuits,  a  standing- 
wave  electric  field  will  be  established  in  a  waveguide  if  the 
impedance  of  the  load  is  mismatched  to  the  characteristic 
impedance  of  the  waveguide.  Therefore,  in  principle,  the 
electro-optic  probing  technique  can  be  applied  to  detect  the 
distribution  of  standing-wave  fields  in  GaAs  microwave 
waveguides. 

Based  on  this  idea,  we  have  successfully  measured,  for 
the  first  time,  the  distribution  of  a  standing-wave  field 
formed  in  a  GaAs  waveguide  using  a  new  technique,  har¬ 
monic-mixing  electro-optic  probing.  The  principle  of  the 
technique,  the  experimental  setup,  and  results  are  presented 
in  this  letter. 

Conventional  slotted  line  techniques  become  impracti¬ 
cal  at  high  microwave  frequencies  in  waveguides.  No  tech¬ 
nique  exists  for  monolithic  microwave  transmission  lines. 
Here  we  employ  the  harmonic-mixing  electro-optic  probing 
to  overcome  these  difficulties.  Briefly,  when  a  laser  pulse 
train  from  a  mode-locked  Y AG  laser,  modulated  by  an  elec¬ 
tric  field,  illuminates  a  photodetector,  the  photocurrent  pro¬ 
duced  by  the  photodetector  is  proportional  to  the  product  of 
the  light  intensity  and  the  modulating  voluge  The  modula¬ 
tor  (GaAs  coplanar  waveguide  sample  m  ’ur  .use  i -photo¬ 
detector  combination  can  be  viewed  as  a  ~r.nct  In  the  fre¬ 
quency  domain,  any  signal  at  frequencv  *.  m  -raeating  in 
vhe  waveguide  will  mix  wirh  all  harmonic'  t  ■‘■c  •undamen- 
tal  repetition  frequency  f,  of  the  mode  .voo  as <r  pulse. 
Sidebands  will  appear  at  frequency  nr..  -  ‘  mixing 
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wym  y  (  sin(fl7i/0r)  \ 
V,  ,t-A  mfQT  ) 


X  [sin  2ir(nf0  +f„  )r  -  sin  2v(nf,  -fm  )t  ],  ( 1 ) 

where  /„,  is  the  average  photocurrent,  Vm  is  the  amplitude 
of  the  modulating  voltage,  VT  is  the  half-wave  voltage  of  the 
GaAs  modulator,  n  is  the  order  of  harmonic  component,  and 
r  is  the  pulse  width  of  the  mode-locked  laser.  Insofar  as  the 
mode-locked  laser  as  a  local  oscillator  is  concerned,  the  only 
requirement  is  to  select  a  proper  frequency  fm  of  rf  signal, 
then  we  can  easily  measure  the  mixing  electro-optic  signal  at 
an  intermediate  frequency  f ,  which  is  given  by 

(2) 

where  is  the  electro-optic  peak  voltage  developed  across  a 
load  resistor  RL,  and  is  proportional  to  the  modulating  vol¬ 
tage  vm- 

The  experimental  setup  is  shown  in  Fig.  1.  A  laser  pulse 
train  of  about  10  ps  duration  at  a  repetition  rate  of/,  =  82 
MHz  is  generated  by  the  compressed  output  of  a  cw  mode- 
locked  Nd:YAG  laser  ( A  =  1.06  ^tm).  The  quarter-wave 
plate  introduces  a  phase  shift  of  ~/2  between  the  “fast”  and 
“slow”  components  of  light,  so  that  the  electro-optic  modu¬ 
lator  (i.e..  the  GaAs  waveguide  sample)  operates  in  the  lin¬ 
ear  region.  The  laser  beam  passes  through  the  GaAs  sample 
and  is  focused  to  a  diameter  of  about  8  /zm  on  the  central 
electrode  of  the  coplanar  waveguide.  The  widths  of  the  cen¬ 
tral  electrode,  side  electrodes,  and  the  spacing  between  them 
are  135,  590  and  65  ptm,  respectively.  The  rf  output  with 
frequency  fm  of  an  oscillator  is  fed  to  the  GaAs  coplanar 
waveguide  (CPW)  via  a  connector  (type  SiMA),  and  the 
other  terminal  of  the  waveguide  is  open,  so  a  standing-wave 
field  is  established  along  the  waveguide.  The  laser  beam 
modulated  by  the  standing-wave  field  is  reflected  from  the 
electrode,  then  passes  through  an  analyzer  to  convert  the 
phase  change  into  amplitude  change.  A  low-speed  germani¬ 
um  photodetector  is  used  to  detect  the  mixing  signal  at  the 
intermediate  frequency  f.  A  preamplifier  is  inserted 
between  the  photodetector  and  a  spectrum  analyzer  to  opti¬ 
mize  the  5  /.V  ratio.  A  TV  camera  is  used  to  monitor  the 
measuring  position  exactly.  By  scanning  the  laser  beam  from 
the  open-circuit  terminal  to  the  other  terminal  along  the  lon- 
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gitudinal  <011)  direction  (zaxis  in  Fig.  1 )  of  the  CPW  elec¬ 
trode,  we  can  detect  the  standing-wave  pattern. 

Figure  2  shows  the  harmonic-mixing  electro-optic  sig¬ 
nal  displayed  on  the  spectrum  analyzer  by  mixing  the  rf  sig¬ 
nal  (fm  =  8.2107  GHz)  with  the  100th  harmonic  compo¬ 
nent  ( n  =  100)  of  the  fundamental  repetition  rate  (/Q  =  82 
MHz)  of  the  mode-locked  pulse.  The  intermediate  frequen¬ 
cy/  of  the  mixer  signal  is  equal  to  10.7  MHz.  The  corre¬ 
sponding  standing-wave  pattern  measured  is  given  in  Fig.  3. 
The  zero  of  the  abscissa  corresponds  to  the  point  where  the 
measurement  began,  but  it  should  be  pointed  out  that  the 
real  zero  point  corresponding  to  the  edge  of  the  open  termi¬ 
nal  of  CPW  locates  at  about  45  mil  (  xz  1 . 1 4  mm )  to  the  left  of 
the  zero  of  the  abscissa.  From  the  pattern  we  can  see  that  the 
distance  A,w  between  valleys  is  about  270  mil  (  ~6.9  mm) 
which  is  equal  to  the  half-wavelength  of  the  rf  signal.  Substi¬ 
tuting  these  data  into  the  following  well-known  formula  for 
evaluating  the  effective  index  n,  of  GaAs  CPW, 

n,  =c/2A,w/m,  (3) 

where  c  is  the  speed  of  light  in  air.  we  get  n,  =  2.65,  which  is 
in  good  agreement  with  the  value  used  in  the  design  of  micro- 


wave  GaAs  integrated  circuits.  Similar  measurement  at  the 
rf  frequency  of  12.310  GHz  (by  mixing  with  the  150th  har¬ 
monic  component  of  the  fundamental  repetition  rate/,  =  82 
MHz  of  mode-locked  pulse)  yields  Alw  =  180  mil  C  =  4.6 
mm).  This  is  in  excellent  agreement  with  the  value  ofA,„  as 
evaluated  using  n,  =  2.65.  From  the  measured  standing- 
wave  pattern,  we  can  evaluate  the  voltage  standing-wave  ra¬ 
tio  ( VSWR)  p  and  the  reflection  coefficient  F  of  the  GaAs 
CPW,  the  result  is  as  follows:  fm  =  8.2107  GHz,  p  =  5.97, 
r  =  0.71  ;fm  =  12.310  GHz,  p  =  7.94,  T  =  0.78.  The  rea¬ 
son  for  p^  oo  can  be  attributed  to  fringing  field  capacitance 
at  the  end  of  the  transmission  line  and  electrode  losses. 

In  summary,  we  have  successfully  measured  the  stand¬ 
ing-wave  pattern  in  a  GaAs  coplanar  waveguide  at  frequen¬ 
cies  of  8.2107  and  12.310  GHz  by  using  harmonic-mixing 
electro-optic  probing.  The  experimental  results  are  in  agree¬ 
ment  with  the  theory  and  the  reproducibility  of  the  measure¬ 
ment  is  excellent.  The  technique  is  reliable  and  noninvasive. 
Since  the  input  bandwidth  of  this  technique  easily  exceeds 
100  GHz  for  picosecond  mode-locked  pulses,  we  believe  in 
the  near  future  the  technique  can  be  developed  to  become  a 
standard  in  making  noncontact  on-chip  measurement  of 
standing  waves  in  GaAs  millimeter  wave  integrated  circuits. 
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Resonant  tunneling  diodes  with  AlAs  barrier:  Guides  for  improving 
room-temperature  operation 
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In  this  communication  experimental  results  of  AlAs/GaAs/AlAs  resonant  tunneling  diodes 
are  compared  to  theory,  and  the  effect  of  higher  resonant  states  on  the  room-temperature 
resonant  tunneling  current  peak-valley  ratio  is  studied.  It  is  shown  that  for  a  well  thickness  of 
70  A,  the  leakage  current  through  higher  resonant  states  of  the  well  can  reduce  the  peak-to- 
valley  current  ratio  by  more  than  400%  when  compared  to  a  similar  diode  of  50-A  well 
thickness.  Through  these  studies  a  summary  of  suggestions  is  listed,  which  is  intended  to 
improve  the  room-temperature  operation  of  resonant  tunneling  diodes. 


The  idea  of  using  the  semiquantized  states  in  ultrathin 
heterostructures  for  obtaining  negative  differential  resis¬ 
tance  was  first  proposed  by  Esaki  and  Tsu1  and  developed  by 
other  workers2  as  early  as  1970.  As  the  molecular-beam  epi¬ 
taxy  and  other  epitaxial  growth  techniques  developed,  re¬ 
newed  interest  appeared  in  the  field,  and  theoretical  and  ex¬ 
perimental  papers  aiming  at  obtaining  room-temperature 
observation  of  negative  differential  resistance  were  pub¬ 
lished.5"5  The  observation  of  Sollner  era/.,4  who  have  report¬ 
ed  2.5-THz  operation  of  resonant  tunneling  diodes  at  77  K, 
has  proven  the  potential  of  these  diodes.  Also,  the  recent 
results  of  Tsuchiya  et  a/.5-4  show  that  the  material  quality  of 
resonant  tunneling  diodes  has  improved  enough  to  allow  the 
observation  of  peak- valley  current  ratio  ( P/V)  of  3  at  296  K. 
Comparing  the  results  of  Tsuchiya  et  al.  to  other  results  of 
AlAs  barrier  resonant  tunneling  diodes7  shows  that  design 
and  material  quality,  and  not  phonon  scattering,  are  the  fac¬ 
tors  limiting  the  peak-valley  ratio  at  room  temperature.  The 
potential  large  peak-valley  current  ratio  and  ultrahigh-fre¬ 
quency  characteristic  of  these  diodes  at  room  temperature 
make  them  ideal  for  high-frequency  generation  and  detec¬ 
tion  systems.  However,  to  realize  the  maximum  capability  of 
the  resonant  tunneling  diodes,  a  systematic  study  of  the  pa¬ 
rameters  affecting  the  current-voltage  characteristic  of  these 
devices  is  necessary.  These  parameters  include  the  effective 
masses  in  the  barrier  and  well,  the  thickness  of  the  well  re¬ 
gion,  and  the  voltage  at  which  the  current  resonance  peak 
occurs.  In  this  paper  the  significance  of  these  parameters  is 
analyzed,  and  conclusions  are  drawn  which  discuss  the 
choice  of  parameters  that  can  improve  the  peak-valley  cur¬ 
rent  ratio  at  room  temperature. 


The  structure  under  consideration  is  shown  in  Fig.  1 .  In 
the  conventional  approach  to  the  calculation  of  resonant 
tunneling  current,  the  Schrodinger  equation  is  solved  in  the 
barrier,  and  the  well  and  the  solutions  are  matched  at  the 
boundary  to  satisfy  the  continuity  requirements.  Then  the 
contribution  to  the  current  from  each  solution  is  weighted  by 
the  Fermi-Dirac  function  and  integrated  to  obtain  the  total 
tunneling  current.  This  procedure  is  very  difficult  to  be  car¬ 
ried  out  self-consistently;  hence  approximations  must  be 
made  to  simplify  the  process.  In  the  approach  used  in  this 
paper,  the  following  two  points  have  been  considered  in  ob¬ 
taining  the  tunneling  current.  First,  in  solving  the  Schro- 
dinger  equation,  the  Wentzel-Kramers-Brillouin  (WKB) 
approximation  has  been  used.  This  approximation  allows  us 
to  incorporate  the  effect  of  electric  field  on  the  solution  of 
Schrodinger  wave  equation.  Second,  to  obtain  a  more  accu- 
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FIG.  I .  General  resonant  tunneling  diode  under  consideration.  The  quanti¬ 
tation  of  electron  energy  in  region  1  is  not  espected  to  influence  our  analysis 
substantially. 
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FIG.  2.  Low-temperature  comparison  of  Tsuchiya  et  al  \  ( Ref.  4 )  experi¬ 
ment  with  current  contribution  from  T  J.JC  valleys.  Note  that  the  f -valley 
current  alone  cannot  give  agreement  with  experiment.  The  dashed  line  in 
the  experimental  curve  i«  due  to  finite  ohmic  contact  resistance.  The  device 
parameters  are  d,  =»  d,  <*  22.fi  A.  d„  «  70  A.  temperature  =  77  K.  and 
W=SxlO,,cm'1. 


rate  picture  of  the  distribution  of  electrons  in  the  conduction 
band,  a  self-consistent  energy-band-bending  model,  taking 
into  account  the  voltage  drop  due  to  the  accumulation  of 
electrons  in  region  1  and  depletion  of  electrons  from  region  2 
of  Fig.  1,  is  used.  The  model  can  be  summarized  as  follows 
(GaAs  base  materials  only): 


V> 


(V)  •  (Efl  +  eV ,  +  Efr)/e  , 


(eV)  =2.65x10' 


(d t  +  d2  -+-  tf, 


-r- 


in 

(2) 


r(eV)  «313JV/  - VJ--  -V  , 

V  (d ,  +  di  +  dw)  ) 


where  Nd  (cm  J)  is  the  doping  outside  the  barrier,  d2, 
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FIG.  3.  Room-temperature  comparison  of  experimental  1  Tmchiva  et  at. 
(Ref.4)|  and  theoretical  resonant  tunneling  currents  The  general  behavior 
between  the  two  curves  agree  well.  The  parameters  used  are  those  of  X  val¬ 
leys.  as  indicated  in  Fig.  2. 


FIG.  4.  Current-voltage  characteristic  of  resonant  tunneling  diode  for  three 
well  thicknesses.  As  the  well  thickness  decrease,  the  room-temperature 
peak-vailey  ratio  increases,  and  the  voltage  at  the  current  peak  (  V, )  also 
increases.  The  other  parameters  are  the  same  as  those  of  Fig.  3. 


and  du  are  in  meters,  and  the  quantities  Efl  ,EmV%,dvdv 
and  d„  are  defined  in  Fig.  1.  It  should  be  noted  that  Efi  is 
calculated  by  linearly  approximating  the  conduction-band 
energy  in  region  1.  and  Err  is  obtained  by  the  usual  depletion 
approximations  in  region  2. 

With  the  use  of  the  above-mentioned  consideration  and 
approximations,  the  effect  of  parameters  on  the  tunneling 
current  will  not  be  discussed.  To  analyze  the  importance  of 
the  effective  mass  of  electrons  in  an  AlAs  barrier,  we  have  to 
decide  which  band  minima  to  use.  AlAs  is  an  indirect  band- 
gap  material,  whereas  GaAs  is  a  direct  one.  Thus,  to  analyze 
the  tunneling  current  contribution  of  each  conduction-band 
minima  ( T  JLJC)  of  AlAs  barrier  to  the  overall  current,  three 
effective  masses  (mf  are  assigned  to  the  corre¬ 

sponding  conduction-band  minima  of  the  AlAs  barrier  re¬ 
gion.  Then  each  band-minimum  state  of  GaAs  is  separately 
matched  to  the  r -band  minima  of  the  GaAs  regions.  The 
corresponding  resonant  tunneling  currents  and  experimen¬ 
tal  results  of  Tsuchiya  et  al*  are  plotted  in  Fig.  2.  Different 
axes  are  used  for  theoretical  and  experimental  values  to  ac¬ 
count  for  series  ohmic  resistance.  As  can  be  seen,  it  is  no 
longer  valid  to  assign  the  T  -band  effective  mass  of  GaAs  to 
that  of  AlAs  barrier.  By  using  the  effective  masses  given  in 
Adachi’s  review  article,8  it  is  concluded  that  the  T  valley  of 
AlAs  alone  cannot  account  for  the  observed  current-density 
peak  of  2  kA/cm2  of  Tsuchiya  et  al.’s*  experiment.  A  linear 
combination  of  all  three  band  minima  are  needed  to  account 
for  the  resonant  tunneling  current.  Note  also  that  the  cur¬ 
rent  is  a  sensitive  function  of  effective  mass  and  barrier 
height.  Without  a  detailed  knowledge  of  the  real  and  com¬ 
plex  Bloch  wave  function  and  energy-crystal  momentum- 
dispersion  relation  of  AlAs  barrier,  it  would  be  useless  to 
calculate  the  magnitude  of  total  resonant  tunneling  current. 

However,  for  studying  the  room-temperature  behavior 
of  resonant  tunneling  diodes,  we  do  not  need  to  go  into  this 
detail,  the  tunneling  current  due  to  theT  valley  will  suffice  as 
our  qualitative  framework.  This  choice  is  not  a  bad  one, 
since  X  minima  of  AlAs  represents  the  lowest  energy  barrier 


3475 


J.  Appl.  Ptrya.,  Vol.  62,  No.  8, 1 5  October  1 987 


0.  Vakhahoori  and  S.  Wang 


3475 


( ~  0. 57  eV )  to  the  electrons  of  the  GaAs  well.  From  conven¬ 
tional  quantum-mechanical  calculations,  it  is  known  that 
the  smaller  the  barriers  of  quantum  well,  the  less  is  the  sepa¬ 
ration  of  the  first  and  second  quantized  energy  levels.  This 
smaller  energy  separation  will  in  turn  contribute  to  higher 
unwanted  tunneling  current  through  a  second  resonant  state 
for  the  X  valley  compared  to  the  T  or  L  valley,  especially  at 
room  temperature  where  some  electrons  are  energetic 
enough  to  tunnel  through  the  second  resonant  state.  Figure  3 
compares  the  qualitative  behavior  of  calculated  Af-valley  res¬ 
onant  tunneling  current  to  that  of  Tsuchiya  et  a/.’s  experi¬ 
ment  at  296  K.  As  can  be  seen,  the  peak-valley  current  ratio 
(P/V)  of  experiment  and  theory  agree  well.  The  thermionic 
current  was  also  included,  but  the  effect  is  negligible. 

To  decrease  the  unwanted  current  through  the  second 
resonant  state  and  thus  to  increase  the  peak-valley  current 
ratio,  the  well  region  should  be  decreased  in  thickness.  Fig¬ 
ure  4  shows  the  current-voltage  characteristics  of  these  di¬ 
odes  for  different  well  thickness  at  room  temperature.  It  can 
be  seen  that  the  peak-to-valley  ratio  can  easily  be  improved 
by  a  factor  of  4  by  decreasing  the  well  thickness  from  70  to  50 
A.  The  result  shown  in  Fig.  4  is  very  encouraging  for  room- 
temperature  operation  of  resonant  tunneling  diodes  of 
AlAs/GaAs/AlAs  type.  However,  one  should  not  expect  to 
observe  such  a  large  improvement  in  practice.  Not  all  un¬ 
wanted  leakage  currents  are  due  to  higher  resonance  states; 
the  leakage  current  due  to  phonon-assisted  processes  cur¬ 
rents  due  to  deep  levels  and  surface  states,  and  conduction 
through  the  L  valley  of  the  GaAs  well  should  also  be  consid¬ 
ered.  The  other  mechanisms  can  be  considered  as  an  effec¬ 
tive  shunt  conductance.  As  the  thickness  of  the  well  region  is 
decreased,  the  voltage  at  which  current  peak  occurs  ( Vr ) 
will  increase  (Fig.  4).  Also  the  leakage  current  through  the 
above-mentioned  effective  shunt  conductance  increases  as 
the  applied  voltage  is  increased.  The  increase  of  leakage  cur¬ 
rent  as  the  well  thickness  is  decreased  is  also  observed  experi¬ 
mentally.  Tsuchiya  and  Sakaki9  observed  that  the  valley  cur¬ 
rent,  which  is  a  measure  of  leakage  current,  increases 
superlinearly  as  the  well  thickness  decreases.  Thus,  as  the 
well  thickness  decreases,  the  leakage  current  at  the  bias  vol¬ 


tage  required  for  resonance  ( V, )  increases.  Therefore,  there 
is  an  optimum  well  thickness  at  which  the  leakage  current 
due  to  other  phenomena  will  become  more  important  than 
the  unwanted  current  due  to  second  resonant  state.  This 
thickness  is  expected  to  be  less  than  60  A. 

To  summarize  the  discussion  above  for  enhancing  the 
nonlinear  characteristic  of  resonant  tunneling  diodes,  the 
following  conclusions  are  made:  (a)  To  build  up  the  reso¬ 
nant  state,  the  electron  scattering  should  be  small.  This 
means  that  the  interface  quality  should  be  good.  To  get 
abrupt,  homogeneous  interfaces,  the  method  of  interrupted 
growth,10  which  has  resulted  in  one  of  the  best-reported 
room-temperature  peak- valley  ratios,  is  recommended,  (b) 
Use  AlAs  instead  of  AlxGa,  _*As  as  the  barrier.  Higher 
barrier  means  more  quantization  of  states  in  the  well,  which 
can  reduce  unwanted  current  leakages,  (c)  Use  well  thick¬ 
ness  of  less  than  70  A  to  maximize  peak- valley  current  ratio. 
Consideration  of  the  above-mentioned  points  is  expected  to 
be  of  use  in  designing  and  growing  resonant  tunneling  diodes 
of  large  nonlinear  current- voltage  characteristics  to  be  oper¬ 
ated  at  room  temperature. 

This  work  is  supported  in  pan  by  Joint  Services  Elec¬ 
tronics  Program  under  Contract  No.  F49620-84-C-0057 
and  by  National  Science  Foundation  under  Grant  No.  ECS- 
8410838. 
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A  novel  surface  emission  of  coherently  generated  second-harmonic  wave  is  reported  for  the 
first  time. The  technique  is  used  for  the  observation  of  the  difference  in  propagation  constant  of 
the  TEu  and  TM0  mode  of  GaAs/AlGaAs  waveguide  cavity  to  a  high  degree  of  accuracy.  In 
this  technique  the  second-harmonic  signal  propagates  out  from  the  top  surface  of  the 
waveguide  structure,  converting  the  modal  phase  difference  between  TEo  and  TM0  modes  into 
intensity  variation  along  the  waveguide  length.  The  second-harmonic  signal  is  easily  observable 
by  the  naked  eye,  and  the  technique  does  not  require  wavelength  tuning  or  mechanical 
movement  for  the  measurement  of  birefringence. 


In  this  letter  we  report  phase-matched,  coherent,  sec-  MHz  and  pulse  width  of  lOOps  was  focused  by  a  40  X  objec- 


ond-harmonic  generation  of  a  mode-locked  Nd:YAG  laser 
beam  propagating  pependicular  to  the  direction  of  flow  of 
the  fundamental  beam.  In  conventional  second-harmonic- 
generation  techniques,  the  amount  of  second  harmonic 
propagating  out  perpendicular  to  the  direction  of  a  forward 
propagating  fundamental  beam  is  negligible,  since  the  phase¬ 
matching  condition  is  violated.  However,  in  the  technique 
presented  here,  a  standing  wave  of  second-harmonic  polar¬ 
ization  is  generated  along  the  waveguide  cavity.  The  second- 
harmonic  signal  generated  this  way  is  show  to  be  coherent 
and  to  satisfy  the  wave  vector  conservation.  To  the  best  of 
the  authors'  knowledge,  this  is  the  first  time  the  second-har¬ 
monic  wave  has  been  generated  in  this  manner  and  used  to 
study  waveguide  parameters. 

The  parameter  studied  in  this  experiment  is  the  modal 
birefringence  of  a  waveguide  structure.  The  difference  in 
propagation  constants  of  the  TE„  and  TM0  mode  of  a  GaAs/ 
AlGaAs  waveguide  is  of  paramount  importance  for  an  effi¬ 
cient  electro-optic  waveguide  TE-— TM  mode  conversion1 
and  wavelength  filtering.2  Waveguide  polarization  conver¬ 
sion  is  a  basic  signal  processing  function  and  will  be  useful  in 
realizing  single-mode  communication  systems.  In  addition 
to  its  integrability,  the  TE*-»TM  mode  convertor  together 
with  integrated  polarizer1  can  be  used  as  a  high-speed  optical 
modulator.  TE*-*TM  mode  convertor  can  be  realized  as  a 
broad-band  or  narrow-band  (full  width  at  half-maximum 
-4.5  A):  polarization  convenor.  Thus  TE— TM  conver¬ 
tors  can  be  used  as  a  polarization  convenor,  an  optical  mod¬ 
ulator,  or  a  wavelength  filter  in  future  opto-electronic  cir¬ 
cuits.  However,  to  realize  any  of  the  aoove  devices,  an  exact 
knowledge  of  modal  phase  difference  A/?  =  0TE  —  0TM  is 
required.  The  observation  of  surface  emission  of  the  second- 
harmonic  signal  gives  the  modal  phase  difference  accurately 
(error  S  1%),  directly,  and  without  any  ambiguity.  Fur¬ 
thermore,  the  requirement  that  the  second-harmonic  signal 
is  not  absorbed  strongly  by  the  top  cladding  layer  of  the 
waveguide  is  not  severe,  since  some  sommon  electro-optic 
polarization  modulators  on  GaAs  are  made  with  transpar¬ 
ent  top  cladding  layer. 1 

The  experimental  setup  used  is  shown  in  Fig.  1 .  Abut  20 
mW  average  power  from  a  mode-locked  NdtYAG 
(A  =  1.064  ,um)  laser  output  with  a  repetition  rate  of  82 


tive  lens  on  a  facet  of  a  waveguide  cavity  1  mm  long.  The 
waveguide  structure  is  shown  in  the  inset  of  Fig.  1.  The  out¬ 
put  power  from  the  other  facet  was  measured  to  be  —  1  mW. 
The  second-harmonic  generation  was  monitored  by  looking 
into  a  low- magnification  microscope  and  optimized  by 
changing  the  polarization  of  the  incident  fundamental  beam 
by  using  a  half-wave  plate  to  equalize  the  power  coupled  into 
TE0  and  TMq  mode  of  the  waveguide.  Once  the  green  light 
( 532  nm)  was  optimized  by  eye,  the  image  of  the  near  field 
was  recorded  on  a  1600  ASA  Fuji  film  exposed  for  3  min. 
The  picture  was  taken  by  using  a  commercial  camera  and 
looking  into  the  eye  piece  of  the  microscope.  The  picture  is 
shown  in  Fig.  2.  There  are  1 6  periods  of  intensity  modulation 
in  1  mm  corresponding  to/?TE  =  50  rad/mm.  This 

result  is  obtained  by  theoretical  consideration  of  the  relation 
between  second-harmonic-intensity  modulation  and  modal 
phase  difference. 

From  theory  we  have 


E,(z)  —  Re[e"*(A,e  ~rf  +  B,er'“ 


FIG.  I .  Experimental  setup  used.  L ,  couples  light  into  the  GaAs  waveguide. 
L -  collimates  the  light  to  the  power  meter  P  L,  and  L,  form  a  low-magmfi- 
ation  microscope.  L .  and  beam  splitter  B  are  used  to  adjust  coupling 
through  lens  Z.,  to  the  waveguide.  The-i  /2  plate  is  used  to  change  the  polar¬ 
ization  of  the  incident  fundamental  beam.  The  waveguide  under  considera¬ 
tion  is  shown  in  the  inset.  It  was  grown  by  molecular  beam  epitaxy  over 
[1001  oriented  GaAs  wafer. 


422 


Appl.  Phys.  Lett.  52  (6).  3  FeOruary  1 988 


0003-6951  /88/060422-03S01  00 


©  1988  American  institute  of  Physics 


422 


FIG.  2.  Picture  of  the  square  of  the  near-field  pattern  of  the  second  harmon¬ 
ic  ofNd:YAG  generated  across  the  waveguide.  There  are  16  periods  of  vari¬ 
ation  over  1  mm  of  waveguide  length,  giving  a  period  of  A  —  63  pm. 


where  i  =  x  and y  corresponds  to  the  Ex  and  Ey  component 
of  TE0  and  TM0  modes,  respectively  (a  is  the  direction  of 
propagation  along  the  waveguide),  E0l  is  the  incident  elec¬ 
tric  field  on  the  input  facet  of  the  waveguide,  t,  and  r,  are 
transmission  and  reflection  coefficients,  and  y,  =  at  —  i@,  is 
the  complex  propagation  constant.  The  wafer  under  consi¬ 
deration  has  [  100]  surface  orientation.  This  fact  combined 
with  the  symmetry  of  the  nonlinear  susceptibility  tensor 
of  GaAs  yields  second-harmonic  polarization  of  the  form 
p'^'ccd  EtE,x  +  E  \y).  The  first  component  of  polariza¬ 
tion  is  along  the  polarization  of  TEy,  and  the  second  term  is 
polarized  like  the  TM0  mode.  The  TMy-like  component  of 
second-harmonic  polarization  has  ^-directed  polarization 
and  hence  cannot  contribute  to  the  radiation  of  second  har¬ 
monic  coming  out  of  the  top  surface  of  the  waveguide.  How¬ 
ever,  the  TE„-like  polarization  of  the  second  harmonic  can 
radiate  out  of  the  top  surface.  If  we  keep  only  the  sum  fre¬ 
quency  term  and  using  Eq.  ( 1 ),  we  obtain 

p<2"'ocRe(£,£:) 


Re 


ehu'(AtB,e-2'yL) 


x(e-»'  :~L' 


&\Az  e±y.:  -  L) 
A  \B-± 


(2) 


where  y  =  Uy,  +  yz)  =  a  4-  iff.  Ay  =  y,  —  y:  =  Aa 
+  /A/?.  The  term  of  the  form  =  cannot  radiate  per¬ 
pendicular  to  z  direction  and  will  be  ignored.  Also  from  Eq. 
( 1 ).  B{A2/A,BZ  —  r,/r,.  Using  this  in  Eq.  (2).  we  obtain 


-Re 


el‘"'iA,B,e 


(3) 


The  time-averaged  second-harmonic  radiation  is  then 
<x  A,B~  1e~XAL 


X  -^a  ■+■  —  j  -1-  cos[2A/?(z  —  L  )  ] 


(4) 


where  a  =  (rl/r,)e;Aa,i  *  L'  =  1  —  A  a.  In  terms  of  Aa,  Eq. 
(4)  reduces  to 


{P:1“'),-x\AtB2fe- 


(A  a);  01  Aa  ) 


-r  cos(2A£(z  -  L) 


(5) 


For  low-loss  birefringence  (  21aL  =2  a,  —  a;Z.<  l),Aa 
becomes  I  —  rs/rv  which  is  less  than  0.2  for  usual  wave¬ 
guides.4  Therefore,  we  can  safely  ignore  Aa  in  Eq.  (5).  Fig¬ 
ure  3  shows  the  theoretical  variation  of  (P':“'),  as  a  func¬ 
tion  of  length  across  the  waveguide.  For  the  purpose  of 
calculation,  the  actual  waveguide  is  approximated  by  a  slab 
waveguide  shown  in  the  inset  of  Fig.  3.  As  can  be  seen,  the 
theory  agrees  well  with  the  experimentally  observed  second- 
harmonic  variation  of  Fig.  2.  The  most  uncertain  parameter 
of  the  waveguide  is  the  aluminum  mole  fraction  of  the  clad¬ 
ding  layer,  which  was  approximately  20%.  Then  according 
to  Eq.  ( 5 ) ,  the  birefringence  at  A  =  1 .0642 /im  is  given  by 

A£  =  /?TE.  —/?TH  =  77-/A  =  50  rad/mm,  (6) 

where  A  —  63  fim  is  the  period  of  the  variation  of  the  second 
harmonic  signal  observed.  The  experimentally  observed 
A  —  63  jj. m  agrees  well  with  the  theoretical  prediction  of 
A  — 60  /im  (Fig.  3).  There  are  several  points  that  need 
further  discussion.  For  small  Aa  (that  is  r,/r: ~  1 ),  the  sec¬ 
ond-harmonic  polarization  has  a  standing-wave  profile, 
varying  sinusoidally  in  z  along  the  direction  of  the  wave¬ 
guide.  This  standing-wave  polarization  acts  as  a  source  for 
the  radiation  of  second-harmonic  signal.  Assuming  that  the 
waveguide  is  charge  free,  the  polarization  profile  generates 
an  electric  field  along  its  direction,  E  =  —  4-P.  Since  the 
lens  of  the  camera  in  our  experiment  maps  the  near-field 
pattern  of  the  electric  field  across  the  waveguide  to  the  image 
plane  on  the  film  negative,  the  resulting  picture  of  Fig.  2 
registers  the  time-averaged  square  of  the  near-field  pattern 
or  the  time-averaged  square  of  polarization  profile  across  the 
waveguide.  This  point  shows  the  relation  of  Eq.  ( 5 )  to  Fig.  2. 

Another  point  to  consider  is  the  coherence  of  the  second 
harmonic  signal  (SHS).  The  far-field  pattern  approxi¬ 
mately  proportional  to  the  Fourier  transform  of  the  near- 
field  pattern.  With  the  use  of  the  near-field  pattern  discussed 
above,  the  far-field  pattern  is  calculated  to  have  intensity 
peaks  at  angles  \6 ,  =i;„,/2A<l  rad.  where  the  angle  is 
measured  in  they-z  plane  (Fig.  1)  with  they  axis  referring  to 
9  —  0.  This  prediction  of  the  theory  can  account  for  the  fact 
that  when  the  lenses  L}  and  Lt  in  Fig.  1  were  tilted  from  the 
direction  perpendicular  to  the  sample  surface  in  the  y-z 
plane,  the  intensity  of  SHS  would  sharply  decrease.  This 
observation  confirms  the  coherent  nature  of  the  second-har¬ 
monic  signal. 


«  •  i 


FIG.  3.  Theoretical  varialion  of  along  the  waveguide  length  The 

actual  waveguide  is  approximated  by  an  equivalent  slab  waveguide  shown 
in  the  inset  of  the  figure.  For  the  parameters  used,  we  obtain  A  —  bO/im. 
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The  third  point  to  consider  is  the  phase-matching  condi¬ 
tion.  A  standing  wave  can  be  decomposed  into  a  forward  and 
a  backward  propagating  wave  with  wave  vector  K  and  —  K . 
Therefore,  the  net  wave  vector  of  the  standing  wave  is  zero. 
Furthermore,  if  we  express  the  fundamental  electric  field 
profile  [  Eq.  ( 1 )  ]  in  terms  of  a  sum  of  a  standing  wave  and  a 
propagating  wave  and  following  the  steps  leading  to  Eq.  ( 5 ). 
we  find  that  the  second-harmonic  radiation  perpendicular  to 
the  waveguide  surface  is  only  the  result  of  interaction  of  the 
standing-wave  component  of  TEo  and  the  standing-wave 
component  of  TM0  modes.  Also,  the  near-field  pattern  of  the 
SHS  is  described  by  a  standing  wave.  This  means  that  the 
sum  of  the  initial  interacting  fundamental  wave  vectors  is 
equal  to  the  sum  of  the  final  second-harmonic  wave  vectors 
( which  equals  zero ) .  It  should  also  be  noted  that  since  the 
SHS  is  not  generated  collinear  to  the  fundamental  beam,  the 
concept  of  phase  matching  is  not  applicable  in  its  conven¬ 
tional  sense.  For  noncollinear  second-harmonic  generation, 
the  condition  for  the  highest  possible  efficiency  is  the  conser¬ 
vation  of  total  wave  vectors. 

In  summary,  a  novel  second-harmonic  generation  tech¬ 


nique  is  theoretically  and  experimentally  analyzed.  The  sec¬ 
ond-harmonic  signal  obeys  wave  vector  conservation,  is  co¬ 
herent.  and  propagates  perpendicular  to  the  fundamental 
beam.  This  perpendicular  propagation  makes  this  technique 
useful  for  the  study  of  waveguide  parameters.  The  experi¬ 
ment  reported  yielded  the  modal  phase  birefringence 
A/?  =  /?TE  -  0 tv,,,  accurately,  directly,  and  unambiguous¬ 
ly  in  the  GaAs/AlGaAs  waveguide  structure  under  consi¬ 
deration.  The  knowledge  of  the  modal  phase  birefringence  is 
crucial  in  design  of  TE«TM  mode  convertors,  wavelength 
filters,  and  TE— TM  optical  modulators.  These  devices  are 
important  part  of  wavelength  multiplexing  systems  and  oth¬ 
er  future  integrated  opto-electronics  circuits. 
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terminal  device  is  its  tunable  transconductance  and  impedance  This  device  can  be  used  as  a 
NAND  gate  iogic  unit  because  it  can  be  turned  off  by  either  the  top-gate  or  the  back-gate. 

In  order  to  characterise  the  GaAs  material  and  device  in  a  non-invasive  way.  we  develop  the 
electro-optic  probing  technique.  Because  the  refractive  indices  of  GaAs  are  modified  by  electric 
fields,  we  can  obtain  the  information  about  internal  fields  in  GaAs  device  or  material  from  the 
phase  retardation  of  a  probing  beam.  The  experimental  results  in  potential  profile  probing  of  vari¬ 
ous  device  structures  are  reported.  The  problem  of  three-dimensional  internal-field  probing  is 
analyzed  and  a  simple  experiment  is  done  to  demonstrate  the  feasibility.  In  addition  to  probing 
low  frequency  characteristics,  microwave  signal  can  also  be  probed  if  the  probing  beam  is  replaced 
by  a  short  optical  pulse  train.  The  standing  wave  pattern  along  a  coplanar  waveguide  is  measured 
at  fiGllz  by  this  technique.  The  measured  effective  index  in  microwave  frequency  is  very  con¬ 
sistent  with  the  value  measured  by  network  analyzer.  However,  the  detection  bandwidth  of 
electro-optic  probing  can  be  over  100GHz  if  subpicosecond  pulsewidlh  and  low  liming  jitter  light 


source  is  used. 
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